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Summary. — The dependence of the branching ratio 
Ru>e + y)/R(u > e+ y+ v) 


upon the ratio of cut-off to boson mass and upon a possible boson anom- 
alous moment is discussed. Some of the properties of the boson are 
inferred from the present experimental data. 


1. — Introduction. 


It has been suggested that the universal V-A Fermi interaction, which 
has proven so successful in correlating the experimental facts of 8 and u decay, 
arises via an intermediate charged vector toson which is coupled to fermion 
pairs. Such a mechanism would not lead to a Fermi interaction involving 
either four neutral fermions or four charged fermions, and consequently the 
decays u — 3e and K* > ++ 2v would be suppressed, in agreement with the 
experimental observations. 


(*) Work supported in part by the U.S. Atomic Energy Commission and in part 
by the Alfred P. Sloan Foundation, Ine. 
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Whereas the decay u >e+y, which has never been observed experimen- 
tally, cannot arise from the usual direct four-fermion interaction, it would 
be possible under the intermediate boson hypothesis if the boson had a non- 
infinite mass. In fact, it has been asserted that the calculated branching ratio 


o=R(pre+y)/R(u>e+v+Y) 


turns out to be much too large (1), for an experimental upper limit on the 
branching ratio has been set by LOKANATHAN and STEINBERGER at 2-10”, 
by Davis, RoBERTS and Zipr at 105, and by BERLEY, LEE and BARDON at 
210800) ; 

The calculation of the u +e + decay rate is not unambiguous, however, 
for the process involves divergent integrals, and hence any calculation requires 
that a high momentum cut-off be supplied. Although the divergence is only 
logarithmic, and hence for cut-offs large compared to the boson mass the 
branching ratio must be a slowly varying function of the cut-off, this need 
not be true for cut-offs of the order of magnitude of the boson mass or less. 
Indeed, it is not to be expected, for the principal contribution to the u > e+7 
decay comes from radiation by the boson itself, a process which depends 
strongly upon the boson momentum. The major contribution to the decay 
rate comes from large momenta of the intermediate boson, and hence the 
process is sensitive to the magnitude of the cut-off to mass ratio. In particular, 
as the mass tends to infinity with the cut-off fixed, the branching ratio 0 goes 
to zero in agreement with the direct four-fermion interaction hypothesis. 

In the present article the dependence of the calculated branching ratio 
on the cut-off to mass ratio is studied in detail. If one assumes that the 
intermediate Eoson is quite massive, e.g. if it has a mass equal to five times 
the cut-off, then the calculated branching ratio is below 2-10-6 for all values 
of the koson anomalous moment which are of the order of magnitude of a 
boson magneton. 

A peculiar observation is that if the Loson should happen to have an anom- 
alous moment of approximately 0.7 boson magnetons, then for any boson 
mass greater than about one half the cut-off, the calculated branching ratio 
would be less than 2:10, This might help explain why u — e + y decay 
is not observed, and thus a study of the possible sources of a boson anomalous 
moment is clearly suggested. 


() G. FEINBERG: Phys. Rev., 110, 1482 (1958). See also L. MicnEeL: Proc. Intern. 
Conference on High Energy Physics 253 (1958), M. GELL-MANN: Gatlinburg Conference 
(1958). / 

(?) S. LOKANATHAN and J. STEINBERGER: Phys. Rev., 98, 240 (A) (1955); H. 
F. Davis, A. Roperts and T. F. Zrpr: Phys. Rev. Lett., 2, 211 (1959); D. BERLEY, 
J. Lee and M. BARDON: Phys. Rev. Lett., 2, 357 (1959). 
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2. — Branching ratio 9 = R(u+>e-+ y)/R(u—-e+v+y). 


The evaluation of the transition rate ÆR(u —e +7) proceeds from the 
Lagrangian density Y= %+ Li + £', where 


e.m, 
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Here y,, y, Yu, A, and 9, are the neutrino, electron, muon, electromagnetic 
and vector boson field operators, respectively, F,, = 0A,/0x,—0A,/dx, and y 
is the boson anomalous magnetic moment (in units of the vector boson mag- 
neton). The only unfamiliar aspect of the vector boson field is in the free-field 
commutation relations, 


C2 
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OT oo Nes a es rar 


which are required by the demand for compatibility with the equations of 
motion derived from the Lagrangian density %. 

In the lowest order of perturbation theory three Feynman graphs con- 
tribute to the matrix element for u — e + Y. 


(1) (2) (3) 


y È e cs ae 


Letting q, g and © be the 4-momenta of the muon, electron and photon 
respectively and ¢ be the polarization vector of the photon, the S-matrix ele- 
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ments corresponding to the diagrams (1), (2) and (3) are: 
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Although individually these contributions are not gauge invariant, the gauge 
invariance of the sum is reflected in the appearance of a factor 
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in the final result. Assuming that m,<m, and (m,/M)?<1 the S-matrix 
element may be found by expanding the denominators in STES Fan SRI 
powers of q and q’. The result is S=1—i9776“(q—q'— ©), where 
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The transition rate R is (27)-!|977|2x density of states, and thus 
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which is to be compared with the corresponding rate for the normal decay 
mode of the muon 


Rfu—e+v+v) = 2 xt(g/M})tm,. 
Consequently the branching ratio is given by 
o (ox/e7) N?: 


Since the integral J, diverges logarithmically, N converges only if the anom- 
alous moment w = 1, in which case one obtains N = +1, corresponding to a 
branching ratio 0 = 0.870-10-3, far above the experimental upper limit. 
However, if u #1 a cut-off factor such as A?/(p?+ A?) must be introduced 
into each of the integrals 7, in order to obtain a finite branching ratio. Then 
it can be shown that 


T, (a) = J,,(a) Cs Ini (a) si 
where a = (A2/M?) and 
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Now in general the branching ratio can be expressed in the form 


where 
Qo = 0.890-10-*[Z, — I + 3I;PP 


is the branching ratio in the absence of any anomalous moment, and where 


is the value of the anomalous moment for which the branching ratio must 
necessarily vanish. 

It is not difficult to determine the variation of 9, and #, with the cut-off 
parameter a = (A2/M2?), for it can be easily verified that J,(a)= (a In a)/(a — 1), 
and J,(a) for n> 0 can be obtained by using the recurrence relation 
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The case a — 1 is trivial, and we have J,(1)= 1/(1+n), from which it fol- 
lows that 0, = 2.95-10-* and y, = 0.7. Thus, whereas the branching ratio is 
abnormally high for the normal magnetic moment, it becomes zero at u = 0.7. 
Although one might suspect that the suppression of the p—e-y decay for 
u = 0.7 is peculiar to this particular choice of cut-off, À = M, the constant / 
is actually very insensitive to the choice of cut-off. 

First consider the case when a= (A?2/M?)<1. In this limit J(a) & —alna, 
from which it can be shown that I,(a) & a/(1 +7). 

Therefore, 

ADE and Lo = 0.10 


whenever a < 1, and y, is not greatly different from its value at a = 1. 

In the opposite limit, 

Gi wit Mays I sa) = 031021078) nd Hola) = 1 — Ware 1: 
For large values of a, 0, diverges and &, approaches unity, as they must. 

A number of individual points other than a — 1 have been calculated, 
and the results are summarized in the accompanying graphs. The first thing 
that should be noted is the slow variation of 4 as a function of A/M. It 
follows from this that if the boson should happen to have an anomalous mo- 
ment u = 0.7 magnetons, the branching ratio would be small for a large range 
of values of A/M; specifically for any value of the boson mass M greater than 
one-half the cut-off A, the branching ratio o would be less than 2-10-. 

Although the case y = 0.7 is clearly the most favorable for the achieve- 
ment of a small branching ratio, the rate for u el can be reduced to 
2-10 in general by letting the boson mass M exceed the cut-off A by a 
factor 5 or more. This hypothesis does not contradict the weak coupling ap- 
proximation, for since (g/M?) & 10-/M?, where M, is the nucleon mass, 
M ~5A corresponds to a value of g considerably less than unity for any 
cut-off À < 10M,. By attributing to the intermediate charged boson a mass 
much larger than the cut-off, in effect it is prevented from emitting efficiently 
photons of energy m,,c?/2. 

Therefore, the conclusion is that one cannot rule out the possibility that 
the decay 1 >e+y goes via a charged vector koson of very high mass; alter- 
natively, if the boson should happen to have an anomalous moment near 0.7 4 
its mass could be relatively small without contradicting the experimental data. 


3. — Discussion. 


It is apparent that the structure of the matrix element for p>e+y¥ 
decay is the same as that of the matrix element for the boson contribution 
to the muon or electron anomalons magnetic moment. Indeed, the only modi- 
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fication necessary in order to obtain the muon anomalous magnetic moment, 
for example, is to replace the electron in the final state by a y-meson. This 
will change the numerical value of N, but it will not modify the general state- 


ments which have been made regarding 
the dependence of the matrix element 
upon the ratio of cut-off to boson mass. 
In particular, since the quadratically 
divergent part of S, is removed by a 
charge renormalization (coming from S, 
and S;), the anomalous magnetic mo. 
ment also is logarithmically divergent. 
Therefore, the selection of a cut-off 
necessary to give a small pu ->e+7y 
branching ratio will automatically gua- 
rantee a small contribution to the muon 
or electron anomalous magnetic mo- 
ment; in fact, this contribution is about 
10-° muon magnetons for the u-meson 
and hence is undetectable. 

However, a different situation oc- 
curs in the case of the boson contri- 
bution to the self-energy of the muon 
or electron. 


10 


Fig. 1. — Branching ratio o for inter- 
mediate boson with normal magnetic 
moment. 
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Here the divergence is quadratic, and hence there is the possibility of 
making this self-energy large and still keeping the branching ratio and hence 


the leptonic anomalous magnetic moments small. 
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In fact, as the preceding 
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Fig. 2. — Cut-off dependence of the boson anomalous moment which must be postulated 
in order to obtain zero branching ratio o. 


section indicates, the Boson-lepton coupling may be made strong, with an 
associated self-energy of the order of magnitude of the muon mass, without 
contradicting any present experimental results. 


7 


1 


180 M. E. EBEL and F. J. ERNST 


It must be emphasized that this is not to be taken as suggesting that the 
u—e mass difference could be accounted for by the interaction with the 
postulated intermediate Boson. Indeed, the present discussion has treated the 
muon and the electron symmetrically, and hence cannot give any mass dif- 
ference. However, it does seem desirable to demonstrate that the different 
momentum dependences of the divergences in these various processes make 
it necessary to use extreme care in inferring the strength of the self-energy 
process from the 1 —e-+y branching ratio. 


The authors are indebted to G. FEINBERG for a discussion of his work, 
and to R. G. SAcHs and H. W. Lewis for many helpful suggestions. 


RIASSUNTO (*) 


Si discute la dipendenza del rapporto di branching 
Ru —=e + y)/R(u =e + v +5) 
dal rapporto tra cut-off e massa del bosone e da un possibile momento anomalo del 


bosone. Dai dati sperimentali attualmente disponibili si deducono aleune proprietà 
del bosone. 


*\ 7 + s i 
(©) Traduzione a cura della Redazione. 
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I. - Binding Energies. 
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S. LIMENTANI, P. E. ScHLEIN and P. H. STEINEERG (**) 
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(riceyuto 11 23 Giugno 1959) 


Summary. — The analysis of 134 uniquely identified mesic decays yields 
increased accuracy in the knowledge of the binding energies of the hyper- 
DUC Hines Lin Hen, He Uninc, Sul, ands ’Bex. Ine additron, 
individual examples of the new species 7He,, By, and !2B, are described. 
The present data are combined with those collected in the EFINS survey. 
The isotopic spin multiplet structure of the light hypernuclei is discussed 
with reference to the information derived from the binding energies. 


1. — Introduction. 


Emulsion groups at the University of Chicago (EFINS) and Northwestern 
University (NU) are collaborating in the study of 7--mesic decays of hyper- 
nuclei produced by K~-capture. The present report covers the results on 
hypernuclear binding energies and is based on 265 mesic decays found in a 
first area scan ot emulsion stacks containing about 4-104 K--stars. Preliminary 
results of this collaboration were reported elsewhere (1). 

(*) Research supported by the U.S. Air Force Office of Scientific Research, 
Contract no. AF 49(638)-209. 

(**) Research supported by the National Science Foundation and by the U. 
Atomic Energy Commission through Argonne National Laboratory subcontract. 

(1) R. Ammar, R. Levi SETTI, S. LIMENTANI, P. E. ScHLEIN, W. E. SLATER and 
P. H. SreNBerG: Annual Intern. Conf. on High Energy Physics at CERN (1958); 
R. H. Datirz and M. KapLon: Reports, pp. 181, 200. 
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2. — Experimental details. 


2°1. Exposure and calibration of the stacks. — Two stacks were exposed to 
the 300 MeV/e enriched K~ beam (ratio of minimum tracks to K7, 800:1) 
at the Berkeley Bevatron. The EFINS stack consisted of 59 K-5 pellicles 
of dimensions 10 em x 15 em x 0.12 em while the NU stack consisted of 59 K-5 
and 71 L-4 pellicles both of dimensions 10 em x 15 em x 0.06 cm. The 0.12 em 
pellicles were processed unmounted, while the 0.06 cm pellicles were processed 
mounted on glass. The mean thickness of each pellicle before processing was 
obtained from measurements in five to nine different regions. The lateral 
shrinkage of each unmounted pellicle was also determined. Details of the 
methods used have been described elsewhere (?). 

A further stack, already processed and calibrated, was made available by 
the University of California Radiation Laboratory (*). An average shrinkage 
factor was used for this stack, which yielded about 6% of our events. 

The densities of the stacks were determined in two ways: a) for the EFINS 
stack, directly; b) for both the EFINS and the NU stacks, indirectly from 
proton ranges in Y* decays (65 proton ranges in the EFINS stack, 33 and 
26 proton ranges in the K-5 and L-4 parts of the NU stack). The mean range 
for each stack was compared with GILES’ (*) value for emulsion of standard 
density (3.815 g/em*); the densities of our stacks were then inferred following 
BARKAS et al. (4). For the EFINS stack the indirect density value so obtained 
was in good agreement with the direct measurement. The corrections in 
reducing the ranges used in the present work to standard density amount to 
about 1%, with a statistical uncertainty of + 0.59. 


22. Analysis of events. — For the purpose of analysis, the events may be 
divided into three classes 


a) two-body decays of the type: (4, Z) >t-+(Z+1, A) (x-r events). 


b) three-body decays of the type: (2, A) >x-+p4+(Z, A—1) (x-p-r 
events). 


c) other decays, including those in which neutrons are emitted. 


(2) W. E. SLATER: Suppl. Nuovo Cimento, 10, 1 (1958). 


(*) We are greatly indebted to Dr. W. H. BARKAS for the loan of this stack 
(Berkeley 1 U). 


(5) P. C. Gites: UCRL Report no. 8007 (1957). 


(1) W. H. BARKAS, P. H. BARRETT, P. CùER, H. Heckman, F. M. Suit and 
H. K. TicHo: Nuovo Cimento, 8, 185 (1958) 
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a) According to the ranges of the decay products, the z-r decays fall 
into three main groups, as can be seen from Fig. 1. Among these, the PELS 
and 4H, groups are readily identifiable, while the third group will be discussed 
later (Section 3°1). 
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Fig. 1. — Pion and recoil ranges for x-r events. 


b) The identities of the recoils in x-p-r decays are assigned on the basis 
of a comparison between their observed ranges and the ranges inferred from 
the vector sum P_, of the x and p momenta. The procedure given by SLATER (?) 
was used; Fig. 2 is a plot of the range-momentum points obtained for events 
in which the recoil direction is collinear with P_, within errors, and corresponds 
to Slater’s Fig. 14. The curves drawn in Fig. 2 are those calculated by 
WILKINS (5), the curve for *He being displaced by 0.5 um toward lower ranges 
to fit the calibration point from *H,—7--+‘He. The curve for *He was 
derived from the adjusted ‘He curve. This yields a very satisfactory fit to 
the ‘He points from *He, decay over the range from 3 to 30 um. 

For ranges below 3 um it is in general impossible to identify recoils uni- 
quely. The separation between the range-momentum curves for different iso- 
topes becomes at this range, even for flat tracks, smaller than the range measu- 
rement errors. Events with slightly larger recoils but large dip angles can 
also in some cases not be uniquely identified. (See Appendix I for further 
discussion.) 


(5) J, J. WitKIns: A.E.R.E. G/R 664 (1951). 
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1807 (2 Body decay 
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F © Unidentified 
2 Body decay 
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10 20 40 60 80 100: 
Ro (uw) 
Fig. 2. — Range momentum curves for recoils from hyperfragment disintegration. 


c) Events in this class, representing a small fraction of our data (~15°%) 
were analyzed with an electronic computer programmed (*") for the permis- 
sible range of Z and A values for each prong. For every permutation of these 
(7, A) assignments the computer printed out the A-binding energy (B,) and 
the momentum unbalance (AP) for events yielding —30 < By, < 40 MeV. 
An event is considered unique if there exists only one set of prong identities 
for which AP is zero within experimental error (*). This experimental error 
was calculated for each event from the assigned uncertainties in the measured 
input quantities (ranges, angles, etc.) by the electronic computer. For events 
assigned to a well-known hypernuclide, we require in addition that the cal- 
culated B, be consistent with the known average B,. 

Non-colinear two-prong and non-coplanar three-prong decays which require 
at least one neutron are considered to be unique only if the analysis with a 


(9) E. M. SILVERSTEIN: Suppl. Nuovo Cimento, 10, 41 (1958). 
(7) F. W. Inman: UCRL Report no. 3815 (1957) 
() The spread of AP which may reasonably be expected can be seen from Sect. 3°2 
where the AP-distributions for unique “HA and 5He, are disc 


ussed. 
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single neutron yields a single positive binding energy. The identity of four- 
prong decays, where the emission of neutrons cannot be ruled out on the basis 
of coplanarity, will always be open to question. No attempt was made to estab- 
lish any decay mode involving more than one neutron. 

Finally, whenever possible, the charge of the hyperfragment was deter- 
mined from profile (8) measurements on the connecting track. 

For several events in which the kinematics yielded only two alternatives 
differing in Z, it was possible to establish a unique identity with the aid of 
such measurements. Details of the method are given in Appendix IT. 


2°3. Conventions adopted. — The binding energy of the A in a hypernucleus, 
B,, is defined as the energy required to remove the A, leaving the residual core 
nucleus in its ground state. 

B,(Z,A)=%—Q, where @ is the total energy release and Q, is the dif- 
ference between the sum of the A and core nucleus rest energies, and the 
sum of the decay product rest energies. 

In the case oi °Be,, where no stable core exists, the ground state of *Be 
(unstable against ‘He emission by 96 keV) is taken as the residual core. 

The standards adopted for the determination of binding energies are the 
same as those described in ref. (7), except for Q,, the energy release in the 
r_-decay mode of the free A. 

The following values of Q,, all obtained from K -induced events, 
have been averaged: ((37.45 + 0.17) MeV) (°), ((37.58 + 0.18) MeV) (©) and 
(37.61 0.16) MeV) (11). The adopted value is then (37.58 + 0.15) MeV. The 
error assigned to this mean is considered more realistic than the one obtained 
from combining the three independent error estimates. Preliminary results 
of a new determination of Y, at Berkeley (1?) are consistent with the one 


adopted here. 


2°4. Hrrors. — With the exception of the uncertainties in nuclear rest ener- 
gies (of the order of 0.02 MeV), the error in @, equals that in Q,. An esti- 


(8) G. Anyrat, A. Bonerrr, C. Dicwortx, M. Lapu, J. Morgan and G. OccHIA- 
Lint: Suppl. Nuovo Cimento, 4, 244 (1956). 

(9) W. H. BARKAs, P. C. Gires, H. H. Heckman, F. W. Inman, C. J. Mason and 
F. M. Smirn: Padua-Venice Conference Report (1957), p. vr-41. 

(19) J. Bogpanowicz, M. Danysz, A. Firrpkowsky, E. MARQUIT, E. SKRZYPCZAK, 
A. WroBLEWSKI and J. ZAKRZEWSKI: Nuovo Cimento, 11, 727 (1959). 

(1) B. Brucker, F. ANDERSON, J. LoDGe and A. PEVSNER: kindly communicated 
to us by A. PEVSNER (1959). 

(12) Kindly communicated to us by W. H. BARKAS. 
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mate of the errors in 153 (LEA OO can be obtained from the following 


compilation: 


| E | 

oa tay | : in | 

Saree Ba dA | È “i q@ Correlated | 

or pan 
| (a) R-E relation 0.10--0.20 | 0.10+0.20 | yes 

| (b) Stopping power | 0.05 —0.10 0.05 0.10 | no 
| (c) Emulsion shrinkage ~ 0.05 | ~ 0.05 | no | 
| (d) Random errors ~ 0.08 0.08 —0.60 | no | 


Since the same range-energy relation is used in determining both Q, and a) 
the error propagated to B, is negligible for most hypernuclear species. The 
entire set of Bes could be shifted by a concomitance of errors b), ¢) and d) 
in Q, and of errors b) and c) on Q by a total of about 0.2 MeV. 

For events not involving a neutron, the random error in B, is derived 
from the standard deviation o of the experimental distributions for uniquely 
identified decays: these are dis- 
played in Fig. 3, where B,-histo- 
grams for the most abundant spe- 


T T “at T TT 
N Efins Co Mu 


3H, (6) 5 cies are plotted. The quoted errer 
Lo im B, for a given hypernuclide, 
4H, (28) | Oy, is 0/y®, where n is the nurn- 


ber of events. The variation of 


di SAN Mai 
“He(14) o with @ is shown in Fig. 4. For 
t the non-abundant species where 


5 it is not possible to obtain o in 
— this way reliably, a value appro- 


5 . 
He,(60) | priate to the energy release was 
id inferred from Fig. 4. In those 
cases where a neutron was emit- 
a : 
7Li, 12) j ted, the B, error was inferred 


from a knowledge of the straggling 


and measurement errors. 
Sources which contribute to 

the spread of the B,-distribution 

include range straggling, measu- 


Non uniquely identified (79) 


È (E SEAN I n RA rement errors, distortion of the 
B, (MeV) emulsion and uncertainties in the 

Fig. 3. — Binding energy histograms for the Shrinkage factors (?). Our overall 
most abundant hypernuclear species. experimental error from these sour- 
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ces probably does not exceed 1%. In fact the range stragglings obtained com- 
pare with the theoretical one as follows: protons from X*-decays, (1.9-+ 0.2) % 
vs. 1.4%; pions from x-r decays of ‘H,, (2.7 + 0.5)% vs. 2.7%. 


16 | | 


è 81; 
14 — A —|| 


0.4 I | | 
30 40 50 60 


Q (MeV) 


Fig. 4. — Variation of o, experimental spread of B, distributions, with the energy 
release, ©. 


3. — Results. 


8 1. Uniquely identified events. — The binding energies presented here were 
all deduced from « uniquely » identified events in which all decay tracks stopped 
in the emulsion (134 out of 214 complete events). With the exception of one 
decay in flight of *H, (6, not included in the average), all the decays occurred 
at rest. 

Although the B, value of an individual event in no case formed the basis 
for identification, those events with B, departing by more than 30 from the 
mean, (B x), were rejected; this was only done for the species represented by 
more than ten events. One ‘H, with B, — 6.1 MeV and decaying by the 
mode 7r:-+‘He was rejected on this basis. 

Table I gives a breakdown of the data used in the binding energy deter- 
minations. Since only complete events are considered, the frequencies do not 
represent the true relative abundances. 

There appears to exist a small systematic difference (of order 0.6% @) 
between the EFINS and NU results, as can be seen from the data for ‘H,, 
‘He, and *He,. This difference however is of the order of magnitude of the 
uncertainties in the stopping power calibrations, although it may arise from 
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Taeze I. Mesic hyperfragments from K~-capture decay modes of uniquely identified events: 


| EFINS | NU | 
| A i € 
| Decay mode Di PRIA | a Total 
| (MeV) | events | (MeV) events 
| ce = na DI — 
| 
PEN re le 0.06 | 3 — 1.28 1 4 
+n +1H+?H 0.05 il _ — il 
ne it 2 0.79 1 — | — | I 
| 4H, +x-+‘He 2.58 12 EN I - 
| —7--+n-+%He 1.61 2 — — 2 
7 +1H+$H 2.05 4 0.66 I 5 
>n+2?*H 2.21 1 — = ih ()} 
‘He\>x-+!H+*He 2.44 8 Dallo 5 13 | 
>r+2!H+?H 3.09 I _ — 1 | 
5He, +n +1H+*He Dols 35 2.95 25 60 
"He, > +°H +*He 3.04 il — — 126) 
"Li, 7 +"Be (prob. 6.02 | 8 6.57 1 9 
7 +1H+61i 5.68 il | 5.44 1 D 
—>n-+%He+4He 5.47 1 a == 1 
SLi, nm +2%He 5.97 2 6.59 3 | 5 
"Li, +n +n+24He 7.95 2 | 6.41 | 1 | ts 
Be 7 +1H+24He 6.07 1 6.48 | 2 es | 
UB, > +%He+24He 9.93 1 — | — i) | 
BB, 7 +34He — — 9.55 | 1 eg | 


(*) Heretofore unobserved. 


other sources. The averages of the results from both laboratories are, in our 
opinion, the best estimates of the binding energies B je These Ba are given 
in Table IT, where the PE from the EFINS survey (#) are included for com- 
parison; good consistency between both sets is apparent. The final column 
of Table II gives the grand averages. 

A preliminary report of this experiment (!) quoted a negative value of B in 
for °H,, suggesting that the value of Q, adopted at that time, that is 
(37.22 + 0.20) MeV, might be too low. The Q, value here adopted on the 
basis of new evidence, together with an increased accuracy in the stopping 
power calibration of the present stacks, now yields a By for *H, which is 
consistent with a small positive value. 


(13) R. Levi SETTI, W. E. SLATER and V. L. TELEGDI: Suppl. Nuovo Cimento, 10, 
68 (1958). 
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TABLE II. — Binding energies from uniquely identified mesic decays. 
The value @, = (37.58-+0.15) MeV is adopted in the present compilation. See text 
for discussion and references. 


EFINS survey EFINS-NU — Average of survey 


May 1958 collaboration and collaboration 
Identity 
Ba gdo, (0) No-rofl By | 0, (A) ING, of] By | Gan (NoN 
(MeV) | (MeV) levents| (MeV) | (MeV) levents| (MeV) | (MeV) events 
DEUX 0.56 0.50 7 |—0.06 | 0.30 6 0.12 0.26 13 
4H, Dele 0.20 21 PB || OS) DT 2:20) 0.14 48 
‘He, 2.35 0.20 9 223020 14 | 2.36 0.12 23 
‘He, 3.18 0.20 17 3.06 | 0.10 60 | 3.08 0.09 Ti 
Hex — — — 3.0 0.7 i | Bo) Onl 1 
"Li, (0) 2-bodyf | 4.8 0.8 1 6.08 | 0.33 95:96 NO ORAZIO 
other || 5.4 0.6 2 5.51 0.40 3 | 5.46 0.33 5 
ID 6.0 0.4 3 6.34 | 0.67 | ini 0.35 8 
Li, Tale Ose 1 TANTO SRO 4 
Bey 6.6 | 0.6 1 = — | — | 66 | 0.6 1 
Be 6.8 0.4 3 6.4 055 3 | 6.60 0.31 6 
EBS | 9.9 0.6 I 9.9 0.6 1 
DA | 9.6 0.6 1 9.6 | 0.6 1 
(a) The systematic errors (~ + 0.2 MeV) have not been included. 
(6) Reasons for separating 2-body and other decay modes are discussed in Sect. 8'1. 


Some remarks are required concerning those 9 x-r events which are con- 
sistent with the decay scheme 


(1) oi Tee ber 


Because of the large relative uncertainty in determining the range of each 
recoil, it is not possible to infer the identity of individual recoils through a 
range-momentum relation. Both the pion and recoil show range spectra (see 
Fig. 1) which are not in consistent with the decay of a single hypernucleus. 
An average recoil range of (1.85 + 0.2) um was found while the average pion 
momentum is (1070.5) MeV/c. If all these events are to be attributed to 
the decay of a single species, then this would rule out an identification of the 
recoil as *Li (from the decay of °He,). These events could also be alternatively 
interpreted as the decays of certain hypernuclei with Z > 3, of which the lightest 
one is Be, . Lacking reliable range momentum curves near this recoil range, 
one can not rule out this alternative. However, as the mesic decay mode should 
be greatly suppressed with increasing Z, and as there is no indication that 
“Be, or any of the heavier hypernuclides be singularly abundant, we believe 


FZ 13 - Il Nuovo Cimento. 
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that the more plausible interpretation is to attribute these events to the decay 
Ole Lik: 

Among the published events with the same decay characteristics as (1), 
one (14) was identified as “Li, from an analysis of the kinematics both of 
production and decay. 

One can, of course, nevertheless not exclude the possibility that among 
these events there may be some which represent the decay of ‘He,, "Be, 
and possibly even heavier hypernuclides, nor is it possible to determine 
whether there is a mixture of two competing decays of “Li,. Thus the 
recoil 7Be could e.g. be left in the first excited state (at 0.43 MeV above 
the ground state). For these reasons the B x Of ‘Li, from these two-body 
decays is indicated separately in Table IT, although it is not inconsistent with 
the value obtained from the other decay modes. 

As indicated in Table I, the decay modes: 


(2) Wey aye oak AE CE 
(3) lp eee on He ee oe 
(4) eB cnr Ee, 


have been identified for the first time. The limitations in establishing the 
existence of a given species or decay mode on the basis of one isolated event 
are self-evident. The reliability of these events may be judged with reference 
to Appendix III where the experimental details are presented. Appendix ITI 
also contains the data for one event representing the decay mode 


3 
OL 
Za 


‘Li, —> 7 +4 *He+*He 
as well as three new examples of the decay mode. 
(6) "Li, «+4 n+ 2 4He 


confirming previous observations (11), 


32. Non-uniquely identified events. — For three body events in this class, 
the binding energies were calculated assuming +-p-r decay modes. The charges 
of the recoils could often be determined by extrapolation even where the range- 


(4) F. C. GILBERT, C. E. VIOLET and R. 8. WHITE: Phys. Rev., 103, 248 (1956). 

(5) P. H. FowLER: Phil. Mag., 3, 1460 (1958) 

(16) S. LIMENTANI, J. H. ROBERTS, 
Cimento, 9, 1046 (1958). 


P. E. ScHLEIN and P. H.F StEINBERG: Nuovo 
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momentum curves do not enable one to discriminate between the isotopes, 
as can be seen from Fig. 2. An average of the mass, compatible with the 
inferred charge, was used in calculating B,. In most events of this category 
the recoil energy is sufficiently small so that B, was insensitive to the choice 
of recoil mass. The B,’s so obtained are included in Fig. 3, where it is seen 
that a large fraction of the events peak at B, values corresponding to A = 4 
and 5 hypernuclides. In addition, a few events with very low B, suggest 
the presence of *H,, and events with higher binding energies imply a sub- 
stantial contribution from hypernuclides of 7 > 2. Neither in this group of 
non-unique events nor among the incomplete events is there any evidence 
for energy releases inconsistent with the mesic decay of a A-hypernucleus, 
namely larger than about 57 MeV. 


33. Identification bias. — The criteria used in assigning a unique identity 
discriminate against events with very short recoils, as well as against those 
involving neutron emission. For example, decay modes such as 


(7) a) *H, x +'H+'H, 


(8) b) (HO, tk-4He-<Ho,, 


if they exist, are likely to remain unidentified. 


a) There are 8 events without a visible recoil, in which the decay simu- 
lates the decay of a A at rest. In view of their very low B,’s, these events 
are attributed to the decay of *H,, although they could be due to *H,, 

Two of these events had sufficiently long hyperfragment tracks to permit 
mass determinations from gap length measurements. These yielded values of 
(3.1 + 0.37)m, and (3.25 + 0.64)m,, indicating very strongly that these are 
5H, decays. Nevertheless, we do not use these two events for a B, deter- 
mination since one cannot decide which of the decay modes 7-+!H+?H. or 
a +2'H +n is involved without recourse to the a priori knowledge of B, for °H, . 


b) No clear-cut examples of *He, have been detected. However, the 
decay mode 


(9) SH edera pre ie 


may readily be confused with the abundant decay 


(10) "He, re AD RUE 0; 


if a slow neutron is emitted. Those decays in which the neutron carried off 
an appreciable momentum could be detected by a lack of coplanarity and 
large momentum unbalance. 


tn] 
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A plot of the residual momentum AP ir 


om the unique ‘He, events is 


shown in Fig. 5 along with a similar plot for the decays ‘H,> am +4He, 


O (a) es mete | 
BE 4 
0 (es 


tt 


Number of events 


me 
(b) SHe, = m- +H #He 
20 
el 77 [=] étins 
Yj ne 
10 YY 
5— 7 
0 Ver LR eee 
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O10 20° 305 co 
AP(MeV/c ) 


Fig. 5.— Histograms of the residual momentum, 
# AME sires 
(a) *H,—>7-+4He; (b) ‘He, +z-+1H+*He. 


in which the emission of neu- 
trons can be ruled out. A com- 
parison of the two indicates 
that it is very difficult to de- 
tect the emission of neutrons 
with momenta < 40 MeV/c 
against the existing «noise 
level» due to a combination 
of instrumental and measure- 
ment errors. It is clear that 
the AP distribution for ‘He, 
decays may be well understood 
without recourse to the hypo- 
thesis of neutron emission. 

To determine whether slow 
neutrons (< 40 MeV/c) could 
indeed have been emitted, it is 
useful to set limits on the 
energy release in the hypothetic 
“He, decay. The condition for 
the stability of “He, against 
decay into ‘He, +n is that Q, 
be less than Q;. Similarly the 
stability of He, against decay 
into *He,-+n, would impose 


a lower bound on Q,: Using the data of Table II these conditions become: 


(11) (84.5 22 0.1) = 0 = (336-2 0:7) Mev. 


Since % includes the energy carried by the undetected neutron (< 0.8 MeV) 


the visible energy must satisfy the relation 


(12) (34.5 + 0.1) > Q, (visible) > (32.8 + 0.7) MeV. 


i (I 6 ‘ i 1 
Thus possible “He, events may be included among those *He A events having 


(13) (3.1401) < B, < (4.8 


f= 0.7) MeV . 


On the other, hand, the binding energy of *He, is 3.1 MeV, with an observed 
of 0.6 MeV. For a normal distribution about 5% of the events are expected 


N 
a 
ri 
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to yield a B, > 4.3 MeV; only one of the 60 presumed °He, events of this 
experiment has a binding in excess of this value (B, =4.6 MeV, AP=12:5 MeV/c) 

Since neither the condition on Q, nor, as seen above, detection of neutrons 
through large |AP| alone gives a sensitive enough test for *He,, a search 
was made for a correlation between |AP| and B, among *He,, with no posi- 
tive results. 

The same difficulties are encountered with non-uniquely identifiable events, 
and are here enhanced by the large uncertainties associated with the presence 
of very short recoils. 


4. — Discussion and conclusions, 


The variation of B, with mass number A is shown in Fig. 6. As previously 
noted in the literature, there is a general tendency for B, to increase with A, 
although there appear to be some departures from simple linear dependence. 
There are also indications for a structure due to the existence of several 
isotopic spin multiplets. 


T pre pr = pet) T 
np, 
12 
10 B 
Pei °Be, tor; 
8 Li, (m-r) \ 
SÙ6 TÉ a 
x 
4 “Hx we 
ry THe, BLi 8 
à SHe, + en “Lin Be, 
0 ——-$ 5H — | 
fines | ee fl Lt L Je EA] (= 


Mass number A 


Fig. 6. — Plot of B, vs. mass number, 4. 


Within the experimental errors, ‘H, and ‘He, are seen to have nearly 
equal binding energies. This suggests that they are members of an isotopic 
spin doublet, consistent with the assumption that the A has isotopic spin 
T— 0, and with the requirements of charge symmetry, as first pointed out by 
DALITZ (17). The appearance of only one hypernuclide each of mass 3 and 5 
suggests accordingly that they are isotopic spin singlets. 


(47) R. H. DALITZ: Phys. Rev., 99, 1475 (1955). 
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The multiplets structure for À — 7 hypernuclides may be understood with 
the help of a level diagram (Fig. 7), in which the ordinate represents total 
binding energy. 


24 
aes LIRE 
Ze e 
28 (T=1) I 
(T=1) Ten 4He +'H+H 

= 6He N 

3 WLLL (Li) (Hypothetical) 

= LPEEIEIE AIA Aa 

sl zz ZT i 

œ GLLLLLLLE SU (720) ET A 

È "He 

2H +°He Sal 

O 34 A 

al 

TS 

Wa) 

3 r36j= = 
Shah 7Li 
40 

Fig. 7. — Isotopic spin multiplet structure of A=7 hypernuclei. 


The ground states of the core nuclei ‘He and *Li (7=1 and T= 0 res- 
pectively) are indicated by full lines. The first 7=1 excited state of ‘Li, at 
3.6 MeV above the ground state, is indicated by a dashed line. The 7=1 
state of ‘Be is unbound, decaying into ‘He +2p with an energy release of about 
2 MeV (!8). This 7=1 state is indicated as a dashed line in Fig. 7, for com- 
parison. 

The levels corresponding to the observed "He, and "Li, binding energies 
are shown together with their experimental errors. The Bator site is 
(2.5 + 0.8) MeV less than that for "Li,, where the latter is calculated ex- 
cluding the two body decays. This suggests that they may not belong to the 
same isotopic spin multiplet. The ground state “He has T— 1 while that of 
‘Li has T—0. If the presence of the A does not greatly affect the structure 
of the core nuclei, it is reasonable to expect “He, and "Li, to have isotopic 
spins 1 and 0 respectively. This is consistent with the indications drawn from 
the binding energies. In fact, if the T— 1 state of "Li,, (*Li), existed (*), 
it would decay rapidly by heavy particle emission, even though the available 
decay modes, shown in Fig. 7, violate isotopic spin conservation. The third 


(78) F. AJZENBERG and T. Lauritsen: Rev. Mod. Phys., 27, 85 (1955). 


(*) We propose the following notation: A bound to A*— (A*)A3 excited A,:(4,)*. 
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member of the triplet, 7Be,, corresponding to the unbound ‘Be core, has as 
yet not been detected. 

The spin of 7He, (and hence of all members of the 7= 1 triplet) should 
be J(’He,) = +, if one takes J(A)=+ and J(*He)= 0, and if one assumes 
that the A is in an s-state relative to the core. The 


"Li, ground state would 
have a spin J=# if the A-nu- 


and °B, as the other members 
of the triplet. However, if (°*Be), 


cleon potential were stronger in AT (T=1) È 
the antiparallel spin orientation. =) a | 
The “Li, and *Be, binding One eBe 
energies are quite compatible with gt BL = 
an expected 7 = doublet for the 
mass 8 hypernuclei. LE 7 
Fig. 8 shows a level diagram 46K WE + 
for the A= 9 hypernuclides. To 3 LLL: 10 fe) | è} 
date, only ‘Li, and ‘Be, have Se FC (Hypothelicay 
been identified. As the ground D sot \ \ 
Statel of li basi Za; Di, È \ \ vi 
should have 7 =1, with (**Be), orth 
= 
pe) 
© \ À 
were produced it would be un- & 567 Pies ra ont 
stable with respect to breakup sg PeaHefte \ v 7 
into heavy particles, in analogy FA dn | 


with the case of (7*Li), discussed 
above. All the observed Be, 62 À 


A ue 


Fig. 8. — Isotopic spin multiplet structure of 
‘Be, 7x +p+2 2Hex A=9 hypernuclei. 


decays are of the form: 


For all events observed the relative energy of the two ‘He in their barycentric 
system indicates that they come from the break-up of the ground state of 5Be. 
This supports the assumption that in ‘Be, the A-particle is bound to this 
latter system and consequently that *Be, has T— 0. 

Taking as core nucleus two 4He nuclei with relative energy 0.096 MeV, 
SuH (1°) has calculated the A-nucleon interaction strength needed to give the 
observed B, of Bo Suh’s calculation gives results consistent with the 
A-nucleon interaction strength required for 5He, (°°). 


(19) K. 8. Sun: Phys. Rev., 111, 941 (1958). 
(®) R. H. Datirz and B. W. Downs: Phys. Rev., 111, 967 (1958). 


195 


196 R. AMMAR, R. LEVI SETTI, W. E. SLATER, S. LIMENTANI, ETC. 


We are indebted to Professors V. L. TELEDGI and R. H. DALITZ for much 
helpful assistance and criticism. We wish in particular to thank Professor 
DALITZ for discussions concerning the theoretical aspects of our problems and 
Professor TELEGDI for help with the writing of the manuscript. 

We wish to express our appreciation to Professors L. M. Brown and J. H. 
ROBERTS for their valuable assistance and for helpful criticism and discussions. 

We wish to thank Dr. F. INMAN for making available to us his IBM 650 
program for hyperfragment analysis; this was adapted for some of our needs 
by Dr. E. SILVERSTEIN; we are deeply indebted to him for this. 

At the Enrico Fermi Institute for Nuclear Studies, J. BAPTIST, T. BARTHA, 
M. McCoy, P. KLIAUGA, E. Moore, J. Morr and P. SATERBLOM; at North- 
western University, D. ABELEDO, B. BUNKER, P. KAEN, S. uan M. 
KAPLOW, L. WATSON, D. TROXEL and R. TREVINIA assisted in the scanning 
and measurements. We are grateful to all of them for their competent efforts. 

We express our thanks to Dr. E. J. LOFGREN, Dr. H. HECKMAN, and the 
Bevatron staff whose co-operation made this experiment possible. 


APPENDIX I 


Range measurements of recoil 
tracks. 


nm wep manie 
F2 n-r Event =] 

One of the limiting fac- 
tors in obtaining precise 
measurements of short tracks 


10 ~Ro=25 pm 


0 not lying in the emulsion 
1 plane is the accuracy with 
che 7 which the vertical compo- 
g nent (Az) of the range can 
‘s 20 5=Ro <10pm be determined. The mea- 
$ sured value of Az is therefore 
SH used only to establish that 
> the event is colinear or co- 
; planar within the experi- 

10 0< Ro =. Sum 
Fig. 9. — Histograms of diffe. 
v ni > È rence between inferred and mea- 
40 sured di ious 
ein ured dip angles, AA, for various 


intervals of recoil range. 
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mental errors. If this is the case, the dip of the recoil is inferred to be that 
value 2, which makes the event colinear or coplanar (2). This value, together 
with the measured horizontal component is used to calculate the recoil 
range. The difference AA between the observed and inferred dip angles is 
shown in Fig. 9 for different intervals of recoil ranges. As may be expected, the 
shorter recoils are likely to yield higher values of A. 


APPENDIX II 


Charge discrimination from profile measurements. 


For the NU stacks, consisting of 600 um pellicles, the mean track thick- 
ness { measured over a projected distance of 9.1 um centered at 31 um of residual 
range was used as a parameter for charge discrimination. A cell length of 
0.5 um was used. The overall magnification was 6500 X. Fig. 104 shows ¢ 
as a function of dip angle for the measured events. 


E : 
ae 3 60 
ra © ie) o 
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$ Bre DS Me de 
~ 05 S (o) 
= 820 © v 
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Track dip angle (degrees) 
a) NU events b) EFINS events 


Uniquely identified hyperfragments e: Z=1, 0: Z=2. Protons from uuique t-p-r decay: A. 


Fig. 10. — Charge discrimination of hyperfragments and their decay products from 
profile measurements. 


For the EFINS stack, consisting of 1200 um pellicles, it was found more 
convenient to use another parameter. The thickness t; was measured at inter- 
vals of 0.5 um over the last 55 um of track. The skewness of the distribution 
of t;, found to be a charge sensitive parameter (21) was then used. Fig. 10b 
shows the results, plotted as a function of dip, for those events which could 
be measured. 


- 


(21) R. Ammar: Suppl. Nuovo Cimento, 15, 181 (1959). 
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APPENDIX III 


Description of individual events. 
One event was observed for each of the following hypernuclear decay modes: 
Q,= 45.12 MeV (Event 2-55-2); 
Q, = 37.13 MeV (Event 3-46-3) ; 
Q, = 46.25 MeV (Event 834). 


a) "Hey >x-+*H-+‘He 


b) 1B, > +*He+2*He 


Event 2-55-2. — Table III summarizes the range and angle measurements. 
This event cannot belong to the class z~-p-r. In fact the minimum energy 
release involved [(39.8+0.7) MeV, on the assumption that prongs 2 and 3 are 
protons] exceeds QA by a significant amount. Analysis of the event on the 
Silverstein’s permutation code (*) yields only one acceptable interpretation, 
i.e. decay scheme a), with Ba = (3.0+0.7) MeV and momentum unbalance 
AP = (24.3+23) MeV/e. Any other interpretation yields negative binding 
energies. Among these, those referring to known species are smaller than 
—2 MeV and involve incompatible kinematics. Interpretations in terms of 
other new species yield Ba smaller than —1 MeV, and AP in excess of 
50 MeV/c. Decay scheme a) is therefore strongly preferred. 


Tapce III. — Event 2-55-2. He, x +°H+'He. 


tv Colatitude Azimuth Range 

8 (degrees) (degrees) (um) 

NE Gide) 25.4 64.2 PAIAIDAL 
2 (EH) SORA 216.1 91.9 
3 (He) 0 4.9 3.8 
HF (“He,) 21.6 249.4 45.6 


Event 3-46-3. — The data for this event are contained in Table IV. The 
establishment Z=2 for track 3 from thickness measurements (21) narrows 


TABLE IV. — Event 3-46-3. 1B,|>7-+3He+24He, 


Prong Colatitude Azimuth Range 
(degrees) (degrees) (um) 

1 (can) 42.1 : 

2 (He) Six: “#0 dito 1 

3 (He) 33.6 248.6 ae 

4 (4He) SI 107.5 91.5 

HF (B,) 1.6 338.1 133 
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the possible interpretations to two 


i) UB, >x-+*He+2‘He, Ba = (9.93+0.6) MeV, AP = (23+19) MeV/c; 
il) “B, >x-+n-+34He, Ba = (12.13+1.6) MeV, P, = (66+20) MeV/c; 


i) is the only interpretation not involving neutron emission for which 
AP < 50 MeV/c and ii) is the only one with neutron emission for which By > 0. 
On the basis of the low value of AP for i) we conclude that it is probably 
a more likely interpretation. 


Event 834. — It is also a four body decay and is described in Table V. In 
this case the possible interpretations reduce to: 


i) "By >x-+34He, Ba = (9.55-40.5) MeV, AP = (26+33) MeV/c; 
ii) *By>n-+n+34He, Ba = (5.9411.6) MeV, P, = (26+33) MeV/c. 


Again, because of the low value of AP, it is concluded that i) is more prob- 
ably correct. 


TABLE V. — Event 834. 1B, 7x +34He. 


Prone Colatitude Azimuth Range 
È (degrees) (degrees) (um) 
Ge) — 66.6 2.0 14511 
2 (4He) 54 290.5 9.3 
3 (He) 44 215.8 5.5 | 
4 (4He) Dore 94.2 10.6 
HF (PB) 15.8 IAE 4.9 


Finally, Table VI gives the original data for one event uniquely identified 
as 7Li, >7-+*%He+‘*He, (Q) = 41.60 MeV), while Table VII describes three 
examples of the decay °Li\,->T-+n+2 He (Q,= 52.87 MeV). 


TABLE VI. — Event 3-35-2. 7Lij,>n-+*He+4He; 
B,= (5.47+0.6) MeV, AP = (27.4+11) MeV/c. 


Colatitude Azimuth Range 

Prong (degrees) (degrees) (um) 

1 (n°) — 16.4 172.0 10300 
2 (He) 127 163.4 23.0 
3 (*He) 0: 347.8 40.1 
HF ("Li,) 81.8 30.1 39.0 
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Colatitude Azimuth Range 
Prong (degrees) (degrees) (um) 
Event 3-11-3. ‘Li, 7 +n+2*He; By = (9.0+2.38) MeV, P,,= (94.9419) MeV/c. 
Wa) — F9 188.7 20 674 
2 (He) 36.5 310.2 9.6 
3 (4He) = 80 37.2 SI 
HF (*Li,) = ROD 356.4 30.2 
Event 3-26-1. %Liy—n-+n+2*He; By = (4.93.9) MeV, P,= (212417) MeV/c. 
ib (Ge) == BOY 312.0 341.7 
2 (4He) = 300.1 53.6 
3 (4He) 9.7 15347 72.0 
HF (®Li, ) 58.4 70.8 | 28.7 
Event 901. °Lix->x-+n+2*He; By, = (6.41 42.88) MeV, P,,= (130.2+19.2) MeV/c. 
1 (x) SU 0.0 12450 
2 (*He) 83.0 245 1237 
3 (He) 20 0 Ton 39.1 
HF (®Li,) 49.0 22.0 5.0 


RIASSUNTO 


Dall’analisi di 134 decadimenti mesonici univocamente identificati si migliora la 
precisione con cui sono note le energie di legame degli ipernuclidi HA, *H,, ‘Hey, 
5Hey, Lin, SLia, °Lia e *Bex. Inoltre vengono descritti esempi individuali delle nuove 
specie "Hey, !!B, e !2B,. I dati presenti sono combinati con quelli raccolti nella ras- 
segna EFINS. La struttura a multipletti di spin isotopico degli ipernuclei leggeri è 
discussa con riferimento alle informazioni dedotte dalle energie di legame. 
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Hyperfine Structure in Paramagnetic Free Radicals. 


G. CINI-CASTAGNOLI (*) 


Istituto di Fisica dell Universita - Roma 


(ricevuto il 17 Luglio 1959) 


Summary. — Paramagnetic resonance of solutions in dioxane of the free 
radical [(C,H,)sN]* ClO, has been investigated. Hyperfine structure has 
been found which can be interpreted in terms of the interaction of the 
unpaired electron with the Nitrogen nucleus and with some of the 
Hydrogens on the phenyl rings. 


1. — Introduction. 


Since the discovery of paramagnetic resonance absorption phenomena in 
free radicals (1), a large number of such resonances have been reported (?). 
Researches in this field have been concerned either with single peak resonances 
in solids (polycrystalline powders or single crystals) or with resonances in 
solutions which show hyperfine structure. The single peaks are often extre- 
mely narrow and are useful for precision magnetic field measurements in regions 
where proton resonance methods can not be used conveniently. The study 
of resonances exhibiting hyperfine splitting shows great promise as a tool for 
the analysis of the structure of molecules in which the unpaired electrons 
interact with a large number of atoms: for a number of molecules resonance 
experiments have provided valuable information on the distribution of electronic 


(*) This work has been supported by the International Federation of University 
Women, Mary E. Woolley Fellowship and was carried out at the Physics Depart- 
ment, King’s College in the University of Durham, Neweastle upon Tyne, England. 

(1) A. N. HOLDEN, C. Kirrez, F. R. Merrit and W. A. YAGER: Phys. Rev., 77, 
147 (1950); C. H. TowNES and J. TurKEvICH: Phys. Rev., 77, 148 (1950). 

(2) See for instance Discuss. Farad. Soc., 19 (1955), in particular the papers by: 
a) G. E. PAKE, 5. I. WEISSMANN and G. TOWNSEND: Pp. 147; b) E. E. SCHNEIDER: p. 158. 
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wavefunctions over atoms having magnetic nuclei and hence on the positions 
of these atoms and on the nature of the bonding structure in the molecule. 

In the present paper the results obtained on solutions in dioxane of the 
free radical [(C,H,),;N]}*ClO, are reported. 


2. — Paramagnetie resonance absorption in solutions of free radicals. 


The absorption of microwaves by free radicals arises from the fact that 
these compounds possess unpaired electrons When a molecule of a free radical 
is placed in a steady magnetic field B, its unpaired electron may assume either 
a parallel or an antiparallel orientation with respect to the applied magnetic 
field. If in the first instance the interaction of this unpaired electron with its 
surroundings is neglected, the energies of the electronic Zeeman levels are 
A, — + 498B, E_ — — 3gPB, E, are the eigenvalues of the « spin Hamiltonian »: 


(1) #= 98 B-S, 


where 5 is the Bohr magneton, g is the spectroscopic splitting factor, which 
has the value 2.0023 for a free spin, and S the electron spin operator. 

If now the magnetic interaction of the unpaired electron with » nuclei 
of spin 1 is taken into account, the Hamiltonian (1) becomes: 


(2) A=gPB:S+) A,1,-S, 
1 


where A, describes the (isotropic) interaction with the i-th nucleus and is given 
essentially by the product of the magnetic moment of the nucleus and the 
magnetic field at the nucleus produced by the unpaired electron. This leads 
to a splitting of the electronic Zeeman levels into hyperfine components and 
consequently to a hyperfine structure of the resonance spectrum. 

In a poly-atomic free radical in liquid solution the h.f.s. interaction of a 
given magnetic nucleus is determined exclusively by the density of the elec- 
tronic spin magnetic moment at this nucleus i. e. by the value of the wave 
function of the unpaired electron at the position of the nucleus (2) 2 Phisais 
a consequence of the molecular motions in the liquid which cause the magnetic 
field arising from the electron spin distribution outside the nucleus to be 
averaged out to zero and it implies that the interaction of the unpaired electron 


with its local nucleus, i.e. with the nucleus of the atom in the molecule on 


which the electron is localized, arises exclusively from that part of the electronic 
orbital which has s-character. 
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The resonance spectra of tri-phenyl-methyl 
PEUT | 3 

\ >| —C 
\_/ li 


which consists of a complex pattern of a very large number of very closely 
spaced lines and of di-mesityl-methyl (8) 


are assumed to arise from the h.f.s. interaction between the unpaired electron 
and the magnetic moments of some of the H-nuclei, since the magnetic moment 
of 2C is zero. The same considerations can be applied to the structures observed 
by WEISSMAN et al. (4) and by VENKATARAMAN et al. (5) and also to the reso- 
nance spectra in X-irradiated plastics investigated by SCHNEIDER (?) which 
have been ascribed to radical end groups of the type 


On the other hand, the hyperfine structures of free radicals containing 
besides protons other nuclei with non zero magnetic moments have been ge- 
nerally associated exclusively with the latter. Among such compounds we 


recall the sodium trimesitilboron (°) 


Ha 


Jou) |a 
CH, 


containing |B, the triarilaminium perchlorates (7) and various hydrazil com- 


pounds (8), containing IAN, 


(3) H. S. JARRET and G. J. SLoan: Journ. Chem. Phys., 22, 1783 (1954). 
(4) S. I. WEISSMANN: Journ. Chem. Phys., 22, 1135 (1954). 

(5) V. VENKATARAMAN and G. K. FRAENKEL: Journ. Chem. Phys., 23, 588 (1955). 
( 


6) S. I. WEISSMANN, J. TownssEND, D. E PAUL and G. E. PAKE: Journ. Chem. 
Phys., 21, 222 (1953). 
(?) O. R. GicLiam, R. 


1540 (1955). 
(8) H. S. JARRET: Journ. Chem. Phys., 21, 761 (1953). 


) 
4) § 
) 


I. Wazxer and V. W. Comex: Journ. Chem. Phys., 23, 


a 
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The present investigation has been undertaken in order to determine 
whether the unpaired electron in the aminium ion 


(OU 


apart from interacting with the nitrogen nucleus, interacts also with protons 
on the phenyl rings which would indicate how far the wave function is ex- 
tended over the molecule. 


3. — Experimental procedure and results. 


The main features of the sensitive paramagnetic resonance spectrometer 
which was available in the Newcastle Laboratories for the experiments on 
radicals, have been described by SCHNEIDER and ENGLAND (°); it consists of 
a microwave bridge and a super-heterodyne receiver for detecting the unba- 
lanced signal caused by the magnetic absorption. The waveguide bridge is 
of the hybrid ring type and is designed to operate at frequencies near 9500 MHz. 
The signal power is obtained from a reflex klystron oscillator, type CV129 
operated from HT batteries. A small amount of power is coupled into a ca- 
librated absorption type cavity wavemeter, which enables the frequency to 
be determined to +1 MHz. The magnetic field is provided by an electro- 
magnet having 5in. pole faces with a Rose rim and a gap ot 1.25 in.; at 
3000 gauss the field is uniform to within less than 2 gauss over the volume 
of the liquid specimens used in the present work. Two additional coils on 
the magnet are supplied with alternating current from the 50 Hz main so 
that the field is periodically swept over the whole of the resonance region and 
the spectrum is displayed in the usual way on an oscilloscope by connecting 
the signal derived from the superheterodyne receiver on the y-plates and a 
voltage proportional to and in phase with the field modulation to the y- and 
x-deflection plates of the oscilloscope respectively. 

The most important improvement which had been developed by Schneider 
(unpublished) since the earlier work, is a system for the frequency stabilisation 
of the signal klystron referring to the tuning frequency of the measuring cavity. 
It is based on a 465 kHz frequeney modulation of the klystron (frequency 
deviation 100 kHz=.001%). The discriminating action of the measuring ca- 


vity leads to a small 465 kHz amplitude modulation if the mean klystron 


frequency deviates from the tuning frequency of the cavity. The consequent 


465 kHz signal appearing at the output of the superhet receiver is converted 


(*) E. E. SCHNEIDER and T. S. ENGLAND: Physica, 17, 221 (1951). 
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in a phase sensitive detector into a d.c. frequency correcting voltage which 
is applied to the reflector of the klystron. The stabiliser reduces frequency 
noise, interference and tuning fluctuations arising from mechanical vibration, 
which are the factors limiting the sensitivity tor the detection of weak reso- 
nances. At the time at which the work described here was carried out, it was 
possible to detect resonance signals on the oscilloscope (band width about 3 kHz) 
from less than 0.2 ug of crystalline DPPH which is equivalent to about 10!* 
electron spins (the signal is inversely proportional to the mean width which 
is about 4 gauss for DPPH). 

For the study of the resonances of liquid solutions of radicals at room 
temperature a rectangular cavity having a tuning frequency of 9526 MHz 
was constructed. It consists of a section of standard wave guide, one guide 
wave-length long, and is excited in the TE,,, mode through a coupling hole 
in the end wall dividing it off from the main wave guide system. 

The specimens were contained in small glass capillaries of 2.8 mm internal 
diameter which were introduced into the cavity through a hole in the centre 
of its narrow side. The specim2ns were thus in a position where the electric 
field is a minimum and the microwave magnetic field a maximum. 

The solution of the aminium perchlorate was prepared by dissolving the 
weighed solid in a known volume of purified dioxane. The concentrations 
investigated ranged from about 0.01 to .32 molar. It is difficult to give the 
actual concentration of radical ions because, similar to some of the neutral 
radicals, the aromatic aminium radical ion tends to take up an extra electron 
from any reducing agent, so as to be converted into the electrically neutral 
and non-magnetic amine. But since this reducing reaction proceeds slowly 
with a reaction time of about 30 minutes while all operations including the 
actual measurements were performed within less than 10 minutes from the 
moment of dissolving the solid, it is reasonable to assum> that during the 
resonance experiments a substantial fraction of the dissolved molecules was 
present in the form of free radical ions. A comparison with the resonance 
signal, from a small crystal of DPPH placed on the outside of the capillary 
produced simultaneously with that of the amine solution showed that, at the 
lower concentrations, about 10! electron spins were present in the solution. 

Table I summarizes the experimental results for different initial concen- 


TABLE I. 
> _ = ì = a ; | 
Concentration Full width at Structure | 
pure solid powder | 14.8 Gauss unresolved | 
| 
dioxane solution 0.32 molar | DI » unresolved | 
dioxane solution 0.07 molar 25 ) partly resolved | 
dioxane solution 0.05 molar 25 » 8 peaks | 


= 14 - Il Nuovo Cimento. 
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trations. Some of the absorption spectra as displayed on the oscilloscope are 
shown in Fig. 1a, b, c. The horizontal scale is calibrated in terms of magnetic 


(a) solid CH, 


CH, N*C 00, 


CH; 
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field strength by means of proton resonance signals displayed on the second trace 
of the double beam oscilloscope (see Fig. 15). The coil with the water spe- 
cimen of the radio frequency nuclear resonance spectrometer is placed close 
to the wave guide cavity in the field of the electromagnet (Fig. 1b is a double 
exposure at the same field with proton resonances at 14553 and 14346 MHz 
respectively). The variation of the overall width of the radical resonance in 
solution does not vary as much with concentration as in other tree radicals 
(e.g. DPPH). At all concentrations the width of individual lines can be ac- 
counted for by spin-spin interaction. 

The pure solid powder gives a single resonance line width a g-factor of 
2.0033 (determined by comparison with the resonance positions of a single 
crystal of DPPH in different orientations and also by means of proton re- 
sonance). In dilute solution this line splits into 8 peaks the four central ones 
having about 4 gauss separation and those in the wings about 6 gauss (see 
Fig. Le). 

The spectrum can be analysed in many different ways. The unpaired 
electron apart from interacting with the nitrogen nucleus definitely interacts 
also with protons on the phenyl rings, but the spectrum is not resolved enough 
to give a full interpretation of the radicals’ structure. 

Two reasonable interpretations of the hyperfine structure are the following: 


a) The electron interacts mainly with the protons in ortho and meta 
positions and, of course, with the nitrogen. If one assumes that the meta to 
ortho ratio is 2:3 and that the meta to nitrogen is 1:2 one gets 39 lines with 
the intensities 1 0 6 6 16 36 41 96 121 176 272 306 436 512 581 712 736 
812 832.... The corresponding spectrum is in reasonable agreement with the 
experimental one if the width is ~ 2.6 gauss and the separation between indi- 


‘ 


vidual peaks ~ 2 gauss. In formula (2) 


A = 6 gauss , Ape = 4 CAUSS AS Causse 


ortho <—— meta 
b) One could assume that hyperconjugation is so effective that the inter- 
action with the protons on the para methyl group is predominant; than one 


gets 14 lines with intensities 
i) it, CBR 0A 21092100. 
if it is assumed that the nitrogen splitting constant is twice that of the methyl 


protons. A, = 8 gauss, methyl constant A, = 4 gauss and width of individual 


lines ~ 8 gauss. 
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ri RIASSUNTO 


Si studia la risonanza paramagnetica di soluzioni in diossano di radicali liberi 
[(CH,};N]:C107. Si trova una struttura iperfine che può essere interpretata mediante 
l'interazione dell’elettrone spaziato con il nucleo di azoto e con alcuni dei nuclei di 


idrogeno dell’anello fenilico. 


208 


CNT 


IL NUOVO CIMENTO NOT CV ENS) 16 Gennaio 1960 


Internal Motions of Relativistic Fluid Masses. 


F. HALBWACHS, P. HILLION and J.-P. VIGIER 


Institut Henri Poincaré - Paris 


(ricevuto il 4 Agosto 1959) 


Summary. — On the basis of certain general results established by Wrys- 
SENHOFF (!), MoLLER (2), Boum and one of us (J.-P: V.) (3) on the be- 
haviour of relativistic fluid masses enclosed within time-like tubes, we 
study internal motions corresponding to simple physical assumptions. 
We introduce a general new Lagrangian and Hamiltonian formalism in 
terms of Einstein-Kramer variables and relativistic Euler angles and apply 
it to some typical physical cases. 


Introduction. 


In a set of recent papers (*) the physical behaviour of relativistic fluid 
masses extended in ordinary space (that is enclosed within a time-like tube 
in Minkowski space) was studied in order to determine the laws of motion 
of certain global quantities. In the present paper we shall attempt a finer 
analysis of what goes on inside such tubes in order to reach a physical under- 
standing of possible types of internal motion. 

Such an attempt is necessarily limited. It is clear that the general motion 
of the fluid inside the tube can be as complex and violent as one pleases: since 


(1) WeyssENHOFF: Acta Phys. Pol., 9, 8 (1947). 
(2) C. Morten: Ann. Inst. H. Poincaré, 11, 251 (1949). 
(3) A. Boum and J.-P. ViciER: Phys. Rev., 109, 1882 (1958). 
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it varies with the internal tensions and initial conditions considered. More- 
over, any departure from the point-like model implies that we dispose of an 
unlimited number of degrees of freedom and variables corresponding to the 
unlimited number of possible motions inside our time-like tube. As it is evi- 
dently hopeless to describe in detail such a motion, we must necessarily make 
some simplifying assumptions to reduce the number of variables needed. In 
order to do this let us remark that our present objective is to determine the 
behaviour of global quantities averaged over the total volume of the fluid 
mass. These correspond in a sense to collective variables describing overall 
properties of the motion. If we accept this point of view our first objective 
must be to introduce such typical variables associated in a clear way with 
some average properties of the fluid mass. After that, with the help of some 
physical restrictions, imposed on the internal motion, we can deduce in prin- 
ciple from a Lagrangian formalism the general laws of motion. Only then it 
will be possible to study simplified cases of stable internal motions which cor- 
respond to relativistic generalizations of Poincaré’s classical theory of non 
relativistic fluid masses (1). 

This brief analysis determines the plan of our paper. In Section 1, we 
shall summarize known results of the general theory of relativistic fluid masses, 
and introduce new ¢ global » variables. In Section 2 we shall study the Lagran- 
gian and Hamiltonian formalism. In Section 3 we impose physically reasonable 
restrictions on the motion and determine the corresponding equations for the 
internal motion with the aid of a suitable Lagrangian and Hamiltonian. In 
Section 4, we shall recall the usual procedure followed in the case of a non 
relativistic rotation, in order to put into evidence the constants of the motion 
by the introduction of Euler angles. In Section 5, finally, we shall generalize 
the above procedure to the internal motion of a relativistic rotator in terms 
of relativistic Euler angles and determine the constants of the IRE 
motion, which paves the way to its future quantization. 


1.- As was stated in a former paper on the subject (*) no specific assump- 
tion on the exact nature of the internal fields constituting the fluid masses 
studied here is needed at this stage. They could be built for example by 
electromagnetic field cone entration, combination of electromagnetic and gravi- 

tational fields (such as the «geons » proposed by WHEELER (°)) or by still 
more complex combination of known and unknown fields. The only thing 
we shall suppose is that they are built out of an energy-momentum density t,,, 


(n =%,) and current density j,, strong within the particle and negligible ae 


(*) H. Poincaré: Acta Math., 7, 259 (1885). 
(9) J. A. WHEELER: Phys. Rev., 97, 511 (1955). 
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side (5), both satisfying the conservation relations: 


(1a) o.t_=0 


(15) o 


inside the tube which limits the particle. 

As stated before, the theory of global variables has already been developped 
by various authors. If we assume that our particle is isolated with no external 
force, MoLLER has shown that it is possible to determine a physical point 
inside the fluid mass which corresponds to the classical picture (?) (when the 
tube is reduced to a line for example). 

Let us resume his results. MoLLER first defines the total energy momentum 
of the particle by the relation: 


Ge = far k 


y 


where dV denotes a volume element of a plane three dimensional section 2° 


of the tube and the | is taken over the whole section. He then shows that 
because of equations (1a) this G, is a four-vector independent of the chosen 
section. In the absence of an exterior force its components are constants of 
the motion: 


ASS 
bo 
= 
SF 
I 


Then with the help of this vector he defines a particularly remarkable frame 
II, (which we shall hence forward coll Moller frame) namely the one in which 


Gy =. 
As we shall see, this frame will play an essential role in all that follows. 


We have now the possibility of introducing a first physically remarkable 
point: the center of gravity (X.,,) defined in //, by the relations: 


xi far = [et dv 5 NE = 0 . 


| : Se) Sas J Cato, mes at PN 
This point moves with a four velocity w,, GO c?) parallel to Gi In 


(6) F. HALBWACHS: Théorie des fluides relativistes (Paris, 1960), Chap. II and 
Appendix A. 
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troducing the scalar quantity M by 


GG a= al 625, 
we get 
ee I 
G,= Mu,. 


Next we can define an angular momentum M, with respect to that center 
by the relations: 


M°, = | HR ere 


This can be also shown to be à skew-tensor. 

MoLLER then demonstrates that if we call 7 the proper time taken along 
any arbitrary world line followed by the center of gravity we have (calling A 
the expression dA/dt): 


M =0, 
= 
MY = 0), 
ù, 21 


The determination of //° can now be made with the help of four unitary 
and orthogonal vectors, which we shall denote af: (the index (£) labelling the 
vectors varies from one to four and is not à tensor index). We write: 


ry (4) % (2) Gas . * (4) 
Og eee WLU ae 
ic Me 


iy Mapa il Bree (3) Le È c 
The third space like axis a, (which can be called the Z axis) can be chosen 
parallel to the «spin vector» corresponding to the angular momentum MP, 
6 5 a n uv 
which is defined by the relation: 


ali = es . 4 x 
SÙ — — 0 SRO 0 (4) Tr 0 Pe v 
Su 2 Mu, = Mia (is, is the dual 5 Enraÿ Ms) x 


This choice is justified by the fact that Min =0 and M°,=0 so that 


3) MO — 1 ,0 (4) 
( Me, — PASTORI, sy Ag 
and 


ow 


cu 
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The last two remaining space-like axes must be taken orthogonal to a 
and an. For example if we denote by À, the components of the four vector 
joining the center of gravity with another physically important point (such 
as the center of matter density which we shall define further) we can deter- 
mine a vector 


ft 


1 
= P Cir PE Oa?) à 


which is evidently orthogonal to af and a’? if we write P— Ra and 
Q=R, à. We can complete our orthogonal and unitary frame //° by the 
vectors 


(1) 


a nd 2 (È 4 
ae and ‘al = te sap dy 0 Vago. 


x AB 


As one knows, the orthogonality conditions on the ai system can be ex- 
pressed by the equivalent relations 


A l 
afaf = Ô 
Lu 


v ur) 
(4) : 
[ee Oe ie 
uo 
The essential point of the preceding discussion is that one can define a 
point X and a frame //, centered on X independently of the detail behaviour 
of the internal structure of the particle. This implies: 


a) that one can reduce the behaviour of the particle to the consideration 
of these quantities when we neglect this structure and approximate it as à 
rotating point; 


b) that they provide a convenient frame of reference to study these 
internal motions. They play in a sense the well known role of the classical 
inertial frame in ordinary mechanics. 


Now, if we want to study the internal motions let us first remark that they 
depend on so many factors (internal stresses, tensions, etc.) that a detailed 
study is a priori hopeless except in simplified cases. This is also true in clas- 
sical mechanics, where studies of detailed fluid motions inside fluid masses 
have only been carried out in particular cases with special geometrical sym- 
metries (4). 

In general, the kinematical behaviour of the matter inside our tube is 
represented by à bundle of world lines tangent at each point to j, and com- 
parable to a world-rope. Its complete description could be attempted in 


two steps; 


TTT es 
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a) determination of an average world-line { considered as the path of 
the ropes’ central geometrical point; 

b) average motion of the bundle around that line. 


Evidently, if the dimensions of the tube are small and if we assume that 
the lines of this bundle are twisted in a continuous way around 1, we can avoid 
step b, and restrict ourselves to the determination of { considered as describing 
the average twist of the global motion around the path of the center of gre svity. 
Now as was shown before (3) the determination of l is very easy if we assume 
it to be the world-line followed by the center of matter density. 

Such a center, as we know, can be defined in //° by the relation. 


Y,{44V= [eu ay. 


Di 


This point moves with the velocity 
ESSI (Y,Yu=—@2) 


and we know (*) that if we call 7 the corresponding proper time along the 


world line followed by Y,, R, the vector Y, — X,, we can write the angular 


momentum with respect to it in the form: 
My = | ru bin — (0, Yy)tu dV = 
== fee, “n Xy)tys Pi (2, ae Xy)t ua] ay +fta- Ya (X,— Y)t ua) | dv Fi 
me 1%, G,— Rk, G,. 
Then, if we take the derivative of both sides with respect to t we get: 
(5) U,=G,I, GY 


S( 7 av È È È se SÌ i 
> we have for Y the set of classical equations (5) 


[ à =), T2 72 
(6) da in | | 
Mega, (kia 


Now the motion of this center of matter density and the corresponding 


0 Le 
global quantities M,, and Ja can be determined by various sets of kinema- 
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or 


tical parameters. Following a well-known method we shall first attempt to 
describe them with the help of the Einstein-Kramer variables, that is we shall 
parametrize / with the help of a second orthogonal tetrad bs reflecting global 
properties of the whole droplet, referred this time to the center of matter 
density. 


We first write: 


b® = kl 
H ic ? 
S 
(3 adi di è 
(ia = a with s,= My Y, 
Wegs 


and determine a vector B, by the relation: 


1 
B,=b0-—= 


Est Fe) 
ai (Ry — 008) 


with P'=R,b, and Q'=R, bY. Then denoting by b the normalized vector 
B,|/vB,B* we have defined four orthogonal and unitary vectors b?, 09, DI 
and bt? with 


a 
CG) (2) 16) 7, © 
Di 3 Envap Py by be 5 


We have also evidently: 


(bibi # Ol 
(7) ue wie 
less 


and will denote by 2° the corresponding frame. As is well known (7) the in- 
stantaneous rotation of this frame is given by the expression &,, such that: 


b, = Ow Bs 
with 
(8) Ony (65 b§ — b5b°) = bb 


Physically, the ©,, can also be interpreted as the projection of b, on the 
db, direction since bei — = D, 

The w,, tensor aiden describes the instantaneous rotation in 2 close 
neighbourhood OL I, 

Following TAKABAYASI (*) we then define the instantaneous rotation vector 


(7) F. Gürsey: Nuovo Cimento, 5, 784 (1957). 
(8) T. TAKABAYASI: Nuovo Cimento, 18, 532 (1959). 


Le) 


è 


Frs Ur. 


216 T. HALBWACHS, P. HILLION and J.-P. VIGIER 


and acceleration by the set of relations: , 


a) 
| Ds RE D 0, 
(4) (4) 
RO ae Dip Oy 
where 
D = 1 € Come 
uv 2 pap xp 


If we now move along / and denote by 7 the proper time along this line, 
the frame x, will undergo a Lorentz transformation with coefficients depending 
on t. As we said before, it is very natural to study the relative motion of 
(Y,) and Z, with respect to (X,) and II, and write 
n | Ve Lyra: 

LE, = Ne 
In a sense this discussion shows that within the frame of the preceding 
physical assumptions (when we neglect step b), the internal motions of a rela- 
tivistic fluid droplet can be described by the relative behaviour of two frames 
centered on two points namely. 


—2° (b;) with its origin on (Y,), 


== ill (47) » » » es 


As we have shown, in the case where no external force is acting, //, and 
(X,,) move in a straight line, while 2, and (Y,) spiral around it according to 
the laws of motion. Our model can thus be compared to a relativistic hydrogen 
atom where the proton and a fixed frame centered on it play the part of (Y,) 
È y’ 7 7 : . . a = 
and ee ANOLE: .) representing the motion of the planetary electron (with 
its classical « spin »). It can be considered as a realistic representation within 
usual space-time, of the bilocal theories introduced by YUKAWA (°), RAYSKI (1°) 
6 i dra re (11 TANAY a 1 1 
and Oper vuthors (!!). However, such theories which rest on the double point 
model (without the two frames) can be considered as a degenerate case of the 
proposed model of relativistic rotator. 

Our next step is evi 7 press i i i i 
vi ee step nou to express a Lagrangian formalism in //° in 

cs 4 pnd the bis. But before we do that we would like to remark that 
our hydrogen-like model already indicates qualitatively how these motions 

(°) H. Yukawa: Phys. Rev., 91, 415 (1953). 

(1°) J. RAYSKI: Proc. Phis. Soc., A 64, 957 (1951) 


(11) O. Hara, T. MARUMORI, Y. OHNUKI 
Mi. NUKI, and H. SHIMODATRA: Prog. Theor. Phys. 
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could be quantized, if we recall de Broglie’s ("*) early interpretation of Bohr’s 
model for such an atom. Indeed, we shall show in a later paper in collaboration 
with Professor BoHM that stable or metastable types of internal motions cor- 
respond to periodic motions of (I) on quantized orbits (corresponding in a 
sense to Bohr’s original orbits) while 2, comes back periodically to the same 
relative orientation (with respect to //,) after an entire number of circuits. 


2. — According to our program, let us pass now to the Lagrange formu- 
lation of the laws of motion. 

This can be done easily if we follow GÜRSEY (7) and choose as variables 
the co-ordinates of the center of matter density Y, and the components bi 
which determine X,, all these variables depending on the proper time 7 along 
the world line 1 followev by Y,. But before we give any explicit form for / let 
us recall briefly two general results. 


A) If we suppose that the internal structure of a relativistic body is 
expressed by a set of vectors q, expressing certain global properties (instan- 
taneous rotation, etc.) of this structure and localized at a definite « center » 
of the body, we can assume the free motion will be described by a Lagrangian 
function (1°): 

L(Y 5 do du) » 


with 
Mi 
(9) o ft ie qin) dt=0, 
M, 


where À denotes as usual the expression dA/dr, t is the proper time asso- 
ciated to the world line followed by the « center ». The summation is taken 
along this world line from a fixed point M, up to a fixed point M,. 

As a consequence of (9) we get Buler’s equations 


ah _ a 
gos dr Og, 04 
d= 01 ae 


10b aa 
00) drcY;, 


Note that (10b) can also be written in the form 
: OL 
C= 0 with K,=—.- 
(11) Sor 


(12) L. DE BROGLIE: Thèse, Paris (1927). 
(13) F. HALBWACHS and J.-P. VicieRr: Compt. Rend. Acad, Sci., 248, 490, 1124 (1959). 


n 


a 
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If we now apply an idea of E. NoETHER (!4) and consider variations SY,» 
dg, resulting from an arbitrary Lorentz transform 4,,, they will no longer be 
iaia and on another hand, they will no cer vanish at the limits 


My 
of integrations. However, the variation à fz dr will still be zero since the 
Lagrangian is a scalar. We thus get: =e 
My My 
i IE OL 
> 
5] bdr =| = SY, | Ba" sates FF LR = =). 
Ova Wu du 
Mi Mi 


which gives, when we integrate by parts: 


d CL OL d OL 
ne | 1 4 # Nor 
i È dr ay, Je | (a dr ci) a 


‘OL SY 
sl vert] Age = di) dr = Om 


The first two terms vanish as a consequence of Euler’s equations so we have: 


È (= svt EE y sg EST 
dc\ay, “où =i; de oie oe 


If we write down explicitely the variations 3Y, dg, resulting from an 
arbitrary infinitesimal Lorentz transform, and set: 


Sa OL OL 
= ia 04; Ul dgr 


we get finally, taking (11) into account 
(12) Se aii, Y,- HAT 


Then if we identify the « canonical conjugated momenta» # and S_ with 

| OT: = chi 
the integrated global quantities G, , and M, , of the droplet theory in coon 1 
» 


Eee (11) and (12) are just ans two wine. equations of motion (2) 
and (5). | 


B) If we have a Lagrangian of the considered type which depends on 
various vector variables g° and is a quadratic form of the 4, we can write the 
equation of motion in a relativistic Hamiltonian form. Introducing the conju- 


(14) E. NoETHER: Gotting. Nachr., 235 (1918) 
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gate variables by the relations: 


OL 
13 es 
( ) Px ed, (das Gx) 5 


we can introduce the scalar Hamiltonian H = pq — L(g,p,) and find, taking 
into account Euler’s equations: 


d ee hr = OL 
dr x) CT 
the usual relations 
: oH oH 
14 = — —— £ o, = —" 
(14) Pa 3a and 4, apt 


Note that H(p/,07) is a constant of the motion, since 


ge Ae AN. 
ê ee dute oa 


because of (14) when we replace in L(g7, 4,) the gs by their values in terms 
of the p’, deduced from (13). As we shall see in the usual relativistie theory, 
H is just proportional to the rest mass term and we have H =— 3m. 
Moreover, introducing the usual Poisson brackets 


[f(4, p), 94% p)] = 


we see that the equations of motion (14) imply that the variation of any 
function A(p, g) along { will be described by the usual relation 


. LA 
Meso A HI, 
dt 


each quantity 4 whose Poisson bracket vanishes being a constant of the motion. 
The preceding formalism can evidently be applied to a Newtonian point 
particle. If we denote by X(t) the mechanical variables characterizing its 


position we can choose the Lagrangian 


SS va 
L=imX,À 


4 
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which gives as Euler equation 
d : 
dz (MX n) = (i), 


The corresponding Hamiltonian equations follow easily. We have: 


which imply 


So 
lal = CL = 3G, DE 


or in terms of canonical variables: 


1 
Sci 


2Mo 


Hea 


The canonical equations of the motion are: 


G=0 and G, == mot, 


iG 


and according to the latter, we get 
Ha iO 


If we now come back to our problem of determining a Lagrange formalism 
for the extended particle the most natural assumption is to add to the point 
particle Lagrangian some terms which will globally take into account the 
internal behaviour of the particle. However, we must keep in mind the es- 
sential fact that any specific form of L represents only certain possible simple 


particular behaviour of the extended body among the unlimited number of 


possibilities (related to the detailed internal behaviour) so that no form chosen 
‘an correspond to a general theory. Any specific choice must therefore be 
guided by physical reasons and justified by its experimental consequences. 


3. — Our specific proposal is to add to the point particle Lagrangian the 
relativistic generalization of the classical term representing the rotation energy 
described by the motion of the fram 27 with respect to //,. 


In three dimensions where Z, is represented by three vectors 87 (bb = 6,,) 
5 Da x . 4 | è ne) ki 
and the angular velocity is the antisymetrical expression @,,— b;b;, we know 


i=) 
a 
i 


oa" 


INTERNAL MOTIONS OF RELATIVISTIC FLUID MASSES 


bo 
bo 
fai 


that such a term can be written in the form 
£1” dio; 

I being the symetrical tensor of inertia with 

O, = Erin Dire 


As we have shown elsewhere (1), the relativistic generalization takes the 
form: 


the relativistic inertial tensor I,,,, being symetric with respect to the change 
of u, v to x, f and antisymetric in both pairs x} and pr. 

Naturally, this term can be simplified in many physical cases and we can 
assume various simple forms corresponding to typical cases. For example 
NAKANO has proposed (!°) to use the form Loto which represents a natural 
relativistic generalization of the non-relativistic spherical rigid body, I being 
taken constant. Introducing the EX variables and taking into account re- 
lation (7), T takes the very simple form: 


We shall treat this simple case completely with the help of two sets of variables 
in order to illustrate the preceding Hamiltonian formalism. 

The complete Lagrangian will contain the ordinary mass term imoY, Jo 
and à set of terms entailing the orthonormality condition (7), which can bo 
written : 

(0505 0, 
with the symetrical set of Lagrange multipliers 2,,. But we must not forget 
that we have chosen the fourth HK vector oF Are to Ÿ 5 that is tangent 
to the world line followed by the center of Du density. This can be ex- 


pressed in two equivalent ways: 


1) We write simply, with the index (r) running from 1 to 3 only: 
N get L Rs 1; e MRO 
(15) i Fra Tu) + dm (bn vi— a al do) : 


(15) F. HALBWACHS, B. Hirrion and J.-P. Vigier: Nuovo Cimento, Letter to the 


Editor, 11, 882, (1959). 
(16) T. NAKANO: Prog. Theor. Phys., 15, 333 (1956). 


15 - Il Nuovo Cimento. 
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bo 
bo 
bo 


This Lagrangian can be treated by a well known generalization (7) of the 
Lagrange formalism with second derivatives and gives suitable Euler equa- 
tions. But certain difficulties arise if one will pass to the Hamiltonian for- 
malism. To do this we can use a well known formal device. 


2) We keep independently both variables ich} and Y, and introduce 
a subsidiary Lagrange condition which makes them to be equal, with a new 


Lagrange multiplier 4: 
(16) L=imy,Y,— HIDE +A,, (020; —6,,) +k, lied? — Y,). 


This form can be treated in Hamiltonian formalism without any difficulty. 
The canonical angular momentum derived from expression (16) is just, as 


we recalled in Section 2,. 


e OL OL il UE DE 
(17) M,,= bi — de = 51 (bib; — b5 dj.) = Loy 
di QUES 


so that .M,, is parallel to ©,,, which corresponds to the well known classical 
behaviour of the spherical rotator. 


Then introducing the canonical variables: 


aL 


| ar NU 
A OL 12 
i (yy abi “i 15 De ’ 


we get the Hamiltonian 


that is, by elimination of the time derivatives: 


vY al E DE à Pa 
H = 3 (Gut hu)? — BEBE — th be — Au (DÈDE — da)» 


CH oH 
= EST O IL 
A TUE a = mg Guth =, 
oH , hs. 2 
FASE Se) if ESSO n o 
RS, pil, 
le 


(7) M. Von Lave: La Relativité (Paris, 1924), 
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cH . CH 9 
D DS ari ve O) MEA 7, (4) 
Soraya dl av VUS Pus r= Du oi 
Ob, 4 ope. ui H wd 
OH Pr D pore aie 
ei a È ay LC) 7 
Ok, = Mo (Gu | 11) ich, = 0 ’ — that is ich, = VO 3 
CHI? 
21 = bibi On 0 

uv 


Naturally, this Lagrangian could be complicated in many ways be adding 
various Lagrange conditions and taking various forms for the inertial tensor. 
We could also add to the Lagrangian some terms depending on our two centers 
of the form V(X,), this potential between the two points corresponding to 
the global effect of the internal stresses of the droplet. Of course, it is not 
necessary to use the bi variables; any equivalent set would do. For example, 
if we study the motion in the frame //, we can introduce the relativistic gene- 
ralization of the Euler angles which transform //, into 2. This will be done 
in the next Section. In fact, we shall see that this new set of variables allows 
a elear physical understanding of the relativistic behaviour of the two frames 
and facilitates the comparison with the non-relativistic behaviour of the clas- 
sical body in three dimensions. 


4. — In order to clarify the exact behaviour of our spherical four dimensional 
rotator, let us first summarize briefly certain results concerning the three di- 
mensional rigid spherical rotator. As one knows such a rotator is described 
by three orthogonal unitary vectors bd; attached to the body and satisfying 


the conditions: 
(18) bib, = One, Dy ars 


The instantaneous rotation of this frame with respect to a fixed laboratory 
frame a’, is described by the pseudo-vector ©, defined by the well known 


relation: 


(19) . i= US (the dot means the derivation d/dt). 
The rotation energy in the spherical case takes the well known form: 


(20) T= 5100s == 5850; » 


where I represents the classical moment of inertia and s, the angular mo- 
mentum which is just Zw,; taking into account relation (19) this becomes: 


as erie ig fee 
P= SLE 551m did di Om = 3Ib; bi; 


ce m 


oO 


al 
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which is just the non-relativistic limit of the internal energy taken into account 


in the Nakamo Lagrangian. 
The total Lagrangian can be written: 


L = 4166" + 2,,(050; — 05)» 


the 2,, being a symetrical set of Lagrange’s multipliers implying condition (18). 
We pass to the Hamiltonian formalism with the help of the conjugated 
momenta: 
LA OL = i . 
Pi ak 


and form the Hamiltonian 
Mor eee 
Ho 
which becomes in our case 


He È BF Bs — A(b5b5 — dai) - 


The canonical equations can then be written: 


CH 1 . 
Me — b” 
de oo 
CH b :, 
ELA => Aii05 DI . 


The total angular momentum s,s;=1I2©0;@;= ff; is evidently a constant of 
the motion since [H, s,s,;| = 0. 

As one knows this formulation can be 0. simplified by the use of the 
Euler angles, which transform aj, into dI. 

Writing this transformation in the classical matrix form: 


by cosy sing 0/|cos6 0 —sin6|| cosy sing 0 ay, 
by |=|—sing cosp 0 OT o | —siny cosy 0 ay 
be 0 O1) sin 0,0 cos 0 | 0 eel Oe 


we obtain after calculation the well known table of the direction cosines: 


vel bt = A" aq? 
ie hose Ue 

cos p cos 0 COS y — sin g sin y cos p cos 0 Sin y+sin pg cosy — cos p sind 

— sin p così cosy— cos p sing — sin y così sin y+ cos p cos y sin g sin 0 

sin 6 cos y sin 0 sin y cos 6 

=H 

a 

ist 
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“ 


With the help of these new variables and the definition (19) we get the usual 
expression for the projections of ©, on the moving frame Dy: 


wt = — cos gp sin) p + sing 6, 
(21) o'? = sing sind p+ cosp Ò, 
w'3= «+ così y. 


The detail of this calculation in a simplified matrix form can be found in 
a paper from one of us (F.H.)(!). As a consequence, the Lagrangian 
expressed in Euler angles takes the well known form: 


L=41(@+@4+ y +29 080), 


condition (18) being satisfied automatically. This gives, with the help of the 
canonical momentum 


Po=3z==I(P+%c080), 


the usual Hamiltonian in canonical form: 


2 2 € 4 
Dot Pa 2pyPy 6088 | a 
oT sin? 0 dé 


vie 


On the other hand, we can compare the usual expression 


H=4P,G¢+ Py} + Pob) » 
with the form resulting from the formula (20) 


alt wait 
H — Ls Wi 


Py— S33 

# . ! 
Po = 8, SING + 8, COSY, 

iP ee = SU E 
Py =— #1 Sin 9 cos g + s, sin0 sing + 8, cosO, 


(18) F. HALBWACHS: to be published in Ann. Inst. H. Poincaré. 
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or conversely 


7 tg 0 a 
A af Oe, CULO) COR) = [Uy ee — 
S1 = Pe Sin y + Pp, ctg ? ~ Pein 9? 
| i _ sing 
Sy = Po cosg— Py Cbg 0 sin gp 1 Pen” 


I 
S3 — Po) 


+ N a r 
as the projection of the angular momentum s on the axes of the moving frame di. 
The components on the fixed frame a; are (15): 


| oe sa cos Y 

4 = — Po Sin yp — p, ctg 0 cos sin 0 

Sy Po SIM Y — Py CLE SIR, sin 0° 
4 Sin y 

S = po cosy — Py cg 6 sin p+ Py =p» 


SEM 


. A . . ! RI È 
These expressions of s, and likewise those of s, allow us to express the 
total angular momenta 


al 
sin? 0 


9 a Mi NU 
Se S781 S78 


(Do + Do — 2Po Py COS 0) + pè = 21H. 


These results, compared to the above canonical form of the Hamiltonian 
shows immediately that we have: 


[H, p,|=0, [A= 0. ist 2200 


That is, we have, as it is well known, three evident constants of the motion, 
namely the projections s, of the angular momentum vector on the fixed a 
axe, its projection s, on the moving b axe and its magnitude s. 

Finally, a short calculation yields the following form for the rotation 
energy in Euler angles: 


T = I(y + Qt + 8 + 2p così). 


5. Our next and last step is to pass to the four dimensional theory. In order 
to do that, let us summarize the relativistic generalization of Euler angles. 
This generalization which we have developed in a recent paper (1?) expands 
a result obtained theoretically by ErnsrEIN and MAYER (2°) and further and 


19 + TTS = Sio = 6 = A . 
CES Hatewacns, P. Hitiion and J.-P; Vier: to be published in Ann. Inst. 
H. Poincaré. 


(20) A. EINSTEIN and R. Mayer: Sitzber. Akad. Berlin (1932). 
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One can show that the 


more precisely by CASIMIR (21) and VAN WINTER (22) 
Lorentz transform which transforms the a into the be can be written into 


the matrix form: 


ei) DA, Ps ae ee Mii dia 
Di cos 9 sin 5 0 0 COSI 0 0 
pes 02) rk 7 ) Qo 
Be sin - cos? 0 0 sin-— cos 0 0 
2 2 2 2 
(22) = : 
: mt 3 + : ) 
OW 0 0 COS a sin iat 0 0 COS ee — sn DE 
ad Li 
È gi gr oe ie 
pa 0 0 sin cos? 0 0 sin? cost 
2 2 2 2 
0% DA (i SOs 
LOS — 0 — sin — 0 COSCE 0 — sin — 0 
2 2 2 2 
0 ò 0 ne 0 2 0 ine 
COS DI Sin a. COS ©. pl id 
OF ; zo à 07 û 
varese gia Su Hie 
sin 9 0 COs 9 ( sin 9 co 9 
i sE + une = 
0 — sin — 0 COS — 0 sin — 0 cos — 
2 2 2 2 
qi je ae e = | 
cos = sin 0 0 cos” sin 0 0 a 
2 2 2 2 i 
75 wr : = = ‘ 
Sin cos = 0 0 sin — cos 0 0 di) 
2 2 2 2 
i 4 = = 
; COS — sin #_ 0 0 COS 1a = sin © Ge 
0 2 2 2 2 
is yr TATA si (a) 
\ 0 0 sin. COS 0 0 bn COs ad, 
a natural generali- 


It can be shown (!°) that this formulation constitutes 
zation of the three dimensional theory with the help of complex conjugated 


angles: 
pi = Pi + ipy 3 6& = 0, + 10; ; pe = wr + Wr . 


21) CASIMIR: Thesis, Leiden (1931). 
22) C. Van Winter: Thèse, Groningen (1957). 
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The real quantities g,, 0, y, correspond to the usual space Euler angles, 
while y, 0, y represent hyperbolic angles (varying from — co to +00) 
expressing pure Lorentz transforms. 

Following a procedure initiated by EINSTEIN and MAYER (?°) we have 
shown in a preceding paper (1°) that, if one introduces two sets of complex 


conjugated self dual tensors: 


Tre = et bi + (6760 — bb) 


pv ey 


and the corresponding tensors associated to the a, 


(4) 


Det rst 8 (Ta at 1 (4) 
aa (aia, a, a, a 


transformation (22) transforms the Af} into the 77 in terms of gt, 67 and 
y* only. 

In the same way the 47, transform into the 7%, in terms of g-, 67 and 
y7 only. 

Now with the help of these / we can build (!°) two sets of three dimen- 
sional complex vectors B’* belonging to a complex three dimensional Eucli- 
dian space. 


We have 


DM REATI (EE) rp(4) TI.(4) 
bi “n bibi va bibi SE bibi . 


eign 


If we form the corresponding fixed vectors 
r (4 L 4 
AT = aja — dat Le, Ga, 


one finds, after some calculation, that transformation (22) induces a simple 
complex Euler rotation from the A4,° to the B7* with Euler angles gt, OF, wr 
and from the A~ to the B7 wi 7 I = Oy dr | 
£ : Si to the B7 with Euler angles g-, 07, y. 

We have for example: 


+ | Une : ; 

Be | cosy’ sing* 0|| cos 4+ 0 — sin 6+|| cosyt sin yt 0| /AM* 
(2)+ Bir 4 nt ‘Mt | 7 

BE* | = |—sing* cost 0 | OM 0 — sin y* cos pt 0 a 
(3)+ | . 

Bi 0 O 1}|—sin0+ 0 così 0 CARLINA 


k 


and the corresponding expression for BI and ATE, 
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Of course, the BI* and A’* both satisfy the usual three dimensional ortho- 
gonality conditions: 


Dis pia 5 assy Ose) BN 
BA Bi dx ? Bab 0 
Att Art = 6 Ars Are — ÿre 
Tee I: ee È 


This procedure corresponds to the well known isomorphism between the 
Lorentz transform and complex Euclidian rotations established by CARTAN (°°). 

It is therefore tempting to utilize, instead of the bi, the new sets of complex 
variables BI* (and express their variations in terms of the complex Euler 
angles) in order to bring into light any analogy with the well known three 
dimensional theory. 

The best way to write the Hamiltonian formulation of the Nakano motion 
in terms of the generalized Euler angles is now evidently to start from the 
form (15) of the Lagrangian. Then in this case no matter if we put by? or 
(Llic)¥, in the matrix formula (22), the expressions in the Lagrangian will 
be finally expressed in terms of the complex angles so that the Lagrangian 
will depend only on the first derivatives of the new angular variables. 

The correspondence between the old and new variables will be expressed 


in the form: 


(1) 
DÈ, 
(2) 
bi, 
(3) 
Di, 
de. 
ie 
cos sin © 0 0 cost sint 0 0 
ce + = Pea 
ain cost. 0 0 une COS — 0 0 
ad = dl dd 
= : i I 
STE 0 30 aa 
J Be S i J Ses — sm — 
0 0 COR as 0 cos 9 9 
sor ae a RE g_ 
0 Of) sii cos a] 0 0 smi COS — 
dl -_ dl dl 


mM 


(23) E. CARTAN: Legons sur la Théorie des Spineurs (Paris, 1938), T. 1. 
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O+ Us Di » 0 
‘ i) = 0 cos — 0 — 10 
cos — 0 sin 9 Û 9 9 
G+ honk hee dim VICE 
0 COS — 0 pin 0 COS > 0 Que 
: ; i ig | kx à 
sin a 0 COs — 0 sin Ex 0 CO: D 
Sa 0+ . = = 
0 == Rita =" 0 COS — 0 sin "i 0 COS D 
di - 1 APR coe @ 
cos sin 0 0 cost sn? 0 0 dh 
SR 2 2 D 
i di E A a (2) 
sin Ld COS = 0 0 sin 2 COS ic 0 0 Uh 
t + y 7 
2 np ven € (3) 
Nel ST oe OS— — sin — a 
0 0 COS 9 sin 9 0 0 COS 9 9 i 
: : à Rees ae Do 4) 
0 0 sin È cos È 0 “sind reat aîî 
\ 2 2 0 2 2 


If we introduce these new variables in the Lagrangian (15) we see imme- 
diately that the orthonormality conditions taken in the form 


ia 
(24) RASO CA) 


wy 


are identically satisfied so that the third term in the Lagrangian identically 
vanishes. Further, because relations (24) are equivalent to 


a NSF 
b A = 0 


the rest mass term in the Lagrangian put in terms of the Euler angles, reduces 
identically to the constant — Lime. 


We are thus left with the bare rotation term: 


With the aid of the matrix form (23) w 
of the angular velocity (Oi bb” I 


n 


energy ©. 


e can express (!5) the components 
VICE and finally the rotation 
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In terms of the B°* variable we get: 


TE 


bo] 


(FRA Ba Be) 


and in terms of the Euler angles: 


GP 2 [(prye+ (OH) (0+)?+ 2ptpt cos 0++ (wr)? (D)? (07)? 29 py cos 07] i 


This is the eventual form we should take for our Hamiltonian. 

This is a very interesting result, since it shows that the Hamiltonian of the 
four dimensional spherical rotation splits into those of two three dimensional 
complex conjugated Euclidian rigid rotators. We can therefore apply to each 
of these rotators the preceding three dimensional theory (since the presence 
of the symbol i does not modify qualitatively the Hamiltonian formalism). 

Although the two complex rotators are not independent, their Hamiltonians 
are expressed respectively in terms of the two sets g*, 0°, y* and g7, Geli, 
together with their conjugated momenta p,., Por, Py+ and Pr Por Py which 
are considered as usually as independent canonical variables. Thus the Poisson 
brackets with the Hamiltonian of each quantity made up with the variables 
of one of these sets are independent from those of the other, so that the con- 
stants of the motion for the whole Hamiltonian are simply those of each 
separate three dimensional top. We may thus simply transfer to the relati- 
vistic rotator the results we have recalled in Section 4. 

Thus we can say that the constant of motion M,, + Y,G, — Y,G, put 
into evidence in the four dimensional formulation [relations (5) and (17) | 
splits into six separate ones which are the same as for the non-relativistic 
top, that is 


(ei ia Se, Say TG 


This essential result relates our physical classical theory to the mathe- 
matical work of VAN WINTER (22) and prepares the wanted quantization. 


5. — Conclusion. 


As we said before, the preceding theory does not pretend to give any 
detailed description of what goes on inside our relativistic fluid masses; its 
main result is to show. 

a) that such masses can be roughly compared to relativistic rotators 
built of two points and two frames spiraling around each other; 
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b) that the relative motion of such rotators can be described in terms 
of those of two three dimensional complex rotators. 


This paves the way to establish a quantum theory for these internal mo- 
tions. It appears already that the usual quantum procedure will introduce 
a set of quantum numbers corresponding to the preceding constants of the 
motion. Indeed we shall see in a further paper that they lead to a classifi- 
cation ot internal states which reproduces the well known Nishijima-Gell Mann 


scheme for elementary particles. 


In conclusion, we wish to express our deepest thanks to Professors L. 
DE BrogLIE, J. L. SYNGE and Dayvrp Boum for many helpful discussions and 
suggestions. 

In fact, this work is part of a common program of research established 
on the basis of common ideas with Professor BoHM, program that will be deve- 
loped in a series of further papers. 


Note added in proof: 


The choice proposed on pages 213 and 215 for the a‘? a and bi? bY? is here given 
for the sake of a simple example. In a series of further papers (in OR Acad Sc.) it 
was shown that we can use the corresponding gauge arbitraryness in view of satisfy- 
ing some important and meaningful conditions. The directions chosen are fixed in 
space for a,’ and aj, and rotate with a definite constant velocity for bf and b?. 
Thus the choice is very different from the one here proposed. 


RIASSUNTO (») 


Sulla base di alcuni risultati generali ottenuti da WEISSENHOFF (1), M@LLER (?), 
Boum, e da uno di noi (J.-P. V.) (4) sul comportamento di masse fluide relativistiche 
rinchiuse in tubi-tempo, studiamo i moti interni corrispondenti a semplici presupposti 
fisici. Introduciamo un nuovo formalismo generale Lagrangiano ed Hamiltoniano in 


funzione delle variabili di Einstein-Kramer e degli angoli relativistici di Eulero eli 
applichiamo ad alcuni casi fisici. 


(*) Traduzione a cura della Redazione. 
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Temporary Capture of Cosmic Ray Particles 
and Their Contribution to the High Intensity Belts. 


R. GALL 


Instituto de Geofisica de la Universidad de México 
Instituto Nacional de la Investigaciôn Cientifica - México 


J. LiFsHitz À (*) 


Instituto Nacional de la Investigacién Cientifica - México 


(ricevuto il 17 Agosto 1959) 


Summary. — The theory of the temporary permanence of charged par- 
ticles in the vicinity of unstable principal periodic orbits in the earth’s 
magnetic field is discussed in relation to the high intensity belts surround- 
ing the earth. The intervals of energies of protons temporarily trapped 
by this mechanism are calculated for various latitudes and distances. 
Several possible sources of charged particles and their contribution to 
the high intensity belts are discussed. The intensity due to the tempo- 
rary capture of primary cosmic ray protons is calculated. The theoretical 
intensity curves are compared with the experimental curves taken aboard 
Pioneer IIT. The secondary albedo protons spectrum is deduced. Theo- 
retical isointensity curves for captured primary and secondary protons 
are plotted for the distances from 2 to 8 earth’s radii. 


1. — Introduction. 


One of the important results of measurements made aboard American and 
Russian satellites and space probes is the discovery of belts of high intensity 
corpuscular radiation surrounding the earth (**). This paper is concerned 


(*) J. Lirsuirz died on May 20th, 1959. The other author takes full responsability 


for the text. 

(1) J. A. VAN ALLEN and L. A. FRANK: Nature, 183, 430 (1959). 

(2) J. A. VAN ALLEN, C. E. Mclzwain and G. H. Lupwic: Journ. Geophys. Res., 
64, 271 (1959). 

(3) S. N. VERNOV, A. E. Coparov, P. V. VAKULOY and Yu. I. Pocaëev: Dokl. 


Akad. Nauk SSSR, 125, 304 (1959). 
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with the temporary trapping of charged particles in the vicinity of unstable 
periodic orbits and their possible contribution to the radiations belts. The 
contribution of both primary and secondary cosmic ray particles is considered. 

The principal periodic orbits in the field of the geomagnetic dipole have 
been studied extensively in connection with the geomagnetic theory of cosmic 
radiation (+1), Their stability has been investigated by several authors and 
the characteristic exponents have been computed (112). For each value of 
Stérmer’s parameter (5) y, > yf (yj = 0.788524), a pair of principal periodic 
orbits exists. These are called the internal and external periodic orbits according 
to their relative positions with respect te the geomagnetic dipole (Figs. 1b 
and 2). For 1<y,<1.313 only one solution exists. The two orbits coincide 
when the constant y, equals 0.788524 and become equatorial for values of 
yi greater than 1.313. The stability of the latter has been studied Ces 


Z Stormers 


> | 


Stormers 


a) 


Fig. 1. — Unstable periodic orbits in the meridian plane. Forbidden (shaded) and 
allowed Stérmer zones for y,>1 and Ya<1. (a) Internal orbit 121; y,=1.1. (6) Ex- 
ternal and internal orbits 6H and 61: =0.83: 


(4) G. STORMER: Zeits. Astrophys., 1, 237 (1930). 
(IG LEMAITRE and M. S. VALLARTA: Phys. Rev., 50, 493 (1936). 
(°) M. S. VALLARTA: An outline of the theory of all 
TT 


Cie owed cones of cosmic radiation, 
University of Toronto Press, Applied M 


athematical Series, no. 3, p. 21 (1938). 


(CO) NTS Si VALLARTA: Theory of Geomagnetic Effects of Cosmic Radiation, Handb. 
d. Phys., 46, (in press). i 


(7) H. Porncarh: Les 
(Paris, 1892). 

(8) G. LEMAITRE and M. S. VALLARTA: 
) G. Lemaitre: Ann. Soc. Scient. de Bruxelles, A 54, 194 (1934) 
(1°) O. GODART: Ann. Soc. Scient. de Bruxelles, 58 27 (1938) 
x A. BaNos, H. URIBE and J. LiFsHITZ: Rev. Mod. Phys., 11, 137 (1939). 
a) 


Méthodes Nouvelles de la Mécanique Céleste, vol. 1 chap. tv 


Ann. Soc. Scient. de Bruxelles, 56, 115 (1936). 


J. Lirsæirz: Journ. Math. and Phys., 21, no. 4 284 (1942) 
R. DE VAGELAERE: Canadian J. of Math., 2, 440 (1950) 
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According to a theorem of Poincaré, the perturbation of the periodic 
orbit (7) is given by 


(1) Sr = exp [Qt] F,(t) + exp [— 21] FL), 


where F, and F, are Fourier series and £ is the characteristic exponent. 


: Sin À 


Internal orbits External orbits 


P=0 


OQ BB nm VJS 


Ore Ol 102s 032904)055 065) 07 08 


Fig. 2. = Principal periodic orbits. Some of the continua of external and internal orbits, 
both stable and unstable, for various values of Stérmer’s parameters y, in the x, 
sin A plane (12). P(x, À) = (2y,)-* exp [2x] — (exp [— a] cos A — 1/cos 7)?. 


A periodic orbit is stable if its Q is imaginary. Such a characteristic ex- 
ponent yields periodic solutions in the neighborhood of the orbit. Any par- 
ticle that remains in the neighborhood of a stable periodic orbit would have 
to originate there; this would be the case with decay protons and electrons 
of albedo neutrons. These orbits will be considered in a future paper. 

Tf Q is real, the orbit is unstable and admits two families of asymptotic 
orbits, one initially towards the orbit and the other initially away from the 
orbit (11). 

A charged particle coming from either infinity or the earth, approaching 
an unstable periodic orbit either along an asymptotic or a nearly asymptotic 
orbit, will remain in the neighborhood of the periodic orbit for a time inversely 
proportional to Q; it will then proceed either towards the earth or towards 
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infinity, depending on its original direction. In this paper we shall refer to 
the temporary trapping of the charged particle in the vicinity of an unstable 
periodic orbit as the temporary capture of the particle by the orbit (™). 


2. — Capture of principal periodic orbits. 


A few of the members of the continuum of principal periodic orbits, both 
stable and unstable, are shown in the a, sin À plane in Fig. 2. In addition 
to these, we also have considered the following orbits: 


Orbit Vi 
(ILS 0.812 
61' 0.815 
| LE 1.070 
| 121! 1.140 
| ZE | 1.150 


The characteristic exponents of the orbits are shown in Fig. 3, where (2/@)? 
in plotted versus y,. Here © is the fundamental frequency of the periodic 


oes a Characteristic exponents (9) and fundamental frequencies (©) of the prin- 
cipal periodic orbits. The unstable orbits lie above the y, axis, the stable below. The 


orbits for y,<1 capture both primary and secondary particles, those for y,>1 only 
] > 1 


secondaries (12), 


(1°) TR. GALL and E. Lirsarrz. 


_ Memoria 2a Asan . 
Fisica (1959); p. 68. a Asamblea Regional, Soc. Mex. de 
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4) 


J) 


Il 


Segretario della S.I.P. 


NORME GENERALI 


Oltre alle lezioni tenute dagli insegnanti ufficiali, seminari e conferenze verranno svolti da 
specialisti su argomenti connessi ai vari Corsi. 


Le iezioni, i seminari, le conferenze, le discussioni saranno normalmente tenuti in inglese 
o in francese. 


L'organizzazione e la direzione dei Corsi sono affidate ai Direttori. 


Chi desiderà partecipare a un Corso dovrà inviare la domanda di ammissione non più tardi 
di 25 giorni prima dell’inizio cel Corso (cioè rispettivamente, a seconda del Corso: entro il 
28 Aprile, il 26 Maggio, il 16 Giugno, il 7 Luglio), corredandola, in modo chiaro e leggibile, 
delle seguenti informazioni: 1) Nome e Cognome; 2) Data e luogo di nascita; 3) Nazionalità; 
4) Indirizzo attuale; 5) Titoli di studio universitari ottenuti e indicazione dell’Università che 
li ha conferiti; 6) Professione attuale; 7) Elenco delle pubblicazioni di Fisica (e per il 3° 
Corso 1960 anche di quelle di Biofisica e Biologia); 8) Quale conoscenza orale e scritta ha 
dell’inglese e del francese; 9) Se verrà a Varenna solo o accompagnato da familiari, e, nel 
caso che venga solo, se è disposto a prendere alloggio in camera in comune con altro allievo. 
Inoltre il richiedente deve allegare alla domanda una lettera di presentazione, rilasciata da un 
professore universitario di Fisica o, eventualmente, nel caso del 3° Corso 1960, da un professore 
universitario di Biofisica o Biologia, il quale assicuri che il richiedente ha interesse per i suoi 
studi a frequentare il Corso indicato e che già possiede adeguata preparazione. 


L'accettazione delle domande è decisa dal Presidente della Società unitamente al Direttore 
del Corso, tenendo conto dei documenti presentati dai vari richiedenti e di un’equa ripartizione 
dei posti tra le varie nazioni cui i richiedenti stessi appartengono, - La ‘comunicazione delle 
deliberazioni prese è inviata agli interessati entro 12 giorni dalla scadenza di presentazione 
delle domande (cioè, rispettivamente: entro il 10 Maggio, il 7 Giugno, il 28 Giugno, il 19 Luglio). 


1 


- Gli allievi sono vivamente pregati di giungere a Varenna il pomeriggio del giorno precedente 
l'inizio del Corso (cioè, rispettivamente, il pomeriggio del 22 Maggio, del 19 Giugno, del 
10 Luglio, del 31 Luglio), e, giunti a Varenna, di presentarsi al nostro incaricato presso il Bar 
dell’« Albergo Stazione », di fronte alla Stazione stessa, per l’assegnazione dell’alloggio, per 
informazioni circa il Corso, ecc. 


Gli allievi saranno alloggiati o nella foresteria della Villa Monastero o in alberghi di Varenna, 
vicino alla Villa, in camere ad 1 o 2 letti. - Nella stessa Villa o negli alberghi sarà organizzata 
la mensa della Scuola. 


- Le somme di cui al n. 10) sono da versarsi, non più tardi di 5 giorni dopo l’inizio del Corso 
(cioè, rispettivamente, entro il 28 Maggio, il 25 Giugno, il 16 Luglio, il 6 Agosto) all’Ammini- 
Strazione della Scuola, a Varenna, in valuta italiana. Borse di studio, in numero molto 
limitato, potranno essere accordate a quegli allievi che ne giustifichino il bisogno. 


La Scuola cercherà, nei limiti delle possibilità locali (che sono però assai limitate causa la 
stagione balneare), di trovare una sistemazione adeguata in alberghi di Varenna per i familiari 
che intendessero accompagnare gli allievi. I familiari possono essere accolti, insieme all’allievo, 
alla mensa della Scuola. - Tutte le spese per alloggio, vitto, ecc., relative ai familiari sono 
da computarsi a parte e totalmente a carico dell’allievo : esse, per trattamento uguale a quello 
fatto all’allievo, si aggirano dalle 2.600 alle 3.600 Lire per giorno e per persona adulta, e 
dovranno essere regolate direttamente dall’allievo con l’albergatore. i 


- Il 1°' Corso 1960 sulle Teorie ergodiche precede un altro Corso dello stesso titolo che si 
terra dal 2 all’ 11 Giugno a Varenna nella Villa Monastero e sarà organizzato dall'Unione 
Matematica Italiana. Sebbene siano indipendenti l’uno dall’altro, questi due Corsi trattano 
il medesimo argomento, il primo sotto l'aspetto fisico e il secondo sotto quello matematico. 


Chi desidera partecipare a entrambi i Corsi i i i i 
: , può rivolgersi per inf i i 
AC A). à EN, g per informazioni al prof. P. Cal 


Milano, 10 Febbraio 1960 


Il Presidente della S.I.F. 


G. C. DALLA NocE 
> G. POLVANI 


I Direttori dei Corsi 
P. CALDIROLA, G, RACAH, N. RASHEVSKY, A. GOZZINI 
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orbits. All the unstable orbits lie above the y, axis, all the stable ones below it. 
Hor.6 lor and 10.0, (Ole) 0. 

As mentioned above, the trapping time in the vicinity of an orbit is in- 
versely proportional to 2; as seen in Fig. 3, the internal orbits are more ef- 
fective in capturing particles than the external ones and contribute more to 
the high intensity belt. 

For the values of the integral of motion y, >1 (Fig. la), the permitted 
zone of motion consists of two parts, one of which is finite (1%), The particles 
in the vicinity of orbits of this finite region (117, 117", 121, 127’) are secondary. 
For y,< 1 (Fig. 1b), there exists only one connected allowed region in which 
both the internal and external orbits can capture particles coming from in- 
finity as well as secondary particles. 

The principal periodic orbits exist between the geomagnetic latitudes + 35°. 
Consequently, the trapping of charged particles by these orbits can only ex- 
plain the contribution to the high intensity belts lying within these latitude 
limits. The contribution of other families of periodic orbits will be discussed 
in a future paper. 

Secondary particles of lower energies that are fed into the high intensity 
belt must originate at higher latitudes. This fact is deduced from the cal- 
culations of the present authors on albedo cones (1*1°). The calculations show 
that particles which originate at low latitudes follow trajectories close to the 
earth, while particles originating at higher latitudes can move away in such 
directions as to be asymptotic or nearly asymptotic to a periodic orbit. 

The possible contributions of the albedo particles (!*) originating at low 
latitudes to the lower layers of the high intensity belt will be studied further 
and discussed in a separate paper. | 


3. — Energies of captured protons. 


For each periodie orbit corresponding to a fixed value of y, a continuum 
of orbits exists (Fig. 4), each of which corresponds to a particular energy. 
The greater the energy, the shorter will be the orbit and the closer it will be 
to the earth’s center (*). An upper limit is thus set to the energy of particles 
that can remain in the vicinity of these orbits, since a periodic orbit of suffi- 
ciently high energy would be tangent to the earth’s atmosphere and, therefore, 


) M. S. VALLARTA: Memoria 2a Asamblea Regional, p. 67 Soc. Mex. de Fisica. 
(14) R. GaLL and J. LIFSHITZ: Phys. Rev., 101, 1821 (1956). 

) R. Gazz and J. LIFSHITZ: Nuovo Cimento, 7, 601 (1958). 
(*) The earth’s center is defined by the central geomagnetic dipole. 


bos 16 - Il Nuovo Cimento. 
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would have no physical meaning. In this paper, only those periodic orbits 
were considered, lying wholly within distances from the earth’s center equal 
to or greater than 11 of the earth’s radius. 


For particles of low energy, no such limit exists. 


However, because of the 


smaller velocity of the particle and its longer trajectory, it is found that, given 


Fig. 4. — Some members of the continuum 
of the 67 orbit in the meridian plane. 
Each orbit corresponds to a different 
proton energy which determines its di- 
stance from the earth’s center. (All di- 
stances are in units of earth’s radii). 


the same number of particles of two 
different energies, the intensity of the 
lower energy particles will be smaller. 
Due to this effect, the contribution 
of such particles to the high inten- 
sity belt becomes vanishingly small 
for sufficiently small energies. 

If we consider any fixed latitude 
A, the distance R in stôrmers of each 
of the unstable periodic orbits from 
the dipole can be obtained by read- 
ing off the values of x (Fig. 2) and 
applying the transformation 


il 
18 = == CP. 
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As one moves towards higher lati- 
tudes, the number of orbits for a 
given À decreases. 

If we now consider the k-th pe- 
riodic unstable orbit at latitude 7 
and at distance r (in em) from the 


dipole, then the energy e, of a particle of rest mass mc? and charge q cap- 


tured by this orbit is given by 


(2) En = VOR) mick — mec, 
where 

2 M2 
(3) CHALET 

(7 à 


M is the earth’s dipole moment. 
For protons and electrons C, = 5.896 72-1038 7-4, 


The energies of protons and electrons captured by 


periodic orbits were computed for latitudes of 0, 5 
? 


the family of unstable 


10, 20 and 30° at each 


of the following 19 distances from the dipole: 14, 14, 14, 13, 2, 24 


Ey Alor) 


3, 3 3, 


DICCI 
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(le) 


1, and 9 earth’s radii. At each À and r, two energy intervals exist, one for 
the internal and the other for the external orbits. Figs. 5a and 5b show these 
intervals for 5° and 30° at all distances. 


100 


Energies of coplured protons in GeV 


DR ASO PRO LIE BO. Cae a 20 
Distance in earth's radii 


a) b) 


Fig. 5. — Energies of protons captured at various distances from the earth’s center by 

the continua of external (upper segment) and internal (lower segment) periodic orbits. 

The energies of the captured secondary protons are indicated by broken lines. 
(Fi So (Ohi dee 30. 


4. — The intensity of captured particles from an arbitrary source. 


Particles from an arbitrary source with a differential energy spectrum D(e) 
that enter the geomagnetic field along certain directions will be temporarily 
captured by the unstable periodic orbits and contribute to the high intensity 
belts. 

In order to evaluate the total intensity I at a distance r and latitude 
A due to the source D(e), the intensity per unit energy, iz, Of particles of 
energy e, captured by the k-th orbit for these values of À and r, must first 
be computed. 

The intensity i, depends on the following factors. 


a) The number of oscillations the particle makes in the vicinity of the 
orbit. This number is proportional to w/Q. We shall refer to it as the weight 
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of the orbit and define a relative weight, px, as the ratio of the weight of the k-th 
orbit to that of the 1H orbit (Fig. 3) 


(4) Pr = (0/2) , 
© (@/2)yz 
b) The volume defined as the vicinity of the orbit. As a first approx- 
imation we shall assume that this vicinity is a portion of a spherical shell. 
Tf the k-th orbit is at a distance r, the intensity of 7, is then inversely propor- 
tional to (r/r,)?, where r, is the earth’s radius. 


c) The intensity of particles of energy e, at the source. 


For particles of any given energy e,, the value of i. at a distance r and 
latitude À is therefore given by 


- TOI 
(5) be Co Pa Tip ya Plex) : 
where 0, is a proportionality constant. 

The total intensity J at a given latitude À and distance r of particles of 
energy spectrum D(e), captured by internal and external unstable periodic 
orbits existing at this latitude, is evaluated by integrating the areas under 
the respective i, vs. e, curves. 

The intensity J was actually computed for primary and secondary cosmic 
rays by an IBM Model 650 calculating machine. The computations were 
carried out for 5 latitudes and 19 distances, and the areas corrected for the 


asymptotic behavior of the i, vs. e, curves at the ends, (6*Z, 12*7, and 10L 
orbits). 


5. — Sources contributing to the high intensity belt. 


Sources that might contribute to the high intensity belts have been dis- 
cussed by various authors (11518), Solar plasma and secondary cosmic rays, 
including albedo protons and neutrons, as well as their decay products, have 
been suggested as possible sources. It has also been suggested that the first 
maximum in the intensity curve is caused by the cosmic ray albedo and the 
second (more distant from the earth) by solar plasma. 


(19) S. N. VERNOY et al.: Special Lecture, 5th Gen. Assembly of OSAGI (Moscow, 
July 30- August 9, 1958). | 

(7) S. F. SINGER: Phys. Rev. Lett., 1, 171. 181 (1958). 

(77) S. F. SINGER: 2nd Symposium on 8 
Texas, November 1958. 


(8) P. J. KELLOGG: Nuovo Cimento, 11, 48 (1959). 


pace Physics and Medicine, San Antonio, 
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Particles sent by any source into the geomagnetic field will be temporarily 
captured and cause an increase in the intensity of the radiation belts. Ac- 
cording to our theory, the distance at which this increase can be detected de- 
pends upon the spectrum of the source. If capture by unstable periodic orbits 
were the only operating mechanism, then, for example, a source of protons 
of energies less than 50 MeV could only be detected at distances greater than 
six earth’s radii (Fig. 5). 

The capture mechanism discussed here is just as effective in capturing 
primary as secondary particles. The fact that the higher energy particles are 
captured closer to the earth suggests that the first maximum of the intensity 
of the radiation belt is caused by primary cosmic ray particles, while the 
second maximum is due to all secondary cosmic ray particles. 

In this paper, we have only considered two sources, namely, primary pro- 
tons and albedo protons of cosmic rays. 

The intensity i, of the primary protons captured by the k’th orbit at a 
given À, r, was computed at all lati- 
tudes and distances previously men- 


tioned. Equation (5) was used, in 
which D(e) was replaced by the pri- 
mary differential spectrum (1°) N(e), 
where 


(6) N(e)de=0.048e-4[ (140.09 e5)]-*de 


for e > 0.500 GeV 
—( for € < 0.500 GeV. 


1500 


0.500 GeV was chosen as a reasonable 
cut-off for periods when the solar 
activity is not too high (221). The 1000 
total intensity I of the captured pri- 
mary protons was then calculated for 
all latitudes and distances by the 


Intensity 


500 


Fig. 6. — The intensity curves of captured 0 a . ed 
primary protons. Intensity vs. distance 1 2: MISE 5 6 7 
at various latitudes. Distance in earths radii 


G) H. V. Neuer, in J. G. Wirson ed.: Progress in Cosmic Ray Physics, p. 300 
(New York-London, 1952). 


(2) H. V. Newer: Phys. Rev., 103, 228 (1956). 
(2) P. Meyer and F. A. SIMPSON: Suppl. Nuovo Cimento, 8, 233 (1958). 
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method described in the previous section. Fig. 6 shows intensity decreases 
with the distance but passes through a maximum at approximately 1.7 7, 
the position and intensity of this maximum being a function of the latitude. 
As the atlitude increases, the 
intensity first increases slightly 
and then decreases conside- 
rably. 

A study of the energies 
(Fig. 5) of the protons, captured 
at distances between 2.5 and 
4.57, suggested the idea that 
the second maximum of the 


50.000 


10.000 


5,000 


1000+ 
aa experimental intensity curves 

2 is due to albedo particles. The 

2 experimental curves could then 
È n be reproduced by the additive 

8 capture of primary and secon- 
coe darv cosmic rave 

ary cosmic rays. An attempt 
was made to deduce a distri- 
5 bution of captured secondary 
i (2/03) 14 506078102300 5607 8 “protons by “comparineweine 
oes in earth's radi theoretical intensity curves 

a) 

i x x for captured primary protons 
Fig. 7. — Experimental and theoretical intensity oN P P ÿ P : 
curves for the (a) out-bound and (6) in-bound (Figs. 7a,b, curve 3 and 4) with 
leg of Pioneer III’s trajectory. Curves 1 and 2 the experimental Intensity (*) 
show À measured aboard Pioneer III; curves 3 curves (1) obtained aboard Pio- 


and 4 show the calculated value of R for cap- 
tured primary protons, and curves 5 and 6 for 
captured primary and secondary protons. 


neer III (Figs. 7a,b, curves 1 
and 2). 
Curves 3 and 4 were obtained 


by plotting the intensity of cap- 
tured primary protons at the geomagnetic latitudes 7 and distances r of the 


outgoing and incoming trajectories of Pioneer TIT. Differences AZ between the 
experimental and the theoretical curves were taken and set equal to the total 
intensity of protons of the secondary source whose distribution was to be found. 
However, no distribution seems to exist which would reproduce both curves 
obtained by space probes, if capture by periodic orbits were the only mechanism 


operating. This negative result indicates that other mechanisms must also 
contribute to the high intensity belts. 


(*) In actual fact not the Intensity 


I, but the counting rate R, w ers 
distance where R— TEINS È Se 
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A distribution function of the form 59(1+e)-5 was tested and seems to 
give satisfactory results when limited to the range 0.45 << 2.20 GeV. If 


Lili 


the cut-off is not applied at the 
lower limit, the intensity becomes 
too high at large distances from the 
earth’s center. The higher energy 
cut-off brings the intensity at shorter 
distances down to the observed level. 
The energy spectrum finally used 
for the secondary protons was the 
above function modified at the two 
ends so that it has a sharp peak at 
0.52 GeV which tends to zero towards 
the lower cut-off value and has a tail 
extending to 2.20 GeV. The intensity 
due to this calculated source was 
computed for 0, 5, 10, 20 and 30° | 
of geomagnetic latitude at all 19 


dist | Fig. 8. — Iso-intensity curves in the geo- 
IAS: magnetic meridian plane. The intensity is 
The intensities due to this distri- due to both primary and secondary pro- 


bution of secondary protons were tons temporarily trapped by unstable pe- 


evaluated for À and r corresponding riodic orbits. Only curves for distances 
larger than 2r; were drawn. (All distances 


to the outgoing and incoming legs à È x 
ese SRO. are in units of earth’s radii). 


of Pioneer III’s trajectory and added 
to the intensities of captured primary 
protons at these 4, r values. The resulting curves are plotted in Figs. Ta,b 
(curves 5 and 6). 
Isointensity curves in the geomagnetic meridian plane were also plotted 


and are shown in Fig. 8. 


6. — Discussion. 


Principal periodic orbits and the capture mechanism suggested here are 
only valid if the dipole approximation to the geomagnetic field is also valid. 

The capture mechanism discussed here can only partly account for the 
existence of the radiation belts. We believe that several mechanisms are re- 
sponsible for the radiation belts (1°). The mirror points mechanism (171$) 
could probably explain the capture of lower energy particles, for which the 
approximation of the constancy of the magnetic moment holds. 

Solar plasma was not considered among the sources discussed in this paper. 
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The injection of solar plasma is accompanied by a perturbation of the geo- 
magnetic field which, especially at large distances, causes a deviation from the 
dipole field. When a small enough perturbation exists, capture of the solar 
plasma particles by unstable periodic orbits would occur. Depending on the 
energy spectrum of the plasma, an increase in the intensity of the radiation 
belt would occur at corresponding distances, and for high enough energy of 
plasma particles the second maximum would be affected. The characteristics 
of the second maximum would be expected to vary more with time than those 
of the first maximum. 

There are discrepancies between the experimental curves and those cal- 
culated on the basis of the temporary capture by unstable periodic orbits, 
which may be partially corrected by using a more adequate energy spectrum 
of secondary protons. An explanation of the remaining basic discrepancies 
will have to be postponed until some idea of the relative importance of other 
mechanisms becomes available. 

The isointensity curves shown in Fig. 8 depend strongly on the distribution 
chosen for secondary protons and are drawn for distances larger than 2r,. 
They will be extended to the region below 27, when the contribution of al- 
bedo trajectories to the lower layers of the intensity belts will be studied. 

The proportionality constant C, is the only free parameter in the present 
theory and was chosen to be 10° in order to bring the calculated value of the 
intensity within the range of measured values. This constant is associated 
with the influx and efflux of particles into and out of the radiation belts. 
A full discussion of this constant must also be postponed until quantitative 
data are available for other trapping mechanisms. 

A more detailed discussion also requires a knowledge of any variations in 
the characteristics of the high intensity belts with time. 

The soundness of the theory of temporary capture by unstable periodic 
orbits could be checked experimentally by determining the energy distribution 
of the particles composing the first and second radiation belts. 


The authors wish to thank Professor MANUEL SANDOVAL VALLARTA for 
valuable discussions on the subject. We would also like to thank Mrs. Luca 
SANCHEZ DE CAMACHO for help with the calculations, Mr. Luis OTTALENGO 
of the Centro Electrénico de Cäleulo of the National University, for compu- 


tations on the IBM calculating machine, and Mrs. Joan KRUGER DE HAUFF 
for careful reading of the English text. 
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RIASSUNTO (*) 


Si discute la teoria della permanenza temporanea di particelle cariche in prossi- 
mità di orbite periodiche principali instabili nel campo magnetico terrestre in relazione 
alle fasce di alta intensità che circondano la terra. Si calcolano per varie latitudini e 
distanze gli intervalli di energia di protoni temporaneamente catturati da tale mec- 
canismo. Si esaminano numerose sorgenti possibili di particelle cariche ed il loro con- 
tributo alle fasce di alta intensità. Si calcola l’intensità dovuta alla temporanea cattura 
di protoni della radiazione cosmica primaria. Si confrontano le curve teoriche dell’inten- 
sità con le curve sperimentali rilevate dal Pioneer III. Si deduce lo spettro dei pro- 
toni secondari d’albedo. Si tracciano le curve teoriche passanti per i punti di uguale 
intensità a distanze comprese fra 2 e 8 raggi terrestri. 


(*) Traduzione a cura della Redazione. 
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IL NUOVO CIMENTO Viole XV Nee 16 Gennaio 1960 


Particle Production in 6.2 GeV p-p Collisions 
Treated by a Statistical Model. 


R. HAGEDORN 
CERN - Geneva 


(ricevuto il 26 Settembre 1959) 


Summary. — The Fermi theory of particle production is used in a rigorous 
and refined form in order to caleulate spectra and mean particle numbers 
for all produced particles (i.e. 97, x, 97 and strange particles) in 6.2 GeV 
(kin. lab.) p-p collisions. 


1. — The Fermi model. 


The probability for a certain final state b is in the Fermi model (1): 


28; +1)5 ; 
(1) P,= Weby(L) (II LITI gerot(E, Mie Mg P= 0), 


RE di 


here W.;,,(7) is the number of linearly independent isospin eigenfunctions of 
total isospin 7 and fixed but arbitrary 7, which can be built up with « par- 
ticles of isospin 3, 8 particles of isospin 1 and y particles of isospin 3 (2). 
The next factor counts the multiplicities of ordinary spin where o; is the 
number of particles of spin S,, the product [I N;! stays for the indistinguish- 
ability of the particles. {2"-1. (27) ig the so called interaction volume 
and 9"(H, m,... mn) is essentially the phase space density of the final state: 


(2) O*(E, m, ... Ma; P) far, ba dp,Ô(E — DV mè) (> p:—p) r 


(2) F. CERULUS and R. HAGEDORN: CERN Report 59-3. 


| (?) Y. YEIVIN and A. DE-SHALIT: Nuovo Cimento, 1. 1147 (1955); V. S. Bara: 
SENKOV and B. M. BaRBASEV: Suppl. Nuovo Cimento, 7. 19 (1957) 


PARTICLE PRODUCTION IN 6.2 GeV p-p COLLISIONS ETC. 247 


Also the momentum spectrum of a particle in the C.M.S. of the collision (CMS) 
is given essentially by o*, e.g. for the particle m,: 


(3) w,(p) dp ~ 4xp?dp - 0*(E —V m2 + p°, Ma ee Mn, P = — Pa) - 


We have a Monte Carlo method (*) which allows us to calculate such integrals 
with any prescribed accuracy and we choose the accuracy such that the results 
given below are 


exact within less than + 1% for the particle numbers , 
» » > » +3% for the spectra. 
In particular, the influence of small differences in the masses, between e.g. x 
and A, shows up clearly far outside the statistical error. 
We have considered all reactions at this energy, which contribute to the 
observable final states and give a list of them below (Table I, Appendix I). 
For all these reactions P, has been calculated with the formula (1) and then 


the kinetic energy spectra w(e), of the different kinds of particles which come 
out from the Monte Carlo calculation in normalized form 


(4) fe) de 1 


are multiplied by P, and by the number of particles n® which are represented 
by w(e) in the reaction considered: 


{ à labels the kind of particle 


(5) WP (e) = Pn 100 (€) | | 
b ) ) » » reaction. 


Finally, they are added up for each kind of particle: 
(6) Wile) = 3 Pin we) , 
b 


which is the not-normalized total spectrum of particle 7. We have then, since 
w(e) is normalized to unity: 


pro de> Bony = <M) D 
b b 


(8) F. CeruLus and R. HAGEDORN: Suppl. Nuovo Cimento, 9, 646 (1958). 


~ 


a 
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Here <n;> is the mean absolute particle number per collision, averaged over 
all reactions with proper weight factors. 

Finally, the normalized total spectrum and the mean kinetic energy are 
calculated for each kind of particle. 


2. — Refinements to the simple Fermi model. 


We have, however, made a few refinements which seem to be essential: 


a 


1) We include the formation of the isobar 97* (T=3, S— 3 mass 1.316) 
which afterwards decays into 97+ 7 (1). 


2) We include final state annihilation of 97-97-pairs by saying that they 
will annihilate, if they come out with a small momentum relative to each other 
and thus stay for some time near to each other. We assumed that the pro- 
bability for annihilation is given by 


[probability for a certain relative momentum, p] times [probability for an- 


Der 


nihilation when the relative momentum is p,]. 


The first factor is calculated by means of the phase-space by simulating a 
M-H-pair with relative momentum p, by one particle of mass 


Otota(Pr) 
calculated with the help of a simple model with complex potential, which 
represents fairly well the data for pp-scattering up to 300 MeV (*). 
The number of annihilating pairs at a given reaction b turns then out to 
be approximately (see Appendix TT) 


AN~N È = (1 | 7 


0 * 
with 
Do 


eal 
Ao* = const 220. p2* : 
0 ast Pr'0*(E, mi... my, m*(p,)) Ap, . 


Ototal 


0 


(4) 8. Z. BELENKI et al.: Usp. Fiz. Nauk. 62 (2), 1 (1957). 
() Private communication by F. CERULUS (not published) 
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This has then to be subtracted from the number of pairs which the primary 
interaction yields. In à similar way the 7-spectrum is changed, in particular 
the low energy part of it is reduced as it is to be expected from this model of 
final state annihilation. We give below the original 97-spectrum, the AGT- 
spectrum (that is: the spectrum of the annihilating AT) and the observable 
spectrum 97 — AI. 

In the yield of 97 this final state annihilation gives à reduction by a factor 
0.26, which is quite important. 


3) We consider a third rather trivial final state interaction by calculating 


the spectra of y and A° from the decays: 


Too Ss 2Y 


Mor AO + Y 
and by adding up all such y-spectra and by adding this A° to those produced 
directly. 


4) The interaction volume is split into three factors: 
Pil a 2. (je. ON 
Q Es 92° £2 3°? 


where now n, = number of 7 + QT +9 + 7", 
nia = number of 2D+A+4, 


Ms = number of K + K. 


We calculated now everything in several versions, namely (*) 


43 ; 
(OT 0a i |; Qx variable, 
3; AT 45 . 
(i). O gino Qy = Qe variable . 


Assumption (i) seems to be related to global symmetry whereas assumption (ii) 
treats the hyperons and K-mesons the same way. 

If one considers © as an expression of the coupling strength, then one 
should expect Q, to be 10--15 times smaller than OR 


(*) We use simultaneously the notation 


Cra as Q, = Qy (hyperons) , Q, = Ox. 


{7 
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Tf, on the other hand, one interprets it as expressing ranges of interactions, 
then Q, should be the same expression as Q_ but with 4, replacing À. Then 
it should be about 40 times smaller than _. 

The meaning of © being obscure, we varied 2, from 2, = 2, down to 
Q, = + ©. and that in both versions (i) and (ii). 

The curves given below show then the effect of this variation. 

Tt is interesting to see that there are certain quantities which depend 
strongly neither on Q, nor on the choice of (i) or (ii), e.g. the total y, 97, x- 
production, whereas others depend very strongly, changing by orders of mag- 
nitude. 

The normalized spectra are hardly influenced and can be said to be in- 
dependent of the version taken. 

Further explanations will be found in the figure captions. 


8. — Angular distribution. 


In this kind of model nothing is known concerning angular distribution; 
in fact, one calculates only an average spectrum over all angles. 


4. — Inelasticity. 


In the Fermi model, as used here, it is assumed that the entire available 
energy is distributed statistically between all possible end products. This can 
of course only be true if the two nucleons suffer a «central » collision. The 
present calculations should therefore strictly speaking be compared only with 
inelastic events with high inelasticity. 


5. — Coneluding remarks. 


The present calculation shows that the Fermi 


model, taken seriously, is 
everything else but a simple model. 


It was possible to perform this calculation 
only with the help of an electronic computer (CERN’s Ferranti Mereury) (5). 

We got the impression that the model may very well be able to reproduce 
many experimental facts—exeept angular distribution—and that the discre- 


(?) Cfr. R. HAGEDORN: Report CERN 59-25, 


cc 
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or 
i 


pancies encountered so far (up to x 50 GeV) lie not in the model but in the 
fact that it has not been used correctly; may that concern very bad approx- 
imations for the phase space or the omission of the other factors, which easily 
become of the order of 1/100 or may it come from considering only a small 
subset of all possible and contributing reactions at a given energy. 

Certainly the final state interactions play an important role also. 

Recently x-r-interactions with a resonant state (which would appear in 
the present work as a particle of about 4 pion masses, isospin 1 and angular 
momentum 1) have been discussed (*). In fact, the inclusion cf such an inter- 
action in the treatment of nucleon pair annihilation with the statistical model (7) 
leads to surprisingly good agreement with experiments. 

Consequently, it seems that this interaction should be included also here. 
We have done such calculations, with the definite result that it leads to dis- 
agreement with experiments. The total pion number comes out to be ap- 
proximately 4, instead of 2.8 as here and in the experiment. Other authors (°) 
come to similar conclusions in the case of x-p inelastic collisions. 

It is very hard to see how one could find a reason, which rules out r-r- 
interactions in nucleon-nucleon collisions, but makes it important in annihi- 
lation. A possible way out would be to relate the x-x isobar to the nucleon 
isobar. Then one hag to include either the nucleon isobar or the pion isobar, 
but not both of them and obviously the nucleon isobar should be taken in 
nucleon-nucleon or pion-nucleon collisions and the pion isobar in annihilation 
processes. 

There might, be another explanation, which is much less ambitious: the 
experimental results include non-central collisions and these, for purely geo- 
metrical reasons are the most frequent ones. The present theory, however, 
applies to central collisions only. One could imagine, that for large impact 
parameters, the number of pions produced goes down, since only part of the 
primary kinetic energy is available for particle production. If we now would 
include the pion isobar (which gives 4 pions per collision) and simultaneously 
could include the effect of non-central collisions, then these two might some- 
how compensate and the 2.8 pions might result again. 

To test this explanation, further considerations on the mechanism of non- 
central collisions are required, together with a large amount of additional cal- 
culations. All one can say at the moment is that if the impact parameter is 
not larger than «& 1/24, equilibrium could be established in the whole inter- 
action region, within the collision time (for the energy considered here). Col- 


(6) W. Fraser and J. FuLco: Phys. Rev. Metti, 2, 360 (1959), 
(7) F. Cervus: Nuovo Cimento, 14, 827 (1959). 
(8) V. S. Baraëenxov and V. M. Mattspv: Dubno Report R-350 (May 1959). 
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isions with an impact paranieter > 1/24 should be considered as peripheral (*). 
They might lead to lower multiplicities (and also be mainly responsible for 
angular anisotropy). At the present moment the influence of the x-r isobar 
is therefore still an open question. 

The results, together with further explanations, are given in Appendix I. 


The author is very grateful to Dr. CERULUS (CERN) for many discussions 
covering almost every part of the present work, in particular the effect ot 
final state annihilation and the x-x isobar. In the spirit of our collaboration 
in earlier publications Dr. CERULUS even took over the tedious job of pre- 
paring some of the data tapes and watching the runs of some Monte Carlo 
calculations on the computer. Finally he read critically the manuscript. 


ASP DE NED Xe 


The initial state is always pp, the final states are listed below together 
with their probabilities in percent of the total cross-section (which cannot 
be calculated by this theory). We give two versions: 


(CA OR KO Oo 0% (« global symmetry »), 


al 
3.5 


(b) Gr = Dr (06, 


(The value 1/8.5 is purely accidental — it just happened that this was one 
of the values for which everything was calculated in order to draw the curves 
of Figs. 13-21. This value, however, is the nearest to 1/10, which we would 
consider reasonable. Anyway, the difference is unimportant). 

We neglected all those reactions which, within a group of similar reactions, 
are very improbable, with respect to the others of the same group, with some 
exceptions. These exceptions have no influence at all, we list them for curi- 
osity only. 


Notation. — 1.094 


5 means 1.09-10=5 ete. 


() See a forthcoming paper on results at 25 GeV by the author, to be published 
in Nuovo Cimento, Vol. 15, n. 3. 
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GENERAL INFORMATION 


A) - In addition to the Courses which are listed above, seminars on questions related to the 
basic subjects of the different Courses will be held by visiting scientists. 


B) - The lectures, seminars, conferences and discussions will generally be held in English or 
French. 


The organization and direction of the Courses is entrusted to the Directors. 


An application should arrive not later than 25 days before the beginning of the Course 
(i.e. resp., not later than April 28th, May 26th, June 16th, July 7th) displaying the following 
information, which should be set out clearly and legibly: 1) Christian name and surname; 
2) Date and place of birth; 3) Nationality; 4) Present address; 5) Degrees and other 
qualifications obtained, with name of University in each case; 6) Present professional activity ; 
7) List of publications in Physics (and for the 3rd Course 1960 also in Biophysics and 
| Biology); 8) Standard of knowledge of English and French - written and spoken; 9) Whether 

È. intending to stay at Varenna unaccompanied or with members of family and, in the former 
| n case, whether willing to share a room with another student. The applicant should send, 
together with his application, a note of reference from a University Professor of Physics 
or, eventually, in the case of the 3rd Course 1960, from a University Professor of Biophysics or 
Biology, testifying the applicant’s interest and preparation for the activities of the School. 


D) 


E) - Each application will be considered by the President of the Italian Physical Society and 
by the Director of the School on the basis of the information submitted, with regard also 
to a fair distribution of the places available among students of various nations. The decisions 

} on the admittance to the School will be made known to the applicants within 12 days 

> from the date fixed for the arrival of applications (i.e. resp., May 10th, June 7th, June 28th, 
July 19th). 


F) - Students will oblige in arriving at Varenna in the afternoon of the day preceding the opening 
+ TT of the Course (i.e. resp., in the afternoon of May 22nd, June 19th, July 10th, July 31st) 
and calling at our delegate at the Bar of the « Albergo Stazione », in front of the Station 
itself, for lodgings, and to receive information and papers concerning the Course. 


G) - Accomodation for the students will be provided in rooms with 1 or 2 beds either in the 
guest quarters of the Villa Monastero or in hotels at Varenna. Meals will be had at the 
Villa itself or at a hotel in Varenna. 


H) - Fees as per 10) should be paid not later than 5 days after the beginning of the Course 
(i.e. before May 28th, June 25th, July 16th, August 6th respectively) to the management 
of the School, at Varenna, in Italian currency. A very limited number of scholarships 


ù may be granted to students whose economic conditions might otherwise prevent them from 
: attending the School. 


Se 


The School will do everything possible to find suitable accomodation in local hotels for 
members of families accompanying the students. It should be noted, however, that, in view 
of the holiday season, local possibilities are limited. Members of students’ families may 
avail themselves of the catering arrangements being organized for the School. All expenditure 
involved in the hotel accomodation, board etc. for students’ relatives will be payable separately. 
These expenses will be, according to the hotel, from 2.600 to 3.600 Lire per person pér day 
and are to be settled directly with the hotel management. 
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J) - The lst Course 1960 on Ergodic Theories precedes the other one which, under the same 


title, is to igs held from 2nd to 11th June in Varenna at the Villa Monastero, and is organized 
by inks Unione Matematica Italiana. Though independent one from the other, the two Courses 
deal in the same subject, the former under the physical aspect, the latter under the 


mathematical one. People who want to attend both these Courses should ask for information 
prof. P. Caldirola (see 4) and 5)). 
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VIETA 


Milan, 10th February 1960 


The Secretary of the S.I.F. 


The Presid sw 
ee e President of the 8.1. F 


G. POLVANI 
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The Directors of the Courses 
P. CALDIROLA, G. RACAH, N. RASHRVSKy, A. GOZZINI 
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Taste I. — Reactions considered. 
Prob. in % 
Reactions | ir n 2 
| (x) (B) 
| | 
297 (elastic, incoherent) (ZA AI aliis 
2704 £59," = 0 488 eee 0) 
297 +27 1.09, +1 | 116, +14 
290+43x | 5.51, + 0 5.80, + 0 
297 +4r Sos cl SIL ET 
297+5r | 310, — 2 | BI, — 2 
TI*+N > 294 | SU NE RL COR oe Al 
q*+Qir > 29420 EM 08 ER 11534 “teal 
WL + 2r > 2901+ 3x N Ido A 2.08: E toa 
M+ 4 3n > 290447 8.70, + 0 9.24, + 0 
WELW +4n > 291-450 I I MI 1 
QT* + QT + 57 > 291 +67 LOL = 2) ANS OS BESS 
29C* > 297+2n 4.70, — 1 | 5.00, — 1 
2*+r >297+3x i 1064 4 FI ee ree, ah 
29T# +27 > 297 +45 len) oe! (SAN 1720 ee Da 
297* + 3x > 297457 | 2.55, + 0 ely, + 0 
2qT* +47 > 297+67 | 1.05, HUIT EEE 
| F 
290*+5r > 29+ 78 RES ME 1.20, — 3 
| | 
397+97 (without final state annihil.) 0 461 03 1:00, "0209 
| 
Ç | .43, — | 4.28, — 2 
A+Q9Q7+K | 3.43 8 
MOT ER MERE: 2101 1.08, —1 
ACCO Kt or | 2.81, — I Bo, = 7 
A+T+K+3r (ME tera 2 DIOR 
| 
A+GT*+K >+A+974+K+7 | 4229; — l 5.26, — 2 | 
A+9M4+K4+n > A+M+K+27 RS CS ARE A =e 
A+90*+K4 2x >A+90+K+3n | 1.04, — 1 1520; — 2 
| 
A+QT+EK +37 -A+MN+K+47x | 2.38, — 3 2.96, — 4 
| 


17 - Il Nuovo Cimenio. 
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TABLE I (continwed). 
Prop ala, OF 
Reactions = Es 
(x) (B) 

yO K BAD tT 6:83. wee 
S4+97+K+2n 1.30 0200; 1°6 [aan 
Z>+97+K+2r SI SMS] AD > 
Z+9(+K+3r DLL ee Dal Whe we eS 
DOT K ed Q1-LR LT 602.1 RI E 
DIS IRA, OA, — Tl oe — 1 
DOTE LKR ET SOT K-+35 (TG RE L'A 
X+9(*+K43r->Z+9+K+47r DAT, — à 3.08, = di 
E+97+2K 1:30 ee Wie. ES 
E+97+2K+xr 1-894 seo ives = Yh 
E+QT+2K+2r 1:09, 05 ie in 
E+97+2K+3r 3.64, == 9 4.54 010 
E+97*+2K > E+97+2K+x DITA SAS NET 
E+Q7*+2K Er >E+97+2K+2r 304 ea ST 22.6 
E+97*#+2K+2rx > E+97+2K+3xr LOG MENTO Cars tile 
20 -+K+K 6.60 er Ol en 
2977+ K+K+x 3145 ae 3°66, 22s? 
277+K+K+ 27 ils — 3 2504 e TR) 
i 290+K+K+3r QU HF 2.64 == 5 
TOUT T K+ K —>297+K+K4+x QU = BY S.65 22 
N*+N+K+K+x +297+K+K+2r I, = © Tua 307 
T*+MLK+K+27>27+K+K+3r 2.98, Da 316) eed 
T*+9+K+K+3rx > 2974+ K4K 445 SH yi 304 nT 
200*+K+K -2797+K+K+2x IN = INS, > 
270*+K+K+r >277+K+K+3r GAS RE I 
290#+K+K+2x +297+K+K44x 3143 = 7 3.65 N 
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50 MeV] 
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| Total N 


<Ekin> 


SS SSS Se Ekin in MeV 
200 400 600 800 1000 
Fig. 1. — Total nucleon spectrum. 
er SOMeV J 


cn? 


TT From isobar decay 


<E Kin? 
L = + T 1——Ekin/N MeV 
200 400 600 800 1000 
_— The x-mesons coming from the decaying 
isobar. 
aS 


NE) Per 50 MeV] 
<n> 
0.25 


020 


Totaln(From direct production 
and isobar decay) 


0.10 5 


0.054 


<E xin? 
T 7 T T t T T TÈ 
200 400 600 800 1000 
Fig. 2. — Total x-spectrum. The peak at low energy 
is mainly due to the decaying 9* (compare 
Fig. 3). The r+ and x spectrum will look similar 
except for the high-energy tail, where the charge- 
state probabilities are different for m+, 7° and x. 
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Total y (From mand © decay) 


0.05 
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Fig. 4. The total y-spectrum coming from 7° 2: 
and A++. Other y-processes are not con 
tained herein. 
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n(e) [per 50MeV 
n(e) [Per 50 MeV] N(E) DERS EME 
<n 
0.254 des 
| 
020: co 
J 
015, i 015 | 
Total AXFrom direct production | 


andËdecay) 


0.05 > 005 
<E kin? 
si ioni ——— EkiniN MeV i ! re order ee ON 
200 400 600 800 1000 200 400 600 800 1000 
Fig. 5. — Total A-spectrum, including the decay Fig. 6. — Total X-spectrum. The spectra fc 
of 4 of the produced E. and © may be somewhat different, in parti 
for the higher energies. 
nie) [Per SQMev] NE) [Per 50 MeV] 
<n> 
025: 
0.204 
0.15 4 
| Total K* (From associated + 


pair production) 


—Exinin MeV = i Piste RATE 
200 400 600 800 1000° 200 400 600 800 1000 


Fig. 8. — Total K-spectrum. To this cont: 
the K-mesons from associated productio1 
those from pair production. 


Fig. 7. — The total E-spectrum. 
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On 


Explanations for Figs. 1-12 (spectra). — All spectra represent kinetic energies 
in the C.M. system of the colliding protons. They are normalized to unity. 
The quantity represented is the ratio: 


number of particles n(e) per energy interval 


mean total number (n> of particles of that kind ‘ 


n(e) [Per 50 MeV] 
<n? 


025! 


0.05 


<E kin? 


Exinin MeV 


mad 


AE —— 

200 400 600 800 1000 

Fig. 9. — K-spectrum. Only KK-pair production contributes. Note the different shape 
as compared with Fig. 8. 


The energy interval was chosen to be 50 MeV except for the antinucleons, 
where it is 1 MeV. 
The points are the resulting values from the above described weighted 


ne) [Per 1MeV ] 
<n? The A 
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ie MO, — Observable anti-nucleon spectrum, assuming final state annihilation. 
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superposition of the Monte Carlo spectra. Some spectra are represented very 
smoothly, namely those to which many of the reactions listed in Table ti 
contribute (e.g. r-spectra). Some others show still large fluctuations (e.g. 97). 


NE) [Per IMev] 
<n> 
I > 
0.025, Original N 


0.0207 ° 


0.015] 


0.010 5 
0005! 
; ? li " 1 r T | T= € kinin MeV 
20 40 60 80 100 120 
Fig. 11. — Original anti-nucleon spectrum, without final state annihilation. 


The solid curves are partly drawn by hand and partly found by a least square 
fit by a Fourier series of sin ve, v=1...5. The curves are statistically more 
accurate than the points and their error is probably of the order of — or less 


n(e) [Per 1 Mev] 


<n? 


0.0255 
| AN (Annihilating N ) 


0020 
0015- 


i y ae T_T Ex, in MeV 
20 40 60 80 100 12008 
Fig. 12. — Spectrum of those anti-nucleons, which have a chance to annihilate in final 
state. This spectrum, with a proper weight factor is subtracted from the original 


9 -spectrum, giving that of Fig. 10. Note that all three spectra are normalized to unity. 


‘ BREA 0/ nog rT AG : 
than about 3%, for some curves much better. The amount of scattering 


of the drawn in points gives a qualitative feeling for the goodness of fit of 
the curve. 


Explanations for Figs. 13-21. — Here the mean v 
numbers per collision and some ratios of 
kind of particle. 
charges. 


alues of absolute particle 
on such numbers are shown for each 
We made no attempt to divide them up for the different 
Roughly, but not exactly, the numbers of Tid ave mL eG peels | 


on 
Ual 
N 
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and analogously for Y, 97. But for cases in which very few particles are 
produced, this is certainly wrong. If in experimental analysis one classifies 
events according to the number of prongs, then one has to be very careful. 


a 
| 
| 
Die 
(o) 
[ Sk Lal a 
Ri | 10 | 
I Se La id SIE Di SAP ea du SS ug, 
0a? ee 2 Sey dE 5 10 
Fig. 13. — Mean number of x and 90 per collision. 


(See F. CERULUS and R. HAGEDORN: CERN 59-3, where the charge distribution 
at another energy is carried through. A calculation to do the same in the 
present case is in progress). 
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Fig. 14. — Mean number of 97 per collision. 
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x 10° 
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155 

150| 
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Qk LES da 
EST e SE E EA ne Pe Tee -Qk 
EME 10e SO SANO 


Fig. 15. — Ratio of n-/p. 


In all figures Qx is the varying variable and the curves show the depend- 
ence of the particle numbers on Qx. In all cases where the assumptions 


(i) MEO 
and 
(ii) Oy le, 
| 
1.0 <n?y 


O 
(o) 


ERGE EE LL I 


Oo 
(co) 


| | 
eet | ae etero 29) 
Jos (OURS SOS. pn 
Ho Gn — Moi y number from x° and Y° decay per collision. 
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10° <P 507215 <N? Ndirect 


10 


(jor 2 OZ G Wo? 2 5 10 


Fig. 17. Total X and A number per collision: the total number of £2*+2°+ 27 equals 
approximately 3 of the number of directly produced A°. If the masses of & and A 
were equal, one would expect a factor 3. but the reduction in phase space of the 
states with ©, because of its higher mass, gives almost exactly a factor 1.5. The curves 


represent therefore simultaneously the mean number of Et and the mean number 
of AS. 4 +20, which equals the total number of observed A. 


show different results, two curves are given. When only one curve is drawn, 
it means that both assumptions lead to the same result. 


o) 
i 2 er 2 Cone A) 55510 


Fig. 18. — Total number of £®- per collision. 
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Obviously both assumptions give nearly the same result for those quan- 
tities, which anyway hardly change with Qx, namely, the numbers of x, 97, 
A and y from 7° and Y° decay. There we show these numbers in a linear scale. 


Berre O 2 sy te 2 S 10 


The numbers of hyperons and K-mesons depend strongly on Q, and in 
a different way for each assumption, they are shown in a logarithmic scale. 


10° <n? & 
i 
E 
Bi 
1 
2 on Pa 
10° 3 ie ae il Q 
[0e 108882 STE Cp Do 


io 90 = GN = È 
Fig. 20. lotal K-production per collision. 
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We would expect that assumption (i) Qy = 2, with a value of Qx between 
1 1 ae pate : = 
160 fn and Fen should give agreement with experiments. This would agree 


10 + SK 
<n? K Va 


10 
10: 
[o] 
SK pen Sn 
o 2 
10 
‘er A re y ae ee 


Fig. 21. — The ratio K/K = K*/K-. Since the K-mesons are produced in associated 

(77+9>Y+K4+nr) and in pair production (97 +91 >-K+K+nr), whereas the K 

are produced only in the latter reaction, this ratio is > 1. In fact, it rather strongly 

depends on Q, in the «global symmetry » version and is constant in the version 

Qx= Qe. If the experiment shows that K/K > 1, then this definitely rules out the 
version (ii) Qy= x in favour of the other one. 


in a loose sense with global symmetry and a coupling strength for the K-mesons 
which is smaller by about one order of magnitude than that of the 7-mesons. 


APPENDIX II 


The total probability for an event in which an 9797-pair is created, is 
proportional to: 


(A.1) 0* (Ey; Ma + Mn} >p=0)— 


= for dpi dr dro (Lo D'e)0(D pp; + +.) 
i=1 


t=1 
Here 7, is the momentum of the 97 and x, that of the 97 of the pair. 


(Ena = eG ee Va: +1; M=1). Correspondingly, 


(A.2) dr, am, (ap, peo (Lp) 0 (D Pete Ti ae 
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would be proportional to the probability of the same event but with the 97 
having a momentum near to 7m, and the 97 having it near to m,. We can 
therefore single out all those events in which these two particles have less 
than a certain relative momentum, simply by integrating x, and x, not over 
all values but only over a restricted region fulfilling the wanted condition. 

To this end we return to (A.1) and write it in a different way, using the 
identity 


n-2 n—2 
(A.1’) o* =[anaP|ap, … dp,»6(E,— E — > ¢:)6(> pi +P): 
o i=1 i=1 
fam, dit, 6(H — en — Ena) O(P — 1 — Tr). 


By integrating over Æ and P this gives back (A.1). After having thus sepa- 
rated off the pair, we may easily carry out the above idea by introducing the 
coordinates one usually takes in two-body systems. That amounts to putting 


oa 


(A.3) = 
0(p'P’) 


UE ER api pari 


Fi- vp Ppt +4P* 41. 


En — Va 
After integration over P’, we have, on account of the 6-function 


(A.4) fam WO sco 


‘ j I ANI 7 ia ae Fed Lù A n E 
= 2a] n° ap aso(u—|/p ° + Pp'E Ap +1 = y ppi i + 1) 5 


RE PE ter Fe 
Now, p' is half the relative momentum and P the total momentum of the pair. 
We wish to single out cases in which the relative momentum p is small. Then, 


for small P, the 1 is the leading term and we may develop the square roots 
by which the term Pp'E drops out: 


er pi TE pe 
VETRATA / | ñ 12 / P? 
lp pepe + 414|/p eg a. Wes 


à 
bo 
+ 
ae 
Ss 
A 
ely 
== 
a 
t 
= 
+ 
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If, on the other hand P is large, 1+ P?/4 is the leading term and 
AN RA) pi 
Vip + ppe+7+1+|(p?-pP8+7+15 
STD: p” eG CTE pa 
ZI, (1 {| —___~_____ wa V I pa L_, 
| ee CETTE - omy) Ai 


So we may write as a good approximation 


Ò (£ _ |)» i Pp'é = Se 1 — ne pps fe + 1) CS 
lar fire ri) 


Inserting this into (A.4) and (A.4) into (A.1’), we find after integrating over H 


Di 
(ATO dpr ta nap’ [a id pied Bs 
0 


n-2 


re PN ee n—2 
-6(Ey— V4 + 4p” + P? Le)ò(X pit P) - 


This, since we integrate only up to po, is proportional to the number of events 

van) . . . lo ! . 
where the 97 77 pair comes out with a relative momentum pa < 2po. Finally 
we put p=2p' and have 


Do 
(A.1”) 60° = | pean far, dp dhe 
<3 
EVE +p) +P? — > e)6(S ps tP)- 
i=1 t=1 


Now the second integral is nothing else than 0* (Eo, Ma +++ Mn—2, MS, > p= vy 
that is the o* for the n —2 first particles plus one more particle (namely 
the 97 Z-pair) with mass 


m*= VA +pre 2 + 
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which has the obvious meaning of a particle consisting of two nucleons plus 
the kinetic energy of their relative motion (remember m= 1). 
Therefore 


Po 


(A.5) 00° = [pare Ma Ma (DI) 


0 


To this we add the further factor (Gannibi/Otota)(P), to account for the experi- 
mental fact that not all low-energy p-p collisions result in annihilation. The 
o* (Hy, M ... Mn, m*) is then computed with the usual Monte Carlo method 
for four different m*(p) and the integral is calculated by Simpson’s formula. 
We call Ao* the integral (A.5) including the o,/o,-factor. 

We choose p, = 1, since this is the order of magnitude which corresponds 
to the 9 and 97 staying together nearer than a pion Compton wavelength 
for longer than one elementary time. 

The value of p, and the factor (Ganninn/Otora)(p) are of course somewhat doubtful. 
However, the contributions to Ao* coming from p < + and p > } are not 
important; for the small p the p? decreases too strongly and for the larger p 
the m*(p) takes away too much energy, hence 0* (Hig, Mas Mn+, m*(p)) , 
becomes small. Consequently, if the rate of the « recombination » of 97-97 pairs 
immediately after the collision is a smoothly varying function, one can account 
for the portion of 97-97 pairs which scatter without annihilating, by intro- 
ducing a suitable factor (<1) multiplying do*. In the present case we took for 
this factor 0,/0,, which is nearly 4 over the whole range of p we have to con- 
sider here. Between p — 0 and the lower limit of the experiment, we assumed 
a 1/0 behaviour of o,, but the contribution of this region to Ao* is quite unim- 
portant. Because of the vanishing of o* for large p, it is also unimportant (in 
the present case where the 97 can have only a small kinetic energy) where we 
put the limit p,, as long as Po> +, which seems a reasonable value. 

We still have to account for the fact that there are three nucleons with 
which the antinucleon may annihilate, whereas our calculation has tacitly 
assumed that only the last one could. Let us visualize Ao* and o* as volumes. 
Then in 9* there are one antinucleon, three nucleons and may be other par- 
ticles. If there were only one nucleon and one antinucleon, then 97-(Ao*/o*) 
would be the number of 97 that would annihilate. In other words: if the 
nucleon is in Ag*, annihilation takes place and Ao*/o* is the probability for 
finding the nucleon in Ao*. 

This is a situation well known in elementary statistical gas theory: we 
have particles (nucleons) in a box of volume o*. Inside there is a smaller 
box of volume Ao* and we just have calculated the probability to find a 
definite nucleon in Ag* no matter where the others are. This, however, is 
not what we need, since annihilation will take place as soon as at least one 
of the n nucleons is in ANG, In complete analog 


1 gy to the elementary statistics 
we find that the probability for this is equal to 


1—[probability that no nucleon is in Ao*] =1 (1 a) 
~ = * . 


0 


266 


PARTICLE PRODUCTION IN 6.2 GeV p-p COLLISIONS ETC. 267 


Hence the total number of annihilated antinucleons will be 


n—afi— (22), 


if n nucleons are present. 

In our case n—3. If there are other particles (pions, etc.) present, then 
these three nucleons may be considered as independent and the above con- 
sideration is correct. That is true for slightly higher energies than ours. Here 
we have no other particles and therefore we will have a correlation between 
the nucleons, because the sum of their momenta must vanish. So only two 
nucleons are really free, though not completely independent. It will be prob- 
ably not a bad approximation to account for the correlations of the nucleons 
by considering two of them as completely free in the volume o* and the third 
to be completely fixed by the position of the first two. That amounts to 
put »=2 instead of 3 as one would do if the mesons could balance the 
momentum. This yields the formula quoted in the main text. 

There remains the correction of the 97-spectrum. We have done this in 
a rather unsatisfactory way and some steps taken are not proved in detail. 
We think, however, that it allows at least a good qualitative guess. 

Suppose again, for a moment, that only one of the three nucleons can 
annihilate with the antinucleon. We wish to calculate the spectrum of those 97, 
which then would annihilate, i.e. that A97 spectrum which had to be sub- 
tracted from the 77-spectrum in order to yield the observable 977,,-Spectrum. 

The 9797-pair with relative momentum p, has been shown to be equi- 
valent to a single particle with mass m*=V4+p?. In the calculation of 
o*(H, ... m*) the spectrum of this fictitious particle m* is also computed and 
using the subprogramme for the decay of this m* particle into two particles 
ofm=1 one obtains a spectrum which is just the spectrum of the 77 and 
of the 97 under the present restrictions on their relative momentum, 1.6. in 
particular the spectrum of the annihilating 97. 

Therefore one had for each p, to calculate the corresponding qgt-spectrum 
and add them up by Simpson’s rule, thus obtaining a total A97-spectrum. 
This has been done. 

Now comes the difficulty: this would be correct for the case that one 
single nucleon is available for annihilation with the 9. In fact there are 
three and our above argument concerning the overall effect of annihilation 
should now be applied to the integrand for every value of Pri. But here the 
correlation must be taken into account and the question becomes much more 
complicated as in the overall effect. We do therefore an approximation: the 
AQ7-spectrum is calculated along the lines given above as if only one nucleon 
were present. Then it is renormalized such that its norm equals the overall 
AQ = F1— (1— Ao*/0*)"]. Then the original -spectrum is normalized to 97 
and the AQt-spectrum is subtracted from it. The resulting spectrum is con- 
sidered as the «observable 97-spectrum ». The philosophy behind this is that 
the most important effect is not so much, that at the same time more that 
one nucleon can be near 97 (inside Ao*), but rather that the probability to 
find just one has become larger because there are three nucleons each of which 
is able to enter in Ao*. Therefore the spectrum of the annihilating 908 is 
given mainly by those states in which only one nucleon comes near to 97 
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(that is, what we have in fact calculated first), but this happens much more 
often than our calculation gives; hence we correct it by renormalizing the 
spectrum of the annihilating 97 to have a norm equal to the total A7. 

If, on the other hand, the probability of two, or more, nucleons being 
in Ao* becomes large, then the effect of there being more than one nucleon 
will turn out to be different for different relative momenta and such an over-all 
renormalization is certainly wrong. We may be just on the border line. 


RIASSUNTO (*) 


La teoria di Fermi per la produzione di particelle viene usata in forma rigorosa 
e raffinata allo scopo di calcolare gli spettri e i numeri medi di particelle per tutte 
le particelle prodotte (cioè 97, x, 97 e particelle strane) in urti p-p a 6.2 GeV (sistema 
cinetico del laboratorio). 


(*) Traduzione a cura della Redazione. 
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Large Amplitude Waves in a Collision-Free Plasma. 
I. - Single Pulses with Isotropic Pressure (*). 


A. BANOS jr. and R. VERNON 


Department of Physics, University of California - Los Angeles, Cal. 


(ricevuto il 13 Ottobre 1959) 


Summary. — We consider an infinite expanse of low density, fully ionized 
plasma, magnetically inmobilized in a constant and uniform field of mag- 
netic induction. In the absence of collisions we assume that the ion and 
electron motions take place in planes perpendicular to the magnetic field. 
We then transform to a «shock » frame of reference moving at constant 
speed at right angles to the magnetic field, and then we search for the class 
of non-trivial, time-independent, one-dimensional, self-consistent solutions 
of the Maxwellian set and of the equations of motion of the charged par- 
ticles; that is, we study the propagation of transverse pulses and waves 
whose characteristic length is much smaller than the collision mean free 
path. Upon making appropriate simplifications, we are led to a system | 
of equations in which collisions are neglected, charge neutrality is as- | 
sumed, and both ions and electrons behave like two-dimensional fluids 
obeying isentropic relations with y=2 and exhibiting rigorously isotropic 
pressure tensors. With these simplifications and assumptions, the system 
of equations can be solved completely by a simple numerical quadrature. 
Both symmetric pulses and periodic waves are obtained, with characteristic 
lengths of the order of the plasma «skin depth », or mean gyromagnetie 
radius for particles traveling with the Alfvén speed. In this paper, Part il. 
we describe the solitary pulses, which are the only solutions satisfying 
the conditions of the undisturbed plasma ahead of the wave train. Two 
basie parameters are required for a complete specification of the problem :, 
the Alfvén Mach number a, which gives a measure of the speed of pro- 
pagation, and the ratio f of the total plasma pressure to the magnetic 
pressure, which specifies the initial state of the plasma. It is shown that, 
for given f, stationary solutions exist only for a limited range of speeds «. 


(‘) Supported in part by the U.S. Atomic Energy Commission through the Uni- 
versity of California Lawrence Radiation Laboratory, Berkeley and Livermore. 
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1. — Introduction. 


The problem of the one-dimensional shock in a collision free plasma fur- 
nishes the principal motivation for the present study. As things stand to- 
day (1), it has not been established, theoretically or experimentally, whether 
or not a shock-like phenomenom exists in the absence of collisions or other 
dissipative mechanisms. But it is known (?), that in the absence of collisions 
there exist non-trivial steady-state solutions which may serve as a guide to 
the whole question of one-dimensional shocks. In this paper we consider à 
two-dimensional plasma magnetically immobilized by a constant and uniform 
field of magnetic induction. We then transform to a «shock » frame of refe- 
rence moving at constant speed at right angles to the magnetic field and inquire 
into the existence of stationary, one-dimensional, self-consistent solutions of 
the Maxwellian set and of the equations of motion of the charged particles. 
Among the steady-state solutions (in the shock frame of reference), we find 
solitary pulses and periodic solutions. 

In this paper, Part I, we describe the structure of large amplitude solitary 
pulses and their domain of existence for a plasma supposed to be endowed 
with rigorously isotropic pressure. In this respect, our results constitute an 
extension of earlier work by ADLAM and ALLEN (*) and by Davis, LÜST and 
SCHLÜTER (4), who consider only the zero temperature case. 

In the sequel to this paper, Part IT, we propose to discuss the structure of 
large amplitude periodic solutions (waves) and their relation to the solitary 
pulses, again for a plasma with isotropic pressure. There is a close connection 
between our work, in both Parts I and IT, with that of MoRAWETZ (?) in the 
absence of the friction term, which will become apparent in the course of this 
discussion. We are at present extending this work to include the non-isotropic 
components of the pressure tensor for the special case of equal masses. These 
results will be reported separately in Part III. 


2. — Fundamental equations. 


The model that we consider here envisions an infinite expanse of neutral 
plasma consisting of singly-ionized ions and electrons in thermal equilibrium, 
embedded in a constant and uniform field of magnetic induction parallel to 

(*) C. S. GARDNER, H. GOERTZEL, H. Gran, C. 8. MorAWwETZ, M. H. Rose and 


H. RUBIN: Hydromagnetic Shock Waves in High Temperature Plasmas, NYO-2538, 
Inst. Math. Sci., N.Y.U., January 30, 1959. 


7 
(*) J. H. ADLAM and J. E. ALLEN: Phil. Mag., 3, 448 (1958). 
(4) L. Davis, R. List and A. SCHLÜTER: Zeits. f. Naturf., 13a, 916 (1958). 


N 


d» 


LARGE AMPLITUDE WAVES IN A COLLISION-FREE PLASMA - I 271 


the z-axis of a laboratory Cartesian frame of reference. Because we neglect 
collision phenomena altogether, it is convenient to assume that the charged 
particles execute circular orbits in planes perpendicular to the magnetic field, 
thus giving rise to what we might term 
a two-dimensional, collision-free plasma. y 
Next, we make a Lorentz transforma- 
tion to a shock frame of reference mov- 
ing at constant speed w, in the negative 
æ direction, with respect to the labora- 
tory frame, and then we look for the 
class of non-trivial time-independent 
(stationary), one-dimensional solutions 
of the Maxwellian set and of the equa- 
tions governing the motion of the 
charged particles. To simplify the no- 
tation we adopt here the convention 
of designating ion quantities with capital 
letters and electron quantities with the corresponding lower case letters. 
Accordingly, the electromagnetic field vectors € and Z (see Fig. 1), and 
the ion and electron fluid velocities V and v exhibit respectively, in the shock 
frame of reference, the components 


E =[(E(x), UoBos 0] > (0, % Bo, 9) 
B—=1[0, 0, Z(x)] — (0, 0, Bo) 
V =[V,(x), V,(æ), 0] > (wo, 0, 9) 
[v,(@), v,(æ), 0] > (wo, 9, 9) 


I 


Dv 


which reduce to the indicated limiting forms as 7 > — co (boundary condi- 
tions), and where % is the speed of the shock frame and B, is the magnitude 
of the original field of magnetic induction. 

In the adopted notation, the fundamental equations for the present one - 
dimensional stationary case can be listed as follows (°): 

Maxwell’s equations, in rationalized mks units, reduce to 


(1) e(dé/dæx) = e(N—n), 
and 
(2) dZ)dæ = pe(nv, — NV,), 


(5) For the derivation of the equations in general form see, for example, H. GRAD: 
Notes on Magneto-Hydrodynamics. I. General Fluid Equations, NYO-6486, Inst. Math. 


Sci., N.Y.U., August 1, 1956. 
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where e is the electronic charge, N and n are the number densities of ions and 


electrons, respectively, and pe din 
The continuity equations (conservation of mass 
which can be integrated at once, yield respectively 


) for ions and electrons, 


(3) NV, = constant = NU» ; 
and 
(4) nv, = constant = NU » 


where », is the number density of ions and of electrons as æ > — oo. 
The equations of motion (conservation of momentum) for ions and electrons, 
become respectively, for the # component of momentum, 


(5) NMV,(dV,jdx) =— dP/du + Ne(& + V,2), 
(6) nmv,(dv,/dx) =— dp/dx — ne(ë + VB) » 


and for the y component of momentum, 


(7) NMYV,(dV,/da) =  Ne(uoBo — Va) ; 
(8) nmv, (do,jda) = — ne(uoBo — iP) . 


Furthermore, assuming that the pressure tensors for ions and electrons 
remain rigorously isotropic, which is equivalent to saying that the respective 
distribution functions behave, in the course of the phenomenon, as local Max- 
wellian distributions with number densities and temperatures which vary only 
as functions of #, we obtain the second moment equations (°) 


(9) a(V,P)/da =: — P(dV./dx), 
and 
(10) d(v,p)/da =— p(dv,/dz) , 


which we can integrate at once to yield 


72 2 
(11) PY, = constant — piu, 
and 
9 Anse 
(12) pr, = constant = pus, 


(9) L. SPITZER jr.: Physics of Fully Ionized Gases (New York, 1956), Appendix. 


a 
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where p, is the pressure of ions and of electrons as x —>— oo. We observe, 
therefore, that both fluids satisfy the isentropic relation for ideal gases with 


y = 2. 
Tn addition to the immediate integrals (3), (4,) (11), and (12), which result 


from the stationary one-dimensional character of the problem, we have also, 
quite generally, the classical integrals of the motion which arise from the con- 
servation laws of momentum and energy for the whole plasma (5). In the 
present instance we have, for the x component of momentum 


m3) (NV, + PIE (nmv? + p) — 386? + (B— B?)/2u = 2Po + NolM + Mu 
and for the y component of momentum 
B14) NMV,V, + nmozv, — EUjpD 6 = 0 x 


Finally, the conservation of energy for the whole plasma yields, in the pre- 
sent case, | 
(15) (ANM(V2+ Vi) + 2P]Vat [Anm(0; + 05) + 2p)]10 + 

Lu BB — Bo) = 4Psvo + MM + Mu, - 


3. — Charge neutrality. 


The set of equations (1) through (10) constitutes a closed system in the ten 
dependent variables PC aN, 13, Vay. Vai Pond Pap. Phe solution of 
this system as listed is intractable by analytic methods and extremely cumber- 
some (if not impossible) by numerical methods because of its complexity and 
because of the recognized smallness of the charge separation, 


(16) IN—ni<n; 

which appears in Poisson’s equation (1). As shown in Section 8, this condition 
is satisfied to a high degree of accuracy in most cases of interest. Therefore, 
to achieve a practical solution of the original system we assume, at the outset, 
that the plasma remains electrically neutral throughout the phenomenon, 
(17) N=nH, 

as a consequence of which (3) and (4) yield 

(18) Vans 


ca 


N 
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while the integrals (11) and (12) then give 


Thus, to this approximation, ions and electrons go hand in hand with identical 
horizontal velocities and identical pressures and temperatures. 

Further reduction of the system is achieved by noting that eqs. (7) and 
(8), making use of (17) and (18), yield at once upon addition 


(20) MV,+ mv, =0, 


which, in conjunction with (14), makes it evident that the assumption of charge 
neutrality (17) is tantamount to letting e + 0 (or ¢ + co) in (1), (13) and (14). 
Thus, with the aid of (18), (19), and (20) we can express the ion quantities 
in terms of corresponding electron quantities, and with the aid of (3) and (4) 
we can eliminate the number densities 


(21) Ni Sato Ons 


In this manner we reduce the original system of ten equations to the follow- 
ing five: 


(22) Mo,(4Z/dæ) = uens(m + M)uyr, , 

(23) Mv,(d/dx) [v, + 2p/no(m + Mu] =— ¢v,8 , 
(24) mv,(dv,/dx) = e(v,8 — Bo) , 

(25) v,(dp/dx) = — 2p(dv,/dx), 

(26) 2Meé = (M — m)[Mv,(dv,/dx) — ev,B] . 


Thus, (22) through (25), constitute a closed set in the four dependent variables 
B, Ver Vyy and p; once this system is solved, eq. (26) yields the longitudinal 
electric field. We note that we could readily make use of (25) to eliminate p 
from (23), but we choose not to do so at this juncture for reasons that will 
become apparent in Section 5. 

In addition to the basic equations listed above we have the integrals cor- 


responding to the conservation laws of momentum and energy. Thus, making 
the same substitutions as before, the 


x component of momentum (13) can now 
be written as 


(27) Moto(M + M)(v, — uo) + 2(p — po) + (B* — B)/2u = 0, 


7 


Dae oe 


bo 
“I 
or 
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= 


while the conservation of energy (15) reduces to 
(28) inouo(M+ M) [v3 + (m/Mjv, — uz] + 4(pv, — Polo) + UBF — Bo) = On 


Finally, we add the integral of (25), which we interpret as the conservation 
of entropy and which reads, in accordance with (11) and (12), as 


(29) pr? = constant = a*p, v4 , 


where a >1 is a convenient parameter that we introduce here for ulterior 
reasons (Section 5) and which must assume the value a=1 for solutions 
which satisfy the boundary conditions as a > — co. 


4. — Reduction to dimensionless form. 


To facilitate the solution of the preceding equations it is useful to express 
the field vectors £ and 4 in terms of the dimension-free variables 


(30) E = €/uBo and Bb 


and to express the horizontal and vertical components v, and v, of the electron 
fluid velocity in terms of the variables 


(31) U = Uy] Uo and D = (v,/Uo)(m/ M) : 


In addition, we replace p, wherever it appears in the preceding equations, by 
the product pp,, where now p is a dimension-free variable and p, is the pres- 
sure of ions and electrons at minus infinity. Finally, instead of the element 
of length dz, we introduce the dimension-free element 


(32) dé = dar, , rm =(mM)} uofeBo ; 


that is, we express lengths in terms of the geometric mean gyromagnetic radius 
r.. for ions and electrons at the mean speed and magnetic field prevailing at 


m 
minus infinity. 

To further facilitate the reduction to dimensionless form we introduce 
Alfvén’s Mach number 


A 


(33) CAPE v = Bilun(m+ M), 
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where ©, is Alfvén’s phase velocity for the conditions of the plasma at minus 


infinity, and the parameter 
(34) p= 2po](Bi/2u) ’ 


which gives the ratio of the total ion and electron pressure to the magnetic 
pressure. Thus, the entire problem is expressed in terms of the two para- 
meters x and 8, which measure respectively the shock strength and the initial 
temperature of the plasma. For numerical purposes we combine the two para- 
meters in the ratio 


(35) v8 = Blo? = 4p,/[n (m+ Mus], 


where vy? is a convenient parameter which expresses the square of the ratio of 
the sound speed to the shock speed. 

Making use of the dimension-free variables and parameters listed above, 
the system of four non-linear, first-order differential equations (22) through (25) 
now becomes 


(36) u(dB/dé)= ov, 

(37) u(d/dé)(u+3v°p)=—vB, 

(38) u(dv/dé) = uB—1, 
(39) d(up)/dé = — p(du/dé) , 


to which we must add the integrals (27), (28), and (29), corresponding res- 
pectively to the conservation laws of momentum, energy, and entropy, that 
now read 


(40) u—1-+ 4°(p — 1) + (B?—1)/2e2 =0, 
(21) Uro TE pt (pu — LE OCB 2120 
(42) pu? = constant = «a. 


De, the longitudinal electric field, as deduced from (26), becomes in dimen- 
sionless form 


(43) DEC m)(mM)*[u(du/dE) — vB] = 


Meme 3%%plu)](1 — rv p/u), 


wherein the second form is obtained by eliminating the term u(du/dé) with 
the aid of eqs. (37) and (39). We observe that, in the special case of equal 


© 
N 
N 
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masses (M = m), the condition of charge neut ‘ality (17) becomes rigorous and 
the longitudinal electric field vanishes identically. 


5. — Singular points of the differential equations. 


The singularities of the system of differential equations (36) to (39) occur 
where (?) 


(44) WB and DEA? 


According to (40), (41) and (42), there are three such singular points, one of 
which occurs for negative values of # and is therefore immediately discarded. 
The remaining two singularities are respectively a saddle point and a focal 
point. 

The saddle point occurs when »=0 and 


I 


(45) (ee eels il, ea hen cones in, 


which are in fact the boundary conditions prescribed as £ > — co. The nature 
of the singularity is readily ascertained by expanding the solution, in its vici- 
nity, in the form 


B=1-+B;exp [Ae] 


u =1+ uu; exp [AE] 


(46) 
p =1+ p, exp [A6] 


v = v, exp [46], 


where A is the characteristic exponent and B,, 1, Pi; and v, are the coeffi- 
cients of the first variation. As is well known, upon substituting the expan- 
sions (46) into the primitive set of differential equations (36) through (39), 
there results a homogeneous system of algebraic equations in the coefficients 
of the first variation whose determinant must vanish, giving rise, in this 
instance, to a quadratic secular equation in A? with roots 


(47) A? = 0) and A Sale 
The pair of zero roots corresponds merely to the trivial solution in which all 
derivatives vanish and the variables retain their values fixed by (45). The 


second pair leads to real roots, which is what we are after in accordance with 


= 
~ 
A 
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(46), provided the Alfvén Mach number « satisfies, for a given » (or, alter- 
natively, for a given B), the inequality 


(48) = Cee. ae (1 =) = 5, 
which fixes the lower bound of « This analysis establishes the saddle-point 
character of the singularity provided (48) is satisfied. 

The focal point occurs where v=0 and 


3012 
(49) HDI, G@=a;, and hence p, = a;/w,, 


where, in accordance with (40), (41), and (42), we must have 


(50) u, =1/B, =4[1 + 2(1 + Blot] <1, 
and 
(51) t(ra,)° = [1 + $(1 + B)/a?]u? — u? — 1/20 ; 


that is, the co-ordinates of the focal point correspond to the same momentum 
and energy as the saddle point, but exhibit a higher entropy (a, >1). In 
the words of a hydromagnetic analogue, this is the state that would be ob- 
tained behind the shock from the familiar Rankine-Hugoniot relations. Pro- 
ceeding as before, we now expand the solution in the vicinity of the focal 
point by using expansions similar to (46), obtaining a corresponding quadratic 
in A? with roots 


(52) ATO and A? = B;[a® — BiJ(u,; — v°p,)]. 


The pair of zero roots again corresponds to a trivial solution, and the second 
pair, it can be shown, leads for all permissible values of the parameters to 
pure imaginary roots. That is, in the infinitesimal neighborhood of the focal 
point (49), our system of equations admits periodic solutions, which thus 
establ'shes the character of the singularity. 
tions in the sequel to this paper, Part II. 

It is evident that the complete numerical solution of the problem can be 
obtained from the integrals of the motion (40), (41), and (42). Thus, making 
use of (42) with a=1, and assigning suitable values to u, as explained in 
Section 7, we can compute B from (40); then, making use of (41) we can com- 
pute the corresponding values of v. In this manner, with a=1, we obtain 
the solitary pulse illustrated by plotting v vs. (B — 1) 


We return to these periodic solu- 


as shown in the pseudo- 


7 


OP 
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phase diagram of Fig. 2. Proceeding similarly for 1<@<4;, one obtains 
the closed curves or periodic solutions (waves) shown in the diagram. Finally, 
putting a= a; merely yields the focal 

point, which may be regarded as the =f 


: : ets a=! (Solitary pulse) 
condensation point of periodic solu- 


tions as a > ay. In this paper, Part I, 
we discuss in detail the structure of 
the large amplitude solitary pulses for 
all permissible values of the parameters, 
and we reserve for the sequel to this 
paper, Part II, a thoroughgoing dis- 
cussion of the structure of the periodic 


1<a <q; (periodic 
solutions) 


a=a,(focal point ) 


Fig. 2. 


solutions or stationary waves, which is 
a somewhat more complicated matter. 
For example, we have discovered that the periodic solutions do not neces- 
sarily exist for all values of the Alfvén Mach number « for which the solitary 
pulses do exist. 


6. — Upper bound on Alfvén’s Mach number. 


As stated in eq. (48), the Alfvén Mach number x must exceed the lower 
bound «_ in order to obtain the solitary pulse solutions described above. We 
will now show that for a given » (or, alternatively, for a given 6), the value 
of « must not exceed the upper limit 


(53) a = 4/(14+3v— 498), 


beyond which the solution of the system breaks down because of multiple- 
valued solutions in €; that is, beyond which solutions of the postulated type 
do not exist. 

For the purpose, making use of (42) with a=1, we rewrite (40) and (41) 
in the convenient forms 


(54) plu; v) = (B?—1)/20? = 1 — 4 — 4y3(u-? — 1), 


(55) C(u, 03 v) = 2(B —1)/o = 1 — w? — 23(u-! — 1) — 0. 


We note that (54) expresses 7 (or B) qua function of uw, for a given vy (0<v<1), 
and that (or B) attains à maximum when 


(56) TI 
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Thus, we argue that the significant range of values of u, for a given », lie in 


the segment 


(57) POULE 


and that, to obtain a solitary pulse solution (see Fig. 2), it is necessary for (55), 
regarded as an equation for v?, to have at least one root (other than u= 1, 
B=1)in the range (57). It then becomes clear that the maximum value of x, 
a=,, must abide when 1? =0 and u=v in (54) and (55) simultaneously ; 
that is, when (55) has a double root for v? at w=». 
To determine this value of « we first compute from (54) and (55), putting 


u“=y and v=0, 


(58) n(v;») = (BE 


(59) C(v, 0; v) = 2(B,, —1)/oe®= — 1 — 3y? + 23 = (1 — v)?(1 + 27), 


from which, by division, we 
readily compute the upper 
limit, for a given », to the 
maximum value of B, 


(60) By, = 3/(1+ 2). 


Inserting this value of B,, 
in (59) we now solve for at 
obtaining 


(61) Ala =(1 — v)(1+27)?= 


— 5 of 4) 


which is equivalent to (53). 
Do È | Thus, we have shown that 
Dio solitary pulse solutions exist, for a given » (or, alternatively, for a given 
6) in the range | 


20 
(62) CS Ch 


VALE SE is acy (48) and x, by (53). This allowed range is illustrated in 
Fig. 3 which exhibits « vs. B, from B=0 to B=}. 
=. 


RS ae In this diagram the 
ocus of v= constant is given by a set of half-lines 


emanating from the origin. 


(=) 
+00 
N 
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7. — Numerical results. 


As already explained (Section 5), the complete solution of the problem 
can be obtained from the conservation integrals (40) and (41), or alternatively, 
from (54) and (55). That is, the procedure of assigning values to uw and com- 
puting B and v, for a =1, leads numerically to the pseudo-phase diagrams 
of Fig. 2. To exhibit the dependence of the several variables qua function 
of È (or x), we need to perform a quadrature. Because « is a convenient 
variable of (numerical) integration we return to (37), and making use of (39) 


to eliminate p with a=1, we obtain 
(63) (i —— 9" lar?) (dude) = vB, 


which we can at once integrate in the form 


(64) CÈ IC — yy?) /v(u) B(u)] du, 


where B(u) is determined from (54) and v(u) from (55). Here the lower limit 
of integration w_, which must lie in the range (57), corresponds to the maximum 
value of B (Fig. 2) which, according to (3€), occurs when v= 0. This value 
of ‘u_, for a given y and a given « within the range (62), must be determined 
numerically as a root of (55), with v=9. However, for x= a+, we have 
simply as we have seen u_=v. At the lower limit, then, we have v(u_) = 90, 
which however, is an integrable singularity behaving as (u — u_)*, except when 


— y, in which case the integrand remains finite. At the upper limit, 
il 
9 


CANE 
on the other hand, as u->1, we have v-> 0, but this time as (1 — y) 
giving rise to a logarithmic infinity in accordance with the exponential be- 
havior (46). Since, for a given u, v assumes either positive or negative values 
(see Fig. 2), it is clear that, depending on the sign of v in (64) we obtain either 
branch of a symmetric pulse starting at — co with the prescribed boundary 
values (45) and returning to these very same values as È > + 00. 

Before describing our numerical results for specific values of « and |, it 
is well to observe that our chosen unit of length defined by (32) depends on 
the magnetic field and hence on the shock strength «. Therefore, to compare 
properly curves corresponding to different values of a, it is desirable to employ 
a unit of length first proposed by ADLAM and ALLEN (3), and independently 
by ROSENBLUTH (7), which is a property of the plasma, but independent of 


(7) M. ROSENBLUTH: Dynamics of a Pinched Gas, in Magnetohydrodynamics (Palo 
5 


Alto, 1957; edited by R. K. M. LANDSHOFF), p. 97. 
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the magnetic field. To this end, we define the normalized dimensionless ele- 
ment of length dé, defined as follows: 


(65) dv = = Von dé = a "y dé n9 nm = (mM)*0,/e Bo = Td 


whence we have dé,= (r/r) dé = «dé, with rn as previously defined by (32). 

We present results for three n corresponding to 6 = 0, B =1/49, 
and p=1/4. For each temperature we display the profiles of the pertinent 
dimensionless variables B, uw, 
v, and [(mM)}/(M —m)]}E as 
functions of the normalized 
distance £,, for the upper 
bound «, of Alfvén’s Mach 
number, as computed from (61), 
and for the two higher tempe- 
ratures we also present the 
corresponding profiles for an 
intermediate value of « in the 
range (62). 

Fig. 4 exhibits the familiar 
results for the zero tempera- 
ture case, with «=a,=2, 
which ADLAM and ALLEN () 
had. already published and 
which a number of other in- 
vestigators had at various times 
obtained. We include these 
results here for comparison 
purposes and because our re- 
sults include the zero temper- 
ature limit as a special case 
(v=0). We note that, in Fig. 4, 
B attains à maximum B,, — 3, and % a minimum w_ = 0, but that in both 
cases the curves have à discontinuous derivative with cu usps which exhibit 
vertical tangents. To show quite generally that, for «—a,, the curve of « 
vs. € exhibits a cusp at &= — 0, we simply compute the derivative du/dé 


from (63) after putting u=v1+o), d<1, and expanding in powers of 6. 
In this manner we obtain, with «=«,, which implies w=», 


30F 


25 


Fig. 4. 


(66) (du/dé),_, ut vm} H( ++ 2v) — = ©) 


v—+>0 


o, 
the double sign arising from the sign of » (see Fig. 4) before taking the limit 
OI) Thus, in the zero temperature case (v= 0), the curve of « vs. CEN. 
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at u = u_=0 a cusp with vertical tangents. The fact that u — 0 for £ = 0, when 
~»—0 means in accordance with (36) that the profile of B also exhibits a cusp 
in this case. These are the familiar results which we associate with the upper 
bound «,=2 at zero temperature, and which are interpreted by saying that 
an increase in x beyond x, would produce looping orbits whith the consequent 
breakdown of the present solutions. 

Figs. 5 and 6 display, for x=x;, the corresponding profiles of B, u, v, 
and Æ vs. €, for the two higher temperatures chosen. This time, however, 
the curves of « vs. É exhibit cusps at u=u_=v>0 with finite tangents in 
accordance with (66). As a consequence because #_ > 0, the corresponding pro- 
files of B vs. £ are continuous with a continuous first derivative everywhere 
and in particular at £ = 0. Increasing the value of « beyond «, again results in 
a breakdown of our mathematical solutions, but this time it is not sufficient 
to merely invoke as a physical reason the appearance of looping orbits, since 
the plasma has an initial temperature and a fraction of the orbits would be 
looping in any case. The breakdown of the solution appears to come from 
the fact that, for x >«,, the horizontal velocity w would become multiple- 
valued in £, which of course has no physical meaning, indicating that in this 
case we would probably set up some turbulent regime leading to solutions 
not included in our one-dimensional stationary flows. 

Another important feature which is evident in Figs. 5 and 6 is the 
fact that the longitudinal electric intensity E is discontinuous at È= 0, in 
contrast to the continuous curve of E vs. É in the zero temperature case 
(Fig. 4). To show how this discontinuity comes about we return to (43), which 
exhibits the longitudinal electric field in terms of the (discontinuous) deri- 
vative du/dé of (66). Thus, putting, u = y(1+6) wtih 6<1 we compute from 
(33) the values of E at È =0 after letting Ô >0. In this manner we obtain, 
for the discontinuity in È, 


(67) AE = E(0+) — E(0—) = (M —m)(mM)*[1— v)v]}}/(1+ 2), 


which vanishes for v=0 as illustrated in Fig. 4 These results indicate, as 
discussed further in the next Section, that for x =, the assumption of 
charge neutrality breaks down completely in the vicinity of the cusp. 

Figs. 7 and 8 correspond to the same values of B as Figs. 5 and 6, respec- 
tively, but display the profiles belonging to intermediate values of x, a<a< 4%. 
The essential feature to observe is that all the curves have continuous deri- 
vatives and that the cusps have disappeared. Thus, for any « within the spe- 
cified range and for a=1 in eq. (29), the solutions are solitary pulses with 
the characteristic symmetric behavior depicted in Figs. 7 and 8. 

Finally, to exhibit the changes brought about by an increase in tempe- 
rature while keeping «=«,, which in turn depends on the temperature in 
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accordance with (53), we have plotted in Fig. 9 the behavior of B and wu vs: &,, 
for the three chosen values of B: B— 0, 6 =1/49 and 6 = 1/4. Similarly, Fig. 10 
shows the behavior of + and Æ for the same values of 6. It is readily seen, 
in agreement with (67), that initially at least, the discontinuity of the electric 


Fig. 9. Fig. 10. 


intensity at the cusp increases with the plasma temperature. This peculiar 
behavior suggests that we should re-examine our assumption of charge 
neutrality. 


&. — Charge separation. 


Returning to Poisson’s equation (1) and expressing it in dimensionless form 
we obtain for the fractional charge separation, making use of (43), 


(68) (N—n)/m = (m+ M)(mM) (0, /¢)?(AH/4€) = 


= + (M m?)/mM |(v,/¢)?(d/dé) [u(du/dé) — vB], 


and it is clear that, as long as du/dé behaves properly, that is, as long as the 
derivative du/dé is continuous, the fractional charge separation is essentially 
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governed by the numerical factor 


(69) [CM? — m?)/mM] (v,/¢)? = (1019/n,)Bî in mks units, 

which, except for very low densities and extremely high magnetic fields, is 
considerably less then unity. Therefore, our fundamental inequality (16) is 
generally well satisfied provided we keep away from the singular behavior of 
E in the vicinity of the cusp for a = «.. 

On the other hand, our assumption of charge neutrality, as already pointed 
out, completely breaks down when ~«=«a,. In this case £ is actually dis- 
continuous at £ — 0, see eq. (67), with the result that the fractional charge 
separation becomes infinite in the vicinity of the cusp. In fact, putting x = a. 
and u=v(1+0), d<1, and going to the limit as 6 — 0, we obtain from (68), 
for the fractional charge separation at the cusp, the limit 


(70) (N — no) a> (MM — m°)fmM\(vje) {(1 — v)/(14-20) + [(1 — v)v}}}/36v?, 


which becomes infinite as 1/0»?; that is, the more so as the initial plasma 
temperature is lowered (v — 0). 


9. — Discussion. 


We have obtained solutions with the desired property that the character- 
istic length 7, can be much smaller than the mean free path. Furthermore, 
owing to the fact that the lower bound of wu is », which decreases toward zero 
as B, increases, there is no upper limit on the possible amplitudes of the den- 
sity and pressure. The magnetic field, on the other hand, never exceeds three 
times its initial value B,. The pulses connect two identical equilibrium states 
and therefore do not correspond to hydromagnetic shock waves. This is a 
consequence of the conservation of entropy (29) which, in turn, follows from 
the neglect of all possible dissipative mechanisms. In Part TI we shall examine 
the consequences of relaxing this restriction, and obtain the conditions under 
which these waves can develop into stationary shocks connecting two unlike 
States. 

Our present results constitute the simplest possible extension of the zero 
temperature work (31) to include thermal motions of the particles. We ob- 
serve that, in general, our equations and results reduce to the zero temperature 
case if we simply set 6 or » equal to zero. The outstanding exception is the 
breakdown of the assumption of charge neutrality in the vicinity of «,. The 
factor (u? — 4v)/(u? — v°), appearing in (43) making use of (29), which governs 
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the behavior of (N — n) near «= », is simply unity if we set y — 0 from the 
start. We also observe that the approach to the limit 6 = 0 is rather steep 
so that the behavior even for fairly small B differs considerably from the 
idealized case B = 0. 

Examination of the effect of charge separation appears to be a fruitful 
subject for further study. Our formulation may be regarded as the zero-order 
approximation from an expansion in powers of (v,/¢)?(M? — m?*)/Mm. Higher- 
order terms can be obtained by inserting the charge separation (68) in our 
original equations. The breakdown of our charge neutrality assumption is 
not responsible, however, for the failure to find solutions for «>a,. This 
failure is a result of the appearance of physically meaningless double-valued 
solutions which do not extend to x —> co. Before returning to the question 
of charge separation, therefore, we shall discuss, in Part III, an extension of 
our equations to include the non-isotropic components of the pressure tensor, 
eliminating the question of charge separation by considering the case M =m. 

In the zero temperature case, failure to find solutions for x > 2 is explained 
by the appearance of «looping » orbits which will make all solutions dependent 
on the time. When 6 >0 the picture is more complex, since a fraction of 
the orbits will always be looping. We simply say that time-independent solu- 
tions do not exist for B, so small that x > «+, and this is not surprising be- 
use in the absence of a magnetic field such solutions do not exist at all. 
Waves will broaden continuously until finally collisions establish a steady 
state with characteristic length of the order of a mean free path. In a given 
plasma we find that waves of the type described in this paper can propagate 
with any arbitrary velocity greater than the sound speed, provided B, exceeds 
some limiting value. We have placed the burden of maintaining the steady 
state upon the magnetic field. It must be strong enough to accomplish the task. 


Part of the work reported here was started while one of us (A.B.) enjoyed 
a sabbatical leave of absence in 1958 at the Istituto di Fisica Guglielmo Mar- 
coni dell’Universita, Roma. We wish to express our thanks to Professors 
Enrico Persico and EpoArpo AMALDI for their generous hospitality at the 
Istituto, and to the J ohn Simon Guggenheim Memorial Foundation and the 
U.S. Fulbright Commission for fellowship grants which made possible the 
leisurely prosecution of these studies. In addition, we are indebted to 
Drs. THEODORE C. NorTHRoP and MARVIN MITTLEMAN, of the U.C. Lawrence 
Radiation Laboratory, and to Drs. HAROLD GRAD and CATHLEEN MORAWETZ, 
of the N. Y. U. Institute of Mathematical Sciences, A.E.C. Computing and 
Applied Mathematics Center, for many illuminating discussions on the subject. 
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RIASSUNTO (*) 


Consideriamo una estensione infinita di plasma a bassa densità, completamente 
ionizzato, magneticamente immobilizzato in un campo costante ed uniforme di indu- 
zione magnetica. In assenza di collisioni, supponiamo chei movimenti degli ioni e degli 
elettroni abbiano luogo in piani perpendicolari al campo magnetico. Poi passiamo 
ad un sistema di riferimento di « urto » che si muova a velocità costante in direzione nor- 
male al campo magnetico, e ricerchiamo la classe di soluzioni non banali, indipen- 
denti dal tempo, ad una dimensione e autocongruenti delle equazioni di Maxwell e 
delle equazioni del moto delle particelle cariche, cioè, studiamo la propagazione degli 
impulsi trasversali e delle onde la cui lunghezza caratteristica sia molto minore del 
cammino libero medio. Facendo delle opportune semplificazioni, arriviamo ad un sistema 
di equazioni in cui si trascurano gli urti, si presuppone la neutralità della carica, e sia 
gli ioni che gli elettroni si comportano come fluidi bidimensionali che seguono le rela- 
zioni isoentropiche con y=2 e presentano tensori di pressione rigorosamente isotropici. 
Con queste semplificazioni e presupposti, il sistema di equazioni può essere risolto 
completamente con una semplice quadratura numerica. Si ottengono impulsi simmetrici 
e onde periodiche, con lunghezze caratteristiche dell’ordine della profondità dello strato 
superficiale del plasma, ossia del raggio giromagnetico medio per particelle che si 
muovono alla velocità di Alfvén. In questo articolo, parte I, descriviamo gli impulsi 
unici, che sono le sole soluzioni che soddisfino le condizioni del plasma non perturbato 
in posizione anteriore al treno d’onde. Per una completa determinazione del problema 
si richiedono due parametri fondamentali: il numero di Mach di Alfvén, che dà una 
misura della velocità di propagazione, ed ilrapporto B fra la pressione totale del plasma 
e la pressione magnetica, che determina lo stato iniziale del plasma. Si dimostra che 
per un dato f esistono soluzioni stazionarie solo per un limitato campo di valori della 
velocità «. 


* q . . 
©) Traduzione a cura della Redazione. 
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The Analytic Properties of Perturbation Theory - I (*). 


J. C. POLKINGHORNE 


Trinity College - Cambridge 
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Emmanuel College - Cambridge (**) 


(ricevuto il 26 Ottobre 1959) 


Summary. —- A general method is given for locating the complex singu- 
larities of the contributions from Feynman diagrams regarded as functions 
of the external scalar products. The method is illustrated by application 
to the third-order vertex function. 


1. — Introduction. 


In recent years it has been noticed that many functions of physical interest 
possess powerful properties when regarded as functions of complex variables. 
As part of the programme for exploiting these properties investigations have 
been made of the analyticity of the terms of perturbation theory regarded 
as functions of energy, momentum transfer, etc. The earlier work on this 
subject (5) was largely concerned with the singularities occurring for real 
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values of the arguments. However, since the motivation for much of this 
work is an interest in Mandelstam’s conjecture (7) on the analytic properties 
of the scattering amplitude, it is also essential to determine the positions 
of the complex singularities. The very informative example of the third-order 
vertex function has been completely discussed by OEHME (8) using the explicit 
expressions obtained by KALLEN and WIGHTMAN (°), but no general method 
has so far been given for investigating these complex singularities. 

It is our purpose in this paper to develop such a method (*). It is based 
on an approach used by EDEN (1) for a discussion of the analytic properties 
of perturbation theory, but the method is applied in our case not to the inte- 
erations over internal momenta but to the integrations over Feynman para- 
meters that remain when symmetric integration has removed these internal 
momenta. The final form of the answer is the same as that given by LANDAU (°) 
for the real singularities. Indeed he seems to state that this form should 
determine the complex singularities also, but no discussion is given. 


2. — Funetions defined by perturbation theory. 


The contribution from an arbitrary graph in perturbation theory is of 
the form 


si 
(1) Lim fax sit da, far soo GEIR, AGE Le = de = 2 5 
ao [E(0;; hz; Pa) + tel” 
where 
(2) Bai; ky; Pr) = D a. (G5 — mi). 


i=1 


3 Pr ne "at CAN i i 
The variables «,...%, are the Feynman parameters used to write the contri- 
butions from the » lines of the graph in the form of a single denominator. 
Jach line is associated with momentum g; and mass m,. The variables k, 


n (*) While this work was being written up we were shown a preprint by Mr. J 
DARSKI in which essentially the same method as that given here is used to atoona. fine 
special case of contributions from graphs with only one internal loop. Mr. TARSKI 
does not, however, discuss the peculiar complications that arise from the A 
integrations over Feynman parameters (see Section 4). 
(7) S. MANDELSTAM: Phys. Rev., 112, 1344 (1958). 
(8) R. OEHME: Nuovo Cimento, 13, 778 (1959). 
) G. KALLÉN and A. S. WIGHTMAN: Mat. Fys. Skr. Dan. Vid. Selsk., 1, no. 6 (1958) 
) R. J. EneN: Proc. Roy. Soc., A 210, 388 (1952). RA si 
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are some set of independent internal momenta associated with the graph. 
They may be conveniently chosen by selecting à set of independent closed 
loops specifying the structure of the graph. The variables p, are the external 
momenta. The case of all Bose particles is discussed for the sake of simpli- 
city: no new difficulties arise from the introduction of Fermi particles. 

It is assumed that some procedure of regularization has made the integral (1) 
convergent for the k-integrations. These are now performed by the technique 
of symmetric integration giving a function of the external scalar products 
Dir = Di Pr, Of the form 


1 


(3) fn) = Lim fa de. sai eg a nr 


e>0+ 


0 


Numerical multiplicative factors have been omitted. The function Ff’ is ob- 
tained from the function F by eliminating the k; by use of the equations sum- 
marized in vector form by 


OF 


== = D = My oc We 
Ok; , 1 9 5} 


(4) 


The fact that F' is obtained from F by taking its stationary values with respect 
to the integration variables has been particularly emphasized by MATHEWS (ae 
The function g(2;) is a simple algebraic function of the «’s which is relevant 
for our considerations only if it vanishes at some critical value of the os. 

The presence of the ie term in (1) permits the use of symmetric integration. 
It tells us that for real values of p,, the physically important function is that 
obtained by letting the masses associated with the internal lines tend to the 
real axis from the lower half plane. Since the method we discuss below is 
easily extended to discuss analytic continuation in these internal masses we 
shall omit the ie term and consider the function of the complex variables 2;, 
defined by 

1 

7 ff (ox) 0(% + el) 
(5) f (2x) ‘| ogee daa wa, ar de 
. 


0 


In order to discuss the singularities of this function it is convenient to establish 
a lemma suggested by the work of EDEN (?). 


(1) J. MATHEWS: Phys. Rev., 118, 381 (1959). 


oO 
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3. — Lemma. 


Consider a function f(£) defined by 


(6) vita) |P, E)du . 


o 


E may stand for several complex variables. Then f(£) is certainly a regular 
function of € in the domain Z such that for all ÊeZ, Flu, È) is a regular 
function for all « lying on the contour C. We 
wish to show that f(Ë) can be continued along 
a path from Z to any point È” lying outside Z 
provided this path does not pass through certain 
exceptional points ¢*, defined below. 

The functions ÆF(u, ¢) with which we shall 
be concerned are all such that in general as è 
varies along a path from some point ¢'<Z to 
the point ¢” then the singularities of F(u, È) 
describe paths in the «-plane in a continuous 
fashion. No singularity crosses C until we 


Fig. 1.— Paths of singularities 


leave Z. When we have reached some point 
in the w-plane. 


E outside Z the singularities, ,(Ë), may have 
described some such paths in the u-plane as 
are depicted in Fig. 1. The contour €, is such that the function 


(7) AOS | F(u, ©) du, 


(pt 
is regular at Ê—£,. Successive points along the path define successive con- 
tours in this way with the property that the contours define functions regular 
at the point in question and coinciding in some neighbourhood with the fune- 
tions associated with neighbouring points along the path. The contours are 
obtained by continuous deformation of the original contour C away from the 
direction of the advancing singularities. | 
continuation of f(£) to the point è”. 

The method may, however, fail for one of three reasons: 
i) A singularity may pass through one of the end points of € in which 


case no useful deformation can avoid it. We shall eall this case 
singularity. 


The associated functions provide a 


an end-point 
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ii) Two singularities coming from opposite sides of the contour may CO- 
incide for the same È, say ¢*, as in Fig. 2. The contour is then « pinched » and 
no further useful deformation is possible. We shall refer to this as a case of 
coincident singularities, though it is important to notice that only singula- 
rities tending to coincidence from opposite sides of the contour give rise to 
exceptional points. The situation illustrated by Fig. 3 is quite harmless. 


UÈ) = Uz(6) 


Fig. 2. - Coincident singularities. Fig. 3. - Harmless coincident singularities. 


iii) For some values of ¢ one or more singularities may recede to infinity. 
This is best discussed by employing a transformation that brings the point 
at infinity into the complex plane. In general our integrands are well-behaved 
at infinity and only obviously coincident singularities give rise to exceptional 
points. These singularities are not, therefore, separately classified in what 
follows. (A particularly striking case of this behaviour occurs when for some È, 
F(u, 6) becomes independent of u. All singularities have receded to infinity 
unless Fu, £) is itself singular at this point when in fact they cover the 
whole plane). 

The exceptional points ¢* are defined as the class of points that give rise 
to singularities in f(¢) in one of these three ways. 

The value of the function at the point ©” that we obtain by this method 
of continuation will depend upon the path of continuation chosen. That is 
to say, the function f(C) is in general multi-valued. It is clear from the freedom 
that exists in the choice of successive contours C, that if two paths in the 
É-space are chosen $0 related that one can be obtained from the other by a 
continuous deformation that never passes through an exceptional point ©, 
then they will give the same value. 

Finally, in our discussion so far we have assumed that the limits of our 
contour 0 are independent of ¢. In the cases to which we wish to apply the 
lemma this is not so. One end point is fixed but the other depends on ¢. 
However, its dependence is in terms of an analytic function (indeed a linear 
function) and so the conclusions of the lemma are not disturbed. 
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4. — Application to multiple integrals. 


We now wish to apply the ideas of the last section to the multiple integral 
defined by equation (5). The ô-function must first be eliminated and so we 


consider the integral 


if LOL Ori een ar Og à 
) 1 1 
(8) f (2sx) =“ da, da, SRO GK, LP": 2 jn) "2 i 
0 0 0 
F’ is obtained from F’ by eliminating the variable «, by the substitution 
An =1— %,—...—a,4. Although the integral (8) selects a particular ordering 


of the x°s we shall find that our final answer can be cast into a form symmetric 
in the n os. 

It is tempting to suppose that f will be singular only if each integration 
in (8) has either an end point singularity or a coincident singularity. This 
in fact turns out to be the case but the argument leading to it is somewhat 
involved. For any particular integration in (8) the yet-unintegrated «’s com- 
bine with the z,, to form the ¢ of the lemma. Moreover, after the «,_, inte- 
gration has been performed the successive integrands are in general multi-valued. 

The method that must be used to deal with this problem is sufficiently 
illustrated by the case of two integrations. Consider the function defined by 


(9) È = [à fa + i 
di È ACTUS 


= and B are the appropriate contours of integration for « and 6, and F(x, B; È) 
Is assumed to be regular in the finite plane. The singularities of the integrand 
are therefore just the zeroes of F. Once again $ may signify a number of 
complex variables. | 

A domain Z is assumed to exist such that for SPAZI ZU 
«E A, BEB. The points outside Z may be singular points of (£) or they may 
Da points lying on the specific cuts defined by the integral representation (9) 
since this representation clearly defines a one-valued function. The purpose 
of this investigation is to Separate these two cases. 

This will be done by defining a continuation of our function by deforming 
oy contours A and B. In the «-plane we shall have to avoid the points 
x (5, Ë) which are such that F(a’, 8,6) =0. By the lemma of Section 3 this 
can always be done to define a (multivalued) function regular at all f and È 
except the exceptional points (6'(2), £) which are such that | 
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either: i) F(a, B'(C); 6) =90, where x is an end point of A (end point sin- 
vularities) ; 


or ii) F = 0F[0ax=0, some a, and the two coincident singularities that 
this implies approach each other from opposite sides of the «contour 
(coincident singularities). 


It will therefore be necessary to deform f to avoid the points B'(£) as È varies 
along its assigned path. 

Until ¢ leaves Z no deformation of A and B is required and so in parti- 
cular no B'(Ë) lies on B. Thereafter the contour B is deformed in a continuous 
way to avoid the points p’(¢), in the manner of Section 8. This continuous 
deformation defines (in a way that is inessentially non-unique) a path linking 
each point on the original contour B with its corresponding point on the de- 
formed contour B,. In this way the original path in the £-plane defines a 
series of paths in the combined f-$ space. These paths select the branch 
of the multi-valued function obtained by performing the g-integration. This 
function is analytic in B and & for every point on the contour B,. Moreover 
the continuity properties of our functions and of the paths defined above 
ensure that the function associated with a particular B and € coincides in 
some neighbourhood with the functions associated with neighbouring B and ¢ (*). 
Thus the contour B, is so arranged that the B-integrand is regular along it 
and so the B-integration yields an analytic function of ¢, which depends, of 
course, on the path of continuation in the £-plane. The method fails, however, 
if B, can not be deformed to avoid the points B'(£). This will happen if one 
of these points coincides with 6), an end point of B, or if two points f tend 
to coincidence from opposite sides of B,. This only happens for certain values 
of £ and these points are the singularities of g. Thus we may say that the 
singularities of g occur at the points ¢* which are such that 


either: 1) Fos Bo; C*) =9 (an end point singularity in x giving an end point 
singularity in f); 


or ii) F(a, Bo; 6*) = (OF 0x) (x, Bo; £) = 0, and the coincident singularities 
in x «pinch » the contour (a coincident singularity in « giving an 
end point singularity in f); 


or ili) Fo, 8; 6*) = (OF /0B) (a, B; ¢*)=0, and the coincident singularities 
in B «pinch» the contour (an end point singularity in « giving 
a coincident singularity in f); 


(*) Singularities receding to infinity require a special discussion here which is best 
made by transforming the point at infinity into a point of the finite plane. 


or 
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or iv) F(a, B; 6*)=(0F} ea) (a, B; ¢*) = (OF /0B) (a, B; ¢*)=0, and the coinci- 
dent singularities in both x and f « pinch » their respective contours 
(coincident singularity in « giving a coincident singularity in B). 


We note that the specification i) to iv) produces a scheme independent of the 
relative x, B order of integration, although that order played a part in our 
definition of the continuation. The final answer is what one expected: there 
must be an end point singularity or coincident singularity in each integration. 


5. — The singularities of perturbation theory. 


The ideas of the previous section can now be applied to the perturbation 
theory integral, equation (7). Each integration must either give rise to an 
end point singularity or a coincident singularity. Thus we must have 


(10) (Ge 25x) Ml) ; 


together with 


either: i) a, =0 (lower limit end point) or 
% = (upper limit end point); 


À ve ap ; 7: 5 pee s 
0) li) CF"/04x; —0 and the resulting coincident singularities in the x, 
integration «pinch » the contour; 


WOE =I ong A 
These conditions may now be written in a more symmetrical form. By 


the use of the condition 


(11) HE... x, —= 1. 


ARE HT 3 ; : 

P"(x;; &x) in (10) may be equivalently replaced by the function F'(x;; 2;,) 
which is homogeneous of degree one in the «’s. Then the derivatives of F” 
are just the derivatives of F’ subject to the condition (11) and may be written 


12 
(12) rae + 1, 1=1,...,n—1, 


where 7 is a Lagrange multiplier. 


| However, Leibnitz theorem for homoge- 
neous functions implies that 


(13) P'(x e» ) == x Xi oF" 


iy “ik 
i=1 dx; 
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and so the vanishing of Y" (or equivalently Æ) implies that under the con- 
ditions (10) À—0 and these conditions may be rewritten as 


either 6 == (Ws 
(14) e OF’ Vel 1 
or — 0, 
OX; 


where the symbol ~ is used to denote the fact that two singularities must 
not only coincide but « pinch » the contour. 

When account is taken of the derivation of F' from F by equations (4) 
and (11), the conditions (14) are seen to be just the generalization of the con- 
dition given by LANDAU (°) for the real singularities and derived by him in 
a different manner. From equations (2) and (4) it follows that 


OF" m OF Ok; 


5 el ee EN 
(15) da, (qi mi) + 2 ah, da, 


(qi — mi) , 
and so the conditions (14) can be interpreted that either a line is omitted from 
a diagram or its momentum vector is put on the mass shell. This leads to 
the dual diagram analysis of J. C. TAYLOR (5) and LANDAU (°) for the real 
singularities. 

The exclusion by the condition of « pinching » in (14) of harmless coinci- 
dences of the type illustrated by Fig. 3 means that not all the points satisfying 
the weaker conditions: 


either ia AU 
(16) oF A160, 7; 
or — = 0, 
Ca; 


are singular points of f on the physical sheet defined by the integral repre- 
sentation (8). However, when the funetion is continued into all its sheets 
then the equations (16) give the possible singular points. This is because this 
continuation may be achieved by considering the totality of contours joining 
the given end points and thus every coincident singularity « pinches » some 
contour. An example of how singularities can appear and disappear on given 
sheets is provided by the third-order vertex function discussed below. Since 
we are concerned in physics with particular sheets it is necessary to supple- 
ment the simple equations (16) by an analysis of whether the coincident sin- 
gularities « pinch » the contour (*). This is simple in principle, but tedious 
in practice, for any given diagram. 


(*) A most interesting illustration of this has been given by Mr. Tarski in his dis- 
cussion of the fourth-order scattering diagram. 
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We may note that the consequences of (16) find a particularly simple expres- 
sion in the case when the diagram under discussion has only one closed loop. 
P'(%;; &x) is then a polynomial of second degree in the x; and linear in the 
21 divided by a linear term in the «;. The linear term is actually the sum of 
the x; and so is just equal to one. The important term is the second degree 
polynomial which may be written in matrix form as 


(17) DA (eas 2. 


The condition that this form together with all its first x-derivatives vanishes. 
for some values of the «°s is that its discriminant vanishes: 


(18) det [A;,(¢)] = 0. 


The other singularities, corresponding to the vanishing of one or more «’s, 
are then given by the vanishing of the principal minors of det [A;;] obtained 
by omitting the appropriate rows and columns. J. C. TAYLOR (12) has already 
conjectured that the location of singularities might be given by these deter- 
minantal conditions in simple cases. 


6. — The third-order vertex function. 


We conclude by applying our method to the third-order vertex function, 
or rather to the differentiated third-order vertex function which has been 
particularly discussed by OEHME (5). To avoid irrelevant algebraic compli- 
cation we will assume the equality of all the masses of the internal particles. 
In the notation of reference (5) we are concerned then with the function 


= Fax dB dy Ô(1 — NE) 


(19) 2,7 
MOV EE = afZ 


This is an example of a function associated with a single loop diagram and 
the conditions described by equation (18) et. seq. may be used to find the 
possible singularities. The vanishing of the determinant gives 


(20) Z,(m?Z, + Z(Z — 4m?)] = 0, 


(©) J. C. TAYLOR: private communication. 
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and the principal minors give 


(21) D(A Any , 
(22) Z,(Z,— 4m?) = 0. 


KALLÉN and WIGHTMAN have integrated (19) explicitly (*) obtaining 


a 
m2Zs + Z(Z — 4m?) 


RS ZZZ ma) 
— (oy A1 — de dj — Ÿ < 5 sh 
(Z— 4m2)3 ° 2m? or) 


Z3— 22 Vi an à la | | 
ZA oye Cie La i) 


(23) {(Z, Z;) == 


Clearly the singularities of (23) lie on the curves (20)-(22). 

It is also important to know when they occur on the physical sheet. The 
physically important function is that branch obtained by letting Z and Z; 
tend to the positive real axis from above. While Z and Z; both have the 
non-vanishing imaginary parts of the same sign the integrand in (19) cannot 
vanish in the region of x and B integration (*). No contour deformation is 
required therefore until the real axis is actually reached. The physical function 
will have singularities at real Z and Z, only if on reaching the real Z and Z; 
axes from above end-point or coincident singularities appear on the relevant 
parts of the « and f real axes. 

For example the doubly coincident singularities associated with (20) are 
found to occur at 


(24) ea 


Since we integrate only over real positive values of x and B such that «+6 <1, 
if Z is real and 


(25) Th, PHO à 

this singularity will not occur for the physical function. If, however, Z is 
real and 

(26) Z > 2m? 


(") We assume now the y-integration done to remove the 6-function. 
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the corresponding values of « and f lie on the contours and a trivial calculation 
shows that they do indeed « pinch» them. Thus we are able to understand 
the sudden appearance on the physical sheet of the anomalous threshold dis- 
covered by KARPLUS, SOMMERFIELD and WICHMAN (ye 


RIASSUNTO (*) 


Si espone un metodo generale per la localizzazione delle singolarità complesse dei 
contributi dei diagrammi di Feynman considerati come funzioni dei prodotti scalari 
esterni. Il metodo viene illustrato con l'applicazione alla funzione vertice del terzo ordine. 


(*) Traduzione a cura della Redazione. 
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(La responsabilita scientifica degli scritti inseriti in questa rubrica è completamente lasciala 
dalla Direzione del periodico ai singoli autori) 


The Effect of Multiple Scattering in the Photoproduction 
of Charged Mesons at Deuterium. 


N. MACDONALD 


Department of Natural Philosophy. The University - Glasgow (*). 


(ricevuto il 14 Dicembre 1959) 


This note reports briefly on a calculation of the multiple scattering correc- 
tion to the impulse approximation treatment of the photoproduction of charged 
pions at deuterium. One of the original aims of this work was to shed light on the 
process of photoproduction of negative pions at free neutrons. It is now apparent 
that a more fruitful approach to this problem would be along the lines of recent 
work by Crew and Low (1) on unstable particles as targets in scattering experiments, 
and such work has been reported by ANDERSON et al. (?). However there is still some 
interest in the multiple scattering approach because of its relation to the problem 
of x-D scattering, recently treated by BRANSDEN and MooRHOUSE (3), ROCKMORE (4) 
and DE ALFARO and STROFFOLINI (5). Also recent experimental work by LAND (°) 
shows that the impulse approximation gives a reasonable description of the process 
in circumstances for which our work predicts that the multiple scattering correction 
is very small. 

The impulse approximation has been used for this process by various authors (7) 
and experimental results have been analysed using their work (98). In this approx- 
imation (°) the transition operator for the process is taken as the sum of transi- 


(*) Present address: Atomic Weapons Research Establishment, Aldermaston, England. 

(1) G. F. CHEW and F. E. Low: Phys. Rev., 113, 1640 (1959). 

(2) J. D. ANDERSON, D. C. GATES, T. L. JENKINS, R. W. KENNY and W. P. SWANSON: Bull. 
Am. Phys. Soc., 4, 356 (1959). 

() B. H. BRANSDEN and R. G. MoorHouse: Nucl. Phys., 6, 310 (1958). 

(4) R. M. ROCKMORE: Phys. Rev., 105, 256 (1957). 
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(8) R. H. LAND: Phys. Rev., 113, 1141 (1959). 

(7) S. MacHIDA and T. TAMURA: Progr. Theor. Phys., 6, 572 (1951); Y. Sarro, Y. WATANABE 
and Y. YamacuoHI: Progr. Theor. Phys. 7, 103 (1952); G. F. CHew and H. W. LEWIS: Phys. 
Rev., 84, 779 (1952); M. Lax and H. FESHBACH: Phys. Rev., 88, 509 (1952). 

(®) R. S. WHITE, M. J. JACOBSOHN and A. G. ScHuLIZ: Phys. Rev., 88, 836 (1952); D. C. HAGE- 
MAN, K. M. Crowe and R. M. FRIEDMAN: Phys. Rev., 106, 818 (1957). 

(*) G. F. CHew and G. C. Wick: Phys. Rev., 85, 636 (1952); G. F. CHEW and M. L. GOLD- 


BERGER: Phys. Rev., 87, 778 (1952). 
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tion operators 7, and T, for photoproduction at nucleons 1 and 2. The effect of 
of nuclear binding is ignored, except in so far as it determines the wave function 
for the initial and final states. The interaction of the meson and the two nucleons 
is also ignored. [We use the term impulse approximation to refer to the treatment 
which ignores both these factors, although strictly it is the failure to include the 
effect of binding that makes it an impulse approximation.] We continue to neglect 
the effect of binding and deal with the meson-nucleon interaction by allowing for 
successive scattering at alternate nucleons. This has been done previously for elastic 
production of neutral pions at deuterium by CHAPPELEAR (1°) and at helium by 
STOODLEY (11). | | 

We use the same, method, essentially, as CHAPPELEAR. Writing the matrix ele- 
ment, between states of one nucleon and one meson, of the transition operator t; for 
the interaction of a meson with nucleon 7 as 


(1) (q\t qi) = &;b(q , &)q1°2 XP [i(qe — qi) Tél ; 


this method involves the approximation that 


= 2ri0(45)b(4 > de) 0(Qr> di) EXP [ig PR]. 


È fe exp [igE]b(4; 4) (4, 1) 
(dz) — @(q) + te 


— oO 


In the above a, is a function of isotopic spin, (the effect of ordinary spin is ignored). 
qi>qz are the initial and final meson momenta; 

r; is the position of nucleon 7; 

R= |r,—r,|; 

(9) is the meson energy for momentum q: and 


de is the meson momentum on the energy shell. 


| If we take b in the form 6(q,, q2)=c(q,)d(q2) we can evaluate the correction 
without using the approximation (2). The form of b used, which is taken from 
work of VELIBEKOV and MESHCHERYAKOV (!2), resembles that used in ref. (5) for 
the 7 - D problem. In this case we find that the approximation (2) is fairly good. 
The use of the approximation (2) enables one to use (qo\t:\qi) in the form 
@b(4z; Az): 2 ex? [i(q, — qz):r;] alone, without requiring knowledge of b off the 
energy shell (*). We take b(q,, dz) as È 


27 


3 en 3 LES o 
(3) n exp [1033(Gz) ] Sin d33(4£) - 


(1°) J. CHAPPELEAR: Phys. Rev., 99, 254 (955): 

‘a & AD) ag STOODLEY: unpublished thesis (1957). 

Sena È n A. MEDS à Dokl. Akad. Nauk SSSR, 105, 941 (1955). 
pproximation does not however correspond to taking scattering on the energy shell 


onl , as has sometimes been st ated DRI LL id VERI E 13 isti 1 etwee hese a proxima t . 
Y BP 2 Ha. anc ‘E ( ) dis ngu sh b WV n i i 
3 $ È > p 1008 
( DST D. DRELL and L. V ERLET: F[ hys. Rev., 99, 849 (1 955) 
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with empirical values of the (3, 3) phase shift. For the operator T; we take the 
electric dipole and dominant magnetic dipole (again corresponding to (3, $) state 
of meson and nucleon) operators of CHEW and Low (14). 

Results were obtained for 300 MeV ys. The impulse approximation energy 
spectrum of mesons at a given angle is made up of two peaks. One is about the 
energy of the meson produced at this angle from a free nucleon, and dominates the 
energy spectrum except at small forward angles. On this peak the multiple scat- 
tering correction is about —5%. The second peak occurs near maximum meson 
energy, and is a consequence of the interaction of the residual nucleons when they 
‘have low relative momentum. We in fact include the nucleon interaction in the 
S-state only, and find that on this peak the correction rises to — 20%. Now this 
situation is the nearest in our process to elastic production, so the larger correction 
is consistent with the fact that Chappelear found a —40% effect in elastic 
r° production. 

In general comparison with experiment is hampered by the fact that most expe- 
riments have been done not with a monoenergetic y beam but with a bremsstrahlung 
spectrum. However recent work by LAND (°) uses a beam with energy (292 +8) MeV. 
He detects mesons at 120°, which means he is dealing with the free nucleon peak, 
and in agreement with our result he finds that the impulse approximation is adequate. 

In the case of x - D Scattering, DRELL and VERLET (#) estimate that double 
scattering provides the major part of the multiple scattering correction. This is 
in agreement with Rockmorn’s results (4) but not with those of pe ALFARO and. 
SrROFFOLINI (5). The corresponding process in our case is photoproduction followed 
by one scattering. In our formalism we cannot estimate the contribution from this 
process separately. However the contributions from photoproduction followed by 
an odd or an even (non-zero)number of scatterings appear separately, and are of 
comparable magnitude. The former would be expected to predominate if the most 
important process were photoproduction with one scattering. 

This work was carried out while the author held a maintenance grant from the 
Department of Scientific and Industrial Research. It forms part of a thesis sub- 
mitted to the University of Glasgow. The problem was suggested by Prof. J. C. GUNN, 
who is thanked for his advice and encouragement. Helpful discussions of the mul- 
tiple scattering formalism with Dr. K. D. C. SroopLEY are also acknowledged. 


(4) G. F. CHEW and F. E. Low: Phys. Rev., 101, 1579 (1956). 
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Cross Sections of Reactions Produced by High Energy Neutrino Beams. 


N. CagrBBo and R. GATTO 


Istituto di Fisica e Scuola di Perfezionamento dell Università - Roma 
Istituto Nazionale di Fisica Nucleare - Sezione di Roma 


(ricevuto il 28 Dicembre 1959) 


1. The expected availability in the next future of high intensity beams of 
high energy particles makes it appealing to think of possible experiments to detect 
reactions initiated by high energy neutrinos (1). Neutrinos are known to have 
only weak interactions and their cross sections at low energy are indeed very 
small (2). Neutrino cross-sections are however suspected to grow-up very fast with 
energy (>). 

One would argue that, in a reaction such as y+p->et+n the phase space 
volume will eventually grow-up with the square of the C.M. energy, while the 
matrix element derived naively from the four-linear @ coupling will remain constant, 
thus giving a cross-section increasing with the square of the C.M. energy. This 
argument is obviously wrong at least because of the presence of form factors and 
of new couplings induced through the virtual fields of particles which interact strongly 
with the nucleons, such as pions, etc. 

In listing the possible reactions which can be produced from high energy neu- 
trinos one has to note that many reactions are forbidden by the lepton conservation 
law. Their absence would actually constitute the most direct proof of lepton con- 


(1) We have learned that various physicists, namely Prof. AMALDI, Prof. BERNARDINI, Prof. Coc- 
coNI, Prof. LEDERMAN, Prof. LEE, Prof. PREISWERK and perhaps many others are thinking at 
present on the possibility of such experiments. 

(AC Di Cowan, F. REINES, F. B. HARRISON, H. W. Kruse and A. D. Mc GUIRE (Nature, 
178, 446 (1956)) gave a cross-section 6-1074 for he reaction V+p — n +e*. 

The experiment was possible because of the high intensity of neutrinos from a thermal neutron 
reactor. 

; (8) Da have altcady been theoretical discussions on the importance of weak interactions at 
high enersies [see, for instance, D. I. BLOKHINSTEV: On the Limits of applicability of Quantum 
Electrodynamics, F.148, Joint Institute for Nuclear Research, (1958)] 
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servation (4). The allowed reactions are 


y+p—et+n VERI Rep 
y+p>et+n+n°0 ytn—-et+n+n- 
eh Pat 
v+p—oe+p+tr* yv+tn->e-+n47*t 
Pr? 
v+poet+A° 
v+p—et+20 yt+n—->et+ x7 


[v+n>e-+27], 


and so on. 

In all the above reactions the electron can be substituted by a muon of the 
same charge. We have enclosed in brackets the last reaction v+n se >, be- 
cause it is one of those reactions which would be forbidden if some of the present 
ideas on the interactions of strange particles are right (°). Such a reaction is forbidden 
if one excludes a weak strangeness-non-conserving current with a strangeness change 
opposite in sign to the sign of charge. Such a forbiddeness is directly related to 
the suggested forbiddeness of 2° >n+e*t+v. Most of the listed reactions can only 
occur incoherently on a nucleus; only reactions in which the same kind of nucleon 
appears in both sides of the equation (such as v+p—>e +p+7", ete.) can also 
proceed coherently. 


9. — To be definite let us consider the first of the listed reactions 
(1) v+p—et+n. 


We shall first derive the expression for its cross section in terms of the nuclear 
form factors. In the next section we shall try a definite estimate of the cross- 
section by adopting the non-renormalization hypothesis for the vector strangeness 
conserving current and making use of extrapolated Stanford data on the nucleon 
form factors, together with some guess on the relevant axial form factor. The 
general expression (7), (8) that we derive is rigorous in the V— A theory (°) under one 
important assumption: that the behaviour of the strangeness-conserving vector 
and axial currents under the product operation of charge conjugation and charge 
symmetry is the same as that of the nucleon terms, yy,t+# and yy,7;T+Y respec- 


(4) The absence of the reaction ¥ +37Cl > A +e7 has been discussed by R. Davis: Bull. Am. 
Phys. Soc., 1, 219 (1956); and by C. O. MBUHLHANSE and $. OLESKA: Phys. Rev., 105, 1333 (1957). 
Evidence against double @ decay is discussed for instance by C. L. Cowan and F. REINES: Phys. 
Rev., 106, 825 (1957). The study of neutrino reactions will also offer a conclusive test of the hypo- 
thesis that the neutrino emitted with the electron is the same as that emitted with the muon. 

(®) A detailed discussion can be found in R. Garro: Fortschr. d. Phys., 7. VAT (1959). 

(°) R. P. FEYNMANN and M. GELL-MANN: Phys. Rev., 109, 193 (1958); R. E. MARSITAK and 
E. C. G. SUDARSHAN: Phys. Iev., 109, 1860 (1958). 
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tively. Such an assumption is certainly verified for the vector current under the 
non-renormalization hypothesis. If it werenot verified the expression for the cross- 
section would become more complicated and involve two more form factors. Of 
course, concerning the numerical estimates that we give, we could say very little 
at present on how much they would be changed without the non-renormalization 
hypothesis for the vector current or for the case of more complicated behaviour of 
the currents under charge symmetry. Considering the ambiguities on the form 
factors any such discussion would be only academic at present. The matrix element 
for (1) is given by 


() VIG a 8 + pe GD, 
where 

(3) ju = (rh | RP (0) in (0) [p> ; 
and 

(4) 1, = (00) 77 3(1+ ys) 0(e*)) . 


In (3) j{P(w) and jf/(x) are the strangeness-conserving vector and axial currents 
respectively. In (3) G is the vector coupling constant for 6-decay, as defined in (5), 
G~1-10° WM, M being the nucleon mass, and », n, p, e are the four-momenta. 
Under the above explained assumption about the behaviour of such currents under 


the charge symmetry operation, j, takes the form (7). 


“FB y(k2) opp ky] (P) , 


2M © 


(5) Ju = (n) CHR?) yyys + tH (k?) hays + F(k?)y, + 


where H,, H,, F; and F, are form factors depending on k?, with ky=(p —n),, pis 
a numerical constant, and u(n), u(p) are free particle spinors. The definition of the 
form factors in (5) is such that 


(6) AO HO Re 2) 


where G, is the axial B-decay coupling constant. The term proportional to H, does 
not contribute if the electron mass can be neglected. From (2), (4) and (5) one 
derives the following expression for the differential cross-section for reaction (1) 
in the laboratory frame 


Y2 


¥ G 
(7) (0) d(cos 0) = — E2X(0, €) d(cos 6) , 
IT 
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with 


26) RON 
(sie 710,6) (14 cnt) Loos 2 À wane) + 


Me HR) LEE te + 2 Con es = 
tao M? 1+ 2(&/M) sin? 0/2 


4 sin? (0/2) 


— Hl?) (F(R) + MNT Ten sine (0/2) 


(us + &? sin? 5) i 


In (7) and (8) & is the neutrino energy in the laboratory system, 0 is the scat- 
tering angle in the laboratory system, and M is the nucleon mass. For the process 
v+n > e-+p the cross section has the same form except for a change of sign of 
the interference term (the last term inside the bracket in (8)). 

One verifies easily that with constant form factors the total cross-section in- 
creases linearly with the laboratory energy &, when E > M, i.e. it is proportional 
to the square of the C.M. energy. 

If, on the other hand, the form factors are essentially equivalent to a cut-off 
for k2—a2 the total cross-section at high energies increases as 


Le Cee CCM energy) Aq , 
where Aw is the «allowed » solid angle, defined by <a, which gives 
Aw cc (C.M. energy). 


The cross-section will then tend to a constant value rather than increase indefini- 
tely. Of course, any definite statement about the high energy behaviour of the cross- 
section is impossible, as it is bound to a knowledge of the form factors for high 
momentum transfer. A study of u-capture reactions can give information on the 
form factors for k? ~(100 MeV)?. If one uses the non-renormalization hypothesis 
the form factors F, and F, can be taken from electron scattering experiments. 


3. — To make an estimate of the cross-section and of its behaviour with energy 
we have used the non-renormalization hypothesis for the vector strangeness-con- 
serving current. In this case u= Mp — Un = 3.7; and F,(k?) and F,(k?) can be taken 
from the analysis of electron and muon scattering on nucleons — it can be easily 
seen that they are given as differences of proton and neutron form factors. The 


Stanford results suggest 


if 
(10) a eee 


In the absence of any information on H,(k?) we make the simplest choice 


A 
A 1 


x (1+ Ka} 
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The total cross-section from (8), (9), (10) and (11) can then be derived analytically 
and it takes the form 


12 


à 2 
Ga? Gi ae H? ~ M G4 
(12) Giga = Lao É si ] + AP(a) = | 1+ i 


@ 
i G 


y2 7 
M? (6 weit Gk | 4 SI 5 GA M, 
Al I) t AS (OF 1682 fl Sta 21 n) = 4 G2 se 


a G A 


a 
A I 4 (1) sE A2 n 
marg IC LA Care | 
In (12) 
AO(&) = 1—(1+ a), 
Aq) = 1+ 2(1+ #)-3— 3(1+ a)’, 
AS) = EM ES LEE) I 
with 
4E? 
A rd 
a 


where Æ, is the C.M. neutrino energy. 

In the low-energy limit the expression for &x vanishes proportionally to &2, and 
it tends to a constant limit for high energies. In (12) the upper sign holds for 
reaction (1), the lower sign for v+n—e-+p. A graph of oj, as given by (12) 


T T Tr T ni T TREAT T T ai T Au T T TETE 
2 49 06) 8.10 2 4 6 80 € 
Fig. 1.— Total cross-section for VY+p—et+n versus neutrino energy in the laboratory as com- 


puted from equation (12). 


is reported in Fig. 1 calculated with u= 3.7: (G,/G)?=1.4 and a?=37.5m2. If the 
three form factors are still taken of the form (10) and (11) but with different a? 
(ai, a, and b? for H,) the high energy limit of the cross-section is 


1 su, 1 
. D 2 be dg - 
(13) co "102 CES — Gb. 
oc. gna 
Such a limit is practically reached already 


for energies k slightly greater than 
the largest of the a’s. 


308 


| 


CROSS SECTIONS OF REACTIONS PRODUCED BY HIGH ENERGY NEUTRINO BEAMS 309 


The high energy cross-section will thus depend very sensitively upon such para- 
meters. One sees from Fig. 1 that the total cross-section for y+p ->e*+n predicted 
from our simple model is ~10-*8 em? for incident neutrinos of laboratory energy 
greater than 1 GeV. The cross section for the process v+n >—e-+p comes out also 


to be of the same magnitude at high energies because of the vanishing in that limit 
of the interference term. 


T T == T 


15 2 ie AST 2 PAT 


Fig. 2.- Total cross-section for v_+ n->e7 +p versus neutrino energy in the laboratory as com- 
puted from equation (12). 


To obtain the total interaction cross-section one should add to such figures the 
value of the cross-section for the analogous processes producing muons in place 
of the electrons. Furthermore, one should also add the cross-sections for reactions 
producing pions and strange particles. We do not know any reliable method to 
evaluate the cross-sections for the processes 


(14) v+QM>MH+ekmx, 


leading to m final pions. When some data on multiple pion production in electron- 
nucleon collisions will be available at comparable energies one will be able to make 
some rough guess on such cross-sections, at least if the non-renormalization hypo- 
thesis holds. Statistical estimates as to the most probable number of pions that 
can be produced in a reaction like (14) are perhaps misleading here because ot the 
presence of the weak leptonic vertex. If, in absence of anything better, one wants 
to adopt statistical models to get an idea on the most probable number of mesons 
that are produced in (14), then one finds at 10 GeV laboratory neutrino energy a 
value ~5 for such a number, which increases to a value ~8 at 20 GeV. We have 
used an interaction radius of the order of the pion wavelength and also taken approx- 
imately into account the indistinguishability for the final pions. Considering the 
different final channels that are available one may hope that the total interaction 
cross-section may turn out to be much larger than the cross-section that we have 
estimated here for v+97 > e+97. Coming back to the question of the high energy 
behaviour of the cross-sections we should perhaps stress the fact that we are not 
making any statement — neither rigorous nor approximate — as to the true high 
energy limit. The energies we are considering are only those practically available 
nowadays. The possible presence of a structure (8) for the weak interactions — for 


(8) T. D. Lee and C. N. YANG: Phys. Rev., 108, 1611 (1957); S. BLUDMAN and A. KLEIN: 
Phys. Rev., 109, 550 (1958); N. CaBrBBo and R. Garro: Phys. Rev., 116, n. 4 (1959). 


309 


310 N. CABIBBO and R. GATTO 


instance if they were mediated by a heavy meson — would also influence the dis- 


cussion of the high energy limit. 
To conclude our discussion we may also remark that other interesting reactions are 


(15) Velos Ven, 
(16) deu 


and similar reactions with v. One would expect that form factors, as arising for 
instance from virtual gammas, do not play an important role in the energy range 
we are considering. If so the cross-sections will increase with the square of the 
C.M. energy. However because of the very low mass of the target in reactions 
like (15) and (16), the C.M. energy is still too small to give a measurable cross- 
section even at (20+ 30) GeV. The cross section for (15), for example, is given, 
already for energies of the order of a few MeV by (°) 


c=8:10-2 & cm?, (& in GeV), 


and is still =10-4 em? at 10 GeV. 
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Handbuch der Physik, Band XVII 
Dielektrika, con 198 figure, 406 pa- 
gine in 8°. Edizioni Springer, 1956. 


Il volume consta di tre lunghi ed 
aggiornati articoli. Nel primo, scritto 
con ammirevole chiarezza da W. FULLER 
Brown jr., l’autore dopo avere definiti 
i principali problemi della teoria dei 
dielettrici e richiamato i concetti di 
elettrostatica, di meccanica quantistica 
e di meccanica statistica di cui intende 
servirsi, affronta lo studio del compor- 
tamento dei dielettrici reali. Si occupa 
dapprima del comportamento dei dielet- 
trici in campi costanti o lentamente 
variabili e la trattazione procede dai 
sistemi più semplici (gas rarefatti) a 
quelli via via più complessi (gas densi, 
fasi condensate e soluzioni). Indi discute 
il comportamento dei dielettrici alle fre- 
quenze ottiche, anche qui cominciando 
dall’esame del comportamento dei gas 
rarefatti per passare in seguito al caso 
più difficile dei gas densi e delle fasi 
condensate. Tratta infine il comporta- 
mento dei dielettrici a frequenze inter- 
medie con particolare riferimento ai pro- 
cessi di rilassamento ed ai meccanismi 
suscettibili di rivelare le proprietà mole- 
colari. La discussione di quest’ultimo 
aspetto del problema viene ripresa e 
sviluppata nell'ultimo capitolo dell’arti- 
colo. Nel penultimo capitolo dell’articolo 
sono esposti i metodi sperimentali di 
misura, con brevi ma sostanziose ed 
aggiornate descrizioni dei circuiti e delle 
tecniche. 
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Nel secondo articolo, opera del Prof. 
Walter Franz, l’autore si occupa del 
meccanismo delle scariche nei dielettrici 
solidi e liquidi e dei processi di ionizza- 
zione interna per urti elettronici. Espone 
dapprima la situazione nel campo degli 
studi sulla scarica « calda », portante a 
zone fuse o gasificate nel dielettrico, ed 
indi passa ad occuparsi più a lungo della 
«scarica elettrica pura» in termini della 
statistica degli elettroni nei materiali 
isolanti. 

Nel terzo articolo l’autore Peter Fors- 
BERGH jr. tratta della piezoelettricità, 
dell’elettrostrizione e della ferroelettri- 
cità ed antiferroelettricità. Ad una intro- 
duzione di carattere generale sui feno- 
meni in oggetto, corredata da una serie 
di bellissime fotografie di domini ferro- 
elettrici e delle relative frontiere, fa 
seguito la trattazione dettagliata della 
piezoelettricità, della ferroelettricità e 
dell’antiferroelettricità. La trattazione è 
completata da un gran numero di figure 
e di grafici. Il penultimo capitolo tratta 
della termodinamica dei cristalli elettro- 
meccanici, ed alla fine di esso viene 
delineata una teoria fenomenologica ca- 
pace di porre in relazione tutte le « osser- 
vabili » macroscopiche, elettriche mec- 
caniche e termiche, che sono le derivate 
parziali dei potenziali termodinamici. 
L'ultimo capitolo abbozza, per quanto 
è a tutt'oggi possibile, una teoria micro- 
scopica statistica della ferroelettricità. 

Leggendo l’opera abbiamo avuto Vim- 
pressione che il primo articolo sia sen- 
y altro al livello dei migliori fra gh altri 
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articoli comparsi nei volumi dell’Hand- 
buch a tutt’oggi pubblicati. Il secondo 
articolo invece ci lascia con la sensazione 
che si sarebbe potuto chiarire meglio 
la portata del problema ed i suoi rap- 
porti coi problemi della fisica dello stato 
solido. Buono il terzo articolo anche se, 
a nostro avviso non vi è sufficientemente 
sviluppata la parte teorica. 


F. GARTA 


Numerische Mathematik. Springer- 
Verlag, Berlin-Gôttingen-Heidel- 
berg. 


Si tratta di una nuova rivista pub- 
blicata dall’Editore Springer, con la col- 
laborazione di un Comitato di redazione 
costituito su base internazionale. Essa 
è destinata a raccogliere lavori originali 
nel campo dell’analisi numerica e dei 
metodi di programmazione per le mo- 
derne macchine calcolatrici. 

Durante il 1959 sono apparsi i 5 fa- 
scicoli che costituiscono il 1° volume di 
320 pagine e contenente 24 lavori; si 
tratta di lavori molto interessanti dovuti 
ai migliori specialisti della materia. Per 
la pubblicazione di lavori nella rivista 
occorre rivolgersi ad uno dei membri 
del sopraddetto Comitato di redazione; 
riteniamo perciò opportuno segnalare che 
il membro italiano è il Prof. A. Guiz- 
zertI dell’Università di Roma. 


A. GHIZZETTI 


W. PAULI — Teoria della Relatività, 
traduzione italiana di Relativitàts- 
theorie, a cura di Paolo Gulma- 
nelli. Editore Paolo Boringhieri, 
Torino, 1958, pag. xVI-327. 


La prima edizione del libro fu stam- 
pata nel 1921 con lo scopo di dare una 
rassegna sui lavori della teoria della rela- 
tività fino allora pubblicati. 


In occasione del cinquantenario della 
teoria della relatività (1955) era proposito 
di Pauli come egli stesso dice nella pre- 
fazione all’ultima edizione, aggiornare il 
libro in modo da estendere il carattere 
di rassegna fino alla nuova data. Tut- 
tavia Egli rinunciò all’idea per la gran 
mole di lavoro che l’aggiornamento 
avrebbe comportato e anche perchè repu- 
tava quasi superfluo, accanto ai numerosi 
trattati sulla Relatività comparsi attorno 
a quell’anno, averne un altro costruito 
mediante la completa trasformazione di 
un libro che già aveva un carattere di 
«documento storico ». 

Sicchè PauLI in definitiva preferì la- 
sciare il vecchio testo nella primitiva 
forma con alcune aggiunte e note che 
contengono gli sviluppi più importanti 
avuti dopo il 1921 e alcune considera- 
zioni sulle teorie unitarie di campi che, 
più che altro, esprimono il parere, invero 
molto scettico, di Pauli sulle possibilità 
di successo di siffatte teorie. 

Il libro si articola su otto capitoli 
alternando e compenetrando fatti spe- 
rimentali e formalismi in modo da avere 
una sintesi molto equilibrata che riesce 
utilissima sia per lo sperimentale che per 
il teorico. Im esso si ha un panorama 
molto vasto che va dalla relatività ri- 
stretta a quella generale, mettendo l’ac- 
cento oltre che sugli aspetti più comu- 
nemente noti di queste teorie (mecca- 
nica, elettromagnetismo, ottica, campo 
gravitazionale e teorie unitarie di Ein- 
stein, Mie e Weyl) anche su altri aspetti 
altrettanto importanti, ma generalmente 
non posti nella giusta luce in altri trat- 
tati sulla relatività; mi riferisco partico- 
larmente ai capitoli dedicati agli sviluppi 
della relatività ristretta concernenti la 
termodinamica e la statistica e alle teorie 
delle particelle elementari cariche. 

Sono da segnalare, per il rigore e la 
chiarezza, alcuni argomenti formali ri- 
guardanti le trasformazioni infinitesimali 
e i teoremi variazionali. Inoltre si pos- 
sono considerare di viva attualità le 
considerazioni sui limiti del concetto 
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classico di campo che si trovano nel 
capitolo sulle particelle elementari. 

Ora che l’Autore è scomparso non 
si può fare a meno di riportare il giudizio 
lodevolissimo dato da EINSTEIN nella 
recensione di questo libro: « Chi legga 
quest'opera densa e dalla solida strut- 
tura non può certamente sospettare che 
l’autore è un giovane di ventun anni. 
Non si sa che cosa ammirare di più: 
se la capacità di cogliere le connessioni 
fra le idee, la sicurezza della deduzione 
matematica, il profondo intuito fisico, 
la chiarezza e la sistematicità dell’espo- 
sizione, la conoscenza della letteratura 
sull’argomento, la completezza della trat- 
tazione o la sicurezza del senso critico ... 
Questo libro andrebbe consultato da 
chiunque faccia lavoro di ricerca nel 
campo della relatività e da chiunque 
desideri un sicuro orientamento su que- 
stioni di principio ). 

Ottima la traduzione e la veste tipo- 
grafica. 

R. Liorra 


E. Persico — Gli atomi e la loro 
energia. Zanichelli, Bologna, 1959. 
Lire 5500, pag. 490. 


Questo libro del Prof. PERSICO viene 
in aiuto a quanti desiderano avvicinarsi 
alla fisica nucleare di oggi, avendo a 
disposizione pochi elementi (o ricordi) 
di fisica generale e qualche elemento di 
calcolo infinitesimale. Perciò finalmente 
abbiamo il libro che consiglieremo agli 
ingegneri ed ai chimici che per necessità 
professionali o per loro diletto ci chie- 
dono spesso un libro del genere. Invece 
esso non ci sembra molto adatto per 
gli studenti, per il carattere un poco 
enciclopedico e di ottimistica facilità con 
cui gli argomenti vengono necessaria 
mente presentati. 

Il libro è diviso in due parti: la prima 
dedicata all’atomistica e la seconda alla 
fisica nucleare. La prima parte ci sembra 
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qualche volta non completamente omo- 
genea come livello degli argomenti trat- 
tati (ad es. la natura del calore e la 
struttura dei liquidi sono trattati ad un 
livello liceale, mentre gli spettri atomici 
sono portati quasi al livello universi- 
tario). La seconda parte è invece più 
omogenea, chiarissima, € veramente pre- 
ziosa per l'attualità degli argomenti e 
la quantità delle informazioni. Infine 
un’appendice sulle fondamentali espe- 
rienze di Fisica Nucleare e soprattutto 
sulle apparecchiature più usate € più 
citate, aggiunge a questo libro il carat- 
tere di vademecum del non-specialista 
in Fisica Nucleare. 

G. CARERI 


Ecole @ été de physique theorique, Les 
Houches, 1958. The many-body 
problem. Ed. Dunod, Paris, pa- 
gine 675, 6900 franchi. 


In questo grosso volume sono rac- 
colte le lezioni tenute a Les Houches 
sulla Fisica Teorica dei sistemi di par- 
ticelle, argomento di grande attualita 
che accomuna campi diversi, come la 
struttura del nucleo, il problema del- 
Lelio liquido e molti aspetti della Fisica 
dei solidi. Si tratta di una quindicina 
di corsi, tutti tenuti da specialisti che 
hanno recentemente contribuito in questi 
stessi argomenti; e tutti i corsi sono 
tenuti ad un livello piuttosto elevato. 
Anzi più che una serie di lezioni si può 
anche parlare di una serie di seminari 
che gli autori hanno tenuto sul loro 
stesso lavoro. 

Un primo gruppo di lezioni, circa la 
metà del volume, si riferisce alla teoria 
della struttura del nucleo ed alle pro- 
prietà della materia nucleare; gli autori 
ne sono: HuGENHOLTZ, BRUECKNER, 
BLocH, De Dominicis, MOTTELSON, 
WEISSKOPF, STRANTINSKY, LIPTIN, BE- 
LIAEV. Un secondo gruppo è invece di- 
retto verso i fluidi costituiti da bosoni 
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e fermioni e gli autori ne sono: Boum, 
PINES, SCHRIEFFER, LYNTON, HUANG. 
Per rendere la raccolta più accessibile 
alle possibilità economiche di molti let- 
tori, l’editore ha anche fornito un’edi- 
zione di prezzo circa metà per le due 
parti separate. 

Non c’è dubbio che si tratta di una 
raccolta di notevole interesse e di grande 
attualità e che lo sforzo della Scuola di 
Les Houches va ancora una volta ricor- 
dato e lodato. Tuttavia questa raccolta, 
considerata come un libro a se stante, 
non è scevra di lacune e risente del 
carattere frammentario dei diversi arti- 
coli. 

G. CARERI 


Advances in Semi-Conductor Science. 
General Editor: H. Brooks, Per- 
gamon Press (pagg. 553, prezzo 
100 s.). 


Il volume contiene i lavori presentati 
al Congresso mondiale sui semicondut- 
tori, tenutosi nello Stato di New York 
dal: 18 al 22 agosto 1958. 

I lavori sono in numero di 132 e 
sono divisi per sessioni ognuna delle 
quali raggruppa ricerche sullo stesso 
argomento riproducendo fedelmente l’an- 
damento del Congresso; ad ogni lavoro 
sono aggiunte le discussioni che ne hanno 
seguito la presentazione. Una sessione 
speciale riassuntiva è dedicata a 5 con- 
ferenze che hanno lo scopo di fare il 
punto della situazione e di sintetizzare 
alcuni dei risultati più salienti presentati. 

Una delle cose che risultano subito 
evidenti è che, mentre il Congresso di 
Amsterdam del 1954 aveva rivelato una 
tecnica completamente nuova per la 
misura delle masse efficaci, questo Con- 
gresso presenta pochissimi lavori sulla 
risonanza ciclotronica, quasi questa tec- 
nica sia stata momentaneamente abban- 
donata. La ragione è probabilmente 
nella limitazione imposta dalla neces- 
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sità di avere per tali esperienze mate- 
riali con un lungo tempo di rilassamento 
per le collisioni dei portatori di carica. 
Il problema attuale sembra quello di 
superare tale limitazione. 

L’enfasi di questo Congresso è posta 
principalmente su argomenti più tra- 
dizionali quali gli effetti dovuti alla 
superficie, l’assorbimento eccitonico e 
studi di mobilità nei semiconduttori più 
noti e nei composti intermetallici. 

I risultati sugli effetti di superficie 
sono ancora di natura piuttosto quali- 
tativa ma permettono di mettere in 
evidenza l’esistenza di livelli elettronici 
accettori alla superficie. Questi sono do- 
vuti all’interruzione della periodicità alla 
superficie e sono stati previsti da TAMM 
nel 1932 sulle basi di un semplice modello 
unidimensionale; essi vengono chiamati 
« stati veloci) a causa del breve tempo 
di rilassamento. Oltre ad essi, altri 
livelli elettronici sono presenti per il 
trattamento che la superficie ha subito 
e sono probabilmente localizzati sulla 
superficie esterna del gas adsorbito: essi 
vengono chiamati « stati lenti). Di qui 
la necessità di distinguere tra superfici 
ideali o « pulite» e superfici « reali ». 
Tra i lavori su superfici pulite sono 
particolarmente notevoli quelli di Fan- 
SWORTH, MISSMAN e HANDLER. 

Due intere sezioni del Congresso sono 
dedicate a lavori su eccitoni. Il concetto 
tradizionale dell’eccitone come uno stato 
legato di elettrone e buca fa pensare 
che ci si dovrebbero aspettare transi- 
zioni del tipo idrogenoide. Una serie di 
misure di Gross e di NIKITINE in alcuni 
composti semiconduttori, quali il Cu,0, 
rivela in effetti tali transizioni, delle 
quali viene anche determinato l’effetto 
Stark e l’effetto Zeeman. Negli alogenuri 
alcalini, tuttavia, si osservano poche 
bande di assorbimento e una loro inter- 
pretazione teorica con il modello di 
FRENKEL perfezionato da WANNIER è 
tuttora incompleto. Osservazioni di assor- 
bimento eccitonico in Germanio e Silicio 
sono presentate da MACFARLANE e col- 
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laboratori e da Haynzs e collaboratori. 
Contributi teorici sugli eccitoni vengono 
presentati e discussi da HAKEN e da 
Muro e collaboratori. 

Tra i molti studi di trasporto vanno 
menzionate le ricerche sulle variazioni 
di mobilità di elettroni in presenza di 
un forte campo elettrico; una rassegna 
completa degli effetti associati a questi 
cosidetti « elettroni caldi » viene presen- 
tata da KOENIG. 

Di estremo interesse sono pure 1 
lavori di Jorré e collaboratori su un 
gran numero di composti intermetallici; 
essi sembrano provare che quando si ha 
a che fare con semiconduttori a bassa 
mobilità, nei quali il libero cammino 
medio degli elettroni e delle buche è 
inferiore al parametro reticolare, i con- 
cetti statici che presiedono alla teoria 
del trasporto nei metalli e nei semicon- 
duttori abituali non possono più essere 
usati. 

Il fondamento teorico per la com- 
prensione generale delle proprietà dei 
semiconduttori rimane la teoria delle 
bande e anche a questo Congresso due 
sezioni sono dedicate a calcoli di questo 
tipo. In alcuni lavori appare lo sforzo 
di svolgere calcoli sempre più precisi, 
e in altri lo sforzo di trovare schemi di 
interpolazione 0 approssimazioni che 
consentano di ragionare anche su cristalli 
complessi. Nello sviluppare la teoria 
delle bande si trascurano le correlazioni 
tra gli elettroni dovute alla loro repul- 
sione coulombiana e ci si pud chiedere 
quale sia la validita dei concetti che in 
tale modo si ottengono. Particolare 
importanza in questo senso hanno alcuni 
lavori di Koun e altri che descrivono 
elettroni e buche come eccitazioni di 
funzioni d’onda a molte particelle; utiliz- 
zando il metodo dei « linked clusters », 
Koun ha giustificato il concetto di massa 
efficace dal punto di vista del problema 
dei molti corpi come pure quello di co- 
stante dielettrica. 

Molto rimarebbe ancora da dire ma 
non è purtroppo possibile commentare in 
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modo completo un volume di questa 
natura. Mi limito a raccomandarlo all’at- 
tenzione di tutti i Fisici. 


E. BASSANI 


Vacancies and Other Point Defects in 
Metals and Alloys. Institute of 
Metals Monograph no. 23 (The 
Institute of Metals, London 1958), 
pag. 248, 40 s. 


Il libro raccoglie i testi di sei lavori, 
aventi sostanzialmente carattere di ras- 
segna, che furono presentati al simposio 
dell’ Institute of Metals tenuto a Harwell 
nel dicembre 1957. Cinque articoli sono 
dedicati ai difetti puntiformi, in parti- 
colare posti reticolari vacanti, in metalli 
e leghe, mentre un articolo riguarda i 
difetti puntiformi nei solidi ionici, con 
particolare riferimento agli alogenuri 
alcalini. 

I primi due articoli, dovuti a Cor- 
TRELL ed a Broom e Ham, discutono 
in dettaglio gli effetti esercitati sulle 
proprieta fisiche dei metalli, a tempera- 
tura di autodiffusione, da difetti punti- 
introdotti in soprassaturazione 
mediante irraggiamento o tempra o 
lavorazione a freddo o (in composti 
intermetallici) producendo deviazioni con- 
trollate dalla composizione stechiome- 
trica. In particolare nell’articolo di CoT- 
TRELL vengono trattati Vindurimento da 
irraggiamento, da lavorazione a freddo 
e da fatica, la determinazione delle 
energie di formazione e di moto di posti 
vacanti nei metalli mediante tecniche 
di tempra, ed i vari tipi di interazione 
fra difetti puntiformi e dislocazioni; 
mentre il lavoro di Broom e Ham con- 
sidera gli effetti dei difetti puntiformi 
su energia interna, densità e parametro 
reticolare ai raggi X, proprietà elettriche, 
termiche e magnetiche, e risonanza 
magnetica nucleare. L'articolo di LOMER, 
e quello di McLean, studiano rispetti- 
vamente il ruolo dei difetti puntiformi 
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nella diffusione e la loro influenza sulle 
proprietà meccaniche, con particolare 
enfasi sull’ascensione delle dislocazioni. 
L'articolo di WILLIAMS e HAYFIELD 
riporta interessanti studi dello stato di 
disordine presso la superficie dei metalli, 
e di fenomeni di notevole importanza 
tecnica quali ossidazione e corrosione: un 
soggetto di recente interesse, e rara- 
mente trattato. Infine il lavoro di PRATT 
costituisce una valida, seppur in alcuni 
punti frettolosa, rassegna della vasta 
informazione raccolta sull’influenza dei 
difetti puntiformi sulle proprietà mec- 
caniche dei cristalli ionici, e presenta le 
ricerche recenti sullo sforzo critico di 
deformazione plastica in alogenuri alca- 
lini soggetti a vari trattamenti: appare 
chiaro dalla lettura dell’articolo che 
accanto allo studio delle proprietà pla- 
stiche sarebbe utile riprendere lo studio 
dell’attrito interno in questi solidi, poichè 
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in essi l’interazione fra difetti punti- 
formi e dislocazioni presenta notevoli 
proprietà caratteristiche, assenti nei me- 
talli. 

L'assenza di un capitolo introduttivo 
rende il libro adatto soltanto a ricer- 
catori già introdotti nel campo dei 
difetti nei solidi. Per questi il libro, 
opera di specialisti nel campo, apparirà 
di notevole interesse, e sarà utile sia 
come rassegna critica di risultati spe- 
rimentali e teorici per lo più dispersi in 
periodici specializzati, sia per la vastis- 
sima bibliografia. Sebbene la rassegna 
ricopra l’informazione raccolta nel pe- 
riodo pre-1957, informazione su ricerche 
di più recente pubblicazione è contenuta 
nelle trentacinque pagine finali, ove sono 
presentati i testi delle discussioni che 
hanno accompagnato la presentazione 
dei lavori. 
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